Next Article in Journal
Inversion of Forest Biomass Based on Multi-Source Remote Sensing Images
Next Article in Special Issue
K-Space Approach in Optical Coherence Tomography: Rigorous Digital Transformation of Arbitrary-Shape Beams, Aberration Elimination and Super-Refocusing beyond Conventional Phase Correction Procedures
Previous Article in Journal
Time-Frequency Multi-Domain 1D Convolutional Neural Network with Channel-Spatial Attention for Noise-Robust Bearing Fault Diagnosis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Tail Artifact Removal via Transmittance Effect Subtraction in Optical Coherence Tail Artifact Images

1
Faculty of Mathematics and Physics, University of Ljubljana, 1000 Ljubljana, Slovenia
2
Jozef Stefan Institute, 1000 Ljubljana, Slovenia
*
Author to whom correspondence should be addressed.
Submission received: 6 October 2023 / Revised: 14 November 2023 / Accepted: 20 November 2023 / Published: 21 November 2023

Abstract

:
Optical Coherence Tomography Angiography (OCTA) has revolutionized non-invasive, high-resolution imaging of blood vessels. However, the challenge of tail artifacts in OCTA images persists. In response, we present the Tail Artifact Removal via Transmittance Effect Subtraction (TAR-TES) algorithm that effectively mitigates these artifacts. Through a simple physics-based model, the TAR-TES accounts for variations in transmittance within the shallow layers with the vasculature, resulting in the removal of tail artifacts in deeper layers after the vessel. Comparative evaluations with alternative correction methods demonstrate that TAR-TES excels in eliminating these artifacts while preserving the essential integrity of vasculature images. Crucially, the success of the TAR-TES is closely linked to the precise adjustment of a weight constant, underlining the significance of individual dataset parameter optimization. In conclusion, TAR-TES emerges as a powerful tool for enhancing OCTA image quality and reliability in both clinical and research settings, promising to reshape the way we visualize and analyze intricate vascular networks within biological tissues. Further validation across diverse datasets is essential to unlock the full potential of this physics-based solution.

1. Introduction

Optical Coherence Tomography (OCT) has emerged as a non-invasive imaging tool, offering high-resolution views of biological tissues and materials [1,2]. Operating on low-coherence interferometry, OCT provides us with precise cross-sectional images of internal structures.
OCT imaging goes beyond static imaging by detecting motion within samples, enabling real-time imaging of dynamic processes. One such OCT technique is Optical Coherence Tomography Angiography (OCTA), a potent technique for visualizing blood vessels and microvascular networks in living tissues [3,4,5,6]. It relies on dynamic contrast, where the movement of red blood cells (RBCs) in vessels causes fluctuations in OCT signals, distinguishing blood vessels from surrounding tissues [5,7,8,9]. OCTA is widely used in ophthalmology, neuroscience, cancer research, and dermatology [8,9,10,11,12].
OCTA faces a significant challenge—occurrence of imaging artifacts known as projection artifacts or tail artifacts. These artifacts occur for various reasons, including variation of transmitted incident light to deeper layers due to the variable transmittance in the shallow layers with the vasculature [13,14] and elongated light paths from multiple scattering of photons interacting with flowing RBCs [15,16]. Variation of transmitted incident light to deeper layers due to the variable transmittance in the shallow layers causes OCT signal variation even in the absence of reflectivity variation, generating vasculature-like signals beneath the vasculature. Elongated light paths from multiple scattering of photons increase the delay of reflected light from the vessels in shallow layers and attribute that signal to deep layers.
Several strategies exist to mitigate those imaging artifacts and often involve image processing techniques and hardware solutions. An example of a hardware solution is a high numerical aperture objective that minimizes multiple scattering effects. However, the high numerical aperture imaging also requires dynamic focusing to reconstruct the OCTA, which may not always be a trivial task and therefore limits widespread adoption [17]. In addition, it only removes the effects of elongated light paths from multiple scattering of photons. Early examples of image processing techniques are Step-down Exponential filtering [18], slab-subtraction [19,20], and projection-resolved algorithm [21]. The slab-subtraction and projection-resolved algorithm are quite popular, especially in ophthalmological imaging. Our experiences with these algorithms on skin OCTA images are less promising. The Step-down Exponential filtering is also partially effective but may cause signal loss, especially in deep vascular regions. A more recent Mean-Subtraction method removes tail artifacts by subtracting the weighted mean value of A-scan but may disrupt deep flows and hinder quantitative blood flow analysis [16]. Deep learning-based OCTA can suppress tail artifacts effectively but demands computational resources and extensive annotated datasets, making manual labeling labor intensive and prone to errors [22].
In this work, we introduce an innovative physics-based OCTA tail artifact removal algorithm. The algorithm relies on a simple OCTA physical model for the main source of the tail artifact—the variation of transmitted incident light to deeper layers due to the variable transmittance in the shallow layers with the vasculature. We demonstrate our algorithm’s effectiveness on hemangioma vasculature OCTA images. As a validated and simple solution backed with a physical model, this algorithm provides a new solution for tail artifact suppression and improved blood vessel network visualization.

2. Materials and Methods

2.1. Tail Artifact Removal via Transmittance Effect Subtraction (TAR-TES) Algorithm

The new OCTA tail removal algorithm, the Tail Artifact Removal via Transmittance Effect Subtraction (TAR-TES), calculates the corrected OCTA image by subtracting the contribution of transmittance variation from the OCTA image. The structural image intensity (also known as OCT image intensity) in the n-th layer, Sn, is proportional to the fraction of incident light transmitted to the n-th layer, Tn, and the reflectivity of the n-th layer, Rn, with w 1 1 being the proportionality factor. We write proportionality factor as the inverse of w1 to ensure that the final equation has it as a linearly proportional factor.
S n = T n · R n w 1
The angiographic image intensity (also known as OCTA image intensity) in the n-th layer, An, is proportional to the variation in the structural image.
A n = δ T n · R n w 1 = δ T n · R n + T n · δ R n w 1
We assume that the corrected OCTA image intensity is proportional to the variation in reflectivity. From this, we can calculate the corrected OCTA by subtracting the variable transmittance effect from the uncorrected angiographic image.
A n c o r r = T n · δ R n w 1 = A n δ T n · R n w 1
The variation in the transmitted light to the n-th layer δTn depends on the variation in reflected light in previous layers, i.e., on the corrected angiographic image. It is calculated as the square root of the sum of squared variations in reflectivity from layers 1 to n − 1, assuming that the variations in different layers are uncorrelated.
δ T n = δ 1 i = 1 n 1 R i = i = 1 n 1 δ R i 2 = w 1 · i = 1 n 1 A i c o r r 2 T i 2
We can insert Equations (1) and (4) into Equation (3) to substitute the variation in the transmitted light to the n-th layer δTn with the corrected OCTA from the previous layers and the reflectivity Rn with the structural image Sn. Thus, we obtain the formula for the correction of the n-th layer in the OCTA image.
A n c o r r = A n w 1 · S n T n · i = 1 n 1 A i c o r r 2 T i 2
We further assume that almost all incident light is transmitted light to the n-th layer, so we set T i = 1 . In that way, we obtain the simplified formula for the corrected OCTA image. It is calculated by subtracting scaled structural image Sn, multiplied by the square root of the sum of squared corrected angiographic images from layers 1 to n − 1 from An.
A n c o r r = A n w 1 · S n · i = 1 n 1 A i c o r r 2
This new algorithm provides a novel approach for removing tail artifacts in OCTA images by accounting for the contribution of variable transmittance to the deep layers in OCTA images. It has one free parameter, w1, which is adjusted to control the strength of the OCTA tail-removal correction.

2.2. Vasculature Model

In this study, we imaged the hand’s birthmark, namely hemangioma, a vascular anomaly characterized by an abnormal overgrowth of blood vessels. Hemangiomas can generally occur in various tissues, including the skin, and are often characterized by a complex network of blood vessels. Hemangiomas are generally harmless, and they offer a good model for blood vasculature [23]. OCTA offers a non-invasive means of imaging hemangiomas, providing high-resolution, depth-resolved images of the vasculature.
The particular hemangioma imaged in this study (Figure 1) served as a representative sample for assessing the efficacy of the Tail Artifact Removal via Transmittance Effect Subtraction (TAR-TES) algorithm in correcting tail artifacts in OCTA images. The complex and intricate vascular structure within the hemangioma allowed for a comprehensive evaluation of the algorithm’s performance in a real-world clinical context. Although analysis of only one vascular image is presented here to demonstrate the performance of the TAR-TES algorithm, similar results were obtained for other vascular lesion OCTA images.

2.3. Speckle-Variance OCT Angiography

Optical Coherence Tomography Angiography (OCTA) was used in this study. The imaging was conducted using a commercial Thorlabs VEG220C1 swept-source OCT system, which featured a central wavelength of 1300 nm. The selection of this particular wavelength offers distinct advantages, notably enhanced tissue penetration and improved contrast for visualizing vasculature structures, making it well-suited for in-depth vascular imaging.
To capture angiographic information, a speckle variance OCTA technique was utilized. The OCTA imaging was acquired with a series of five image repetitions. The image was reconstructed with the software that was provided with the OCT system (ThorImage OCT Version 5.5.5.0). The imaging voxel size was set at 2.3 μm × 2.2 μm × 8.1 μm, and the image size was 357 × 373 × 960. OCT and OCTA images were filtered with a Gaussian filter with sigma equals one voxel size. The noise level was reduced by subtracting the average imaging signal deep inside the tissue (depth 4–7 mm), where only homogeneous noise is present in the image.

2.4. Tail Artifact Correction in OCTA Image

Tail artifact correction via the TAR-TES algorithm was performed by applying Equation (5) for each consecutive layer. Weight w1 was adjusted manually; it was increased until tail artifacts visually disappeared (default w1 value). We also tested 2-times lower w1 value and 2.5-times higher w1 value.
The TAR-TES algorithm was compared to two most similar competitive algorithms: the Mean-Subtraction algorithm [16] and the Step-down Exponential filtering method [18]. They were implemented using Equations (6) and (7). One free parameter in those algorithms, w2 or w3, was adjusted manually to achieve roughly the same level of signal decrease as with the default w1 value.
In comparison to the TAR-TES algorithm, the Mean-Subtraction algorithm reduces tail artifacts in OCTA images by subtracting the weighted mean of the depth profile.
A c o r r = A w 2 N i = 1 N A i
Like the TAR-TES algorithm, it has one free parameter, w2, which is adjusted to control the strength of the OCTA tail-removal correction. Unlike the TER-TES algorithm, the Mean-Subtraction algorithm reduces tail artifacts by subtracting the same value, the mean of the uncorrected depth profile, from all layers. In contrast, the TAR-TES algorithm subtracts the contribution of transmittance variation to the angiographic image separately for each consecutive layer and uses corrected images from previous layers.
An older Step-down Exponential filtering method, which is another well-established method for OCTA tail removal, attenuates the OCTA signal of each consecutive layer by a factor proportional to the exponential sum of values from the corrected images from previous layers.
A n c o r r = A n × e w 3 × i = 1 n 1 A i c o r r
Like the TAR-TES algorithm, the Step-down Exponential filtering method performs correction for each consecutive layer and uses corrected images from previous layers. Similarly, it has one free parameter, w3, which is adjusted to control the strength of the OCTA tail-removal correction. In contrast to the TAR-TES algorithm, the Step-down Exponential filtering method performs correction by scaling each layer, not by subtraction.

3. Results

The structural OCT image S of the hemangioma in all three views is shown in Figure 2. Each view is marked with a cross symbolizing the position of the section corresponding to the other two views. These three sections collectively define a single point within the vessel that is not discernible in the OCT image.
The uncorrected OCTA image exhibits substantial tails beneath the larger superficial vessels and lingering shadow artifacts in deep layers. This is a known OCTA artifact, and a proper vasculature image would show a sharp decrease in the OCTA signal after the vessels—as we see a sharp increase in the signal at the front edge of the vessel. A schematic representation of the OCTA tail is shown in Figure 3.
The TAR-TES algorithm performance for tail artifact removal was examined with weight constants (w1) set at 0.01, 0.02, and 0.05 (Figure 4 and Figure 5). The correction with an adequately adjusted weight constant for our proposed TAR-TES algorithm must effectively suppress tail artifacts while keeping the vasculature image.
The removal of tail artifacts for weight w1 = 0.02 is shown in Figure 4. The vessel marked with the cross is clearly visible on the en face view, and the B-scan cross-sectional and C-scan cross-sectional views show the expected shape of the vessel—a shape that resembles the circle. Tail artifacts are slightly visible but far less than on the OCTA with no correction.
The removal of tail artifacts for different weight conditions (w1 = 0.01 and w1 = 0.05) is shown in Figure 5. At the weight value of w1 = 0.01, we observe partial tail artifact removal, indicating that the weight is insufficient. At the weight value of w1 = 0.05, we observe complete tail artifact removal and missing OCTA signal in some parts of the image, which indicates that the image is over-corrected. Also, the cross-section of the vessel is a flattened circle, which is another indicator that the weight might be too high.
To evaluate the TAR-TES algorithm performance, we compared it to the Mean-Subtraction algorithm [16] and the Step-down Exponential filtering method [18]. Weights for both correction methods were adjusted manually; they were increased until the tail artifacts visually disappeared. As with the TAR-TES algorithm, we also did the correction with 2-times smaller weight and with the 2.5-times larger weight. The results are shown in Figure 6, Figure 7 and Figure 8. Due to the big difference in results for the Mean-Subtraction algorithm with the primarily selected weight and the 2.5-times larger weight, we added intermediate weights.
In essence, the Mean-Subtraction algorithm subtracts a baseline value for each OCTA A-scan. The subtracted baseline is different for each OCTA A-scan, so the method is similar to the window setting. In some parts of the corrected OCTA image, the signal is missing because the original OCTA was lower than the subtracted value. When the weight is below the optimal value, the tail artifact is less effectively corrected—seen when comparing the first and the second row. When these images are compared (the first row with the w2 = 3 that is considered as an optimal and the second row with the w2 = 1.5), the differences appear to be small. Therefore, we consider that the range of optimal weight values is broad and manual adjustment is sufficient. If the weight is too high (third row in Figure 6), almost the whole image disappears, which is easy to notice. Due to this substantial difference, the correction with intermediate correction weight values is presented in Figure 7. As we see, the tail artifact is disappearing when the correction weight is increasing, and also the amount of vasculature visible on en face image is reduced.
The Step-down Exponential filtering method carries out a filtering along the OCTA A-scan direction. It reduces the OCTA tail artifacts, but the reduction is not complete. For any weight value, the image value stays above zero, which is expected from Equation (8). The difference is that all three images are hardly noticeable, so we consider that the range of optimal weight values is broad and manual adjustment is sufficient.
Additional insights into the TAR-TES algorithm performance, as compared to the Mean-Subtraction algorithm and the Step-down Exponential filtering method, can be gathered from Figure 9, which compares en face images at various depths. We can see that the vessel is clearly visible at depth of 1.524 mm on all images. Figure 4, Figure 5, Figure 6 and Figure 7 show en face image at depth of 1.540 mm, and we can see the selected vessel on all images very well. Going to a depth of 1.573 mm, we can see the differences between the correction algorithms; TAR-TES corrected OCTA already has relatively low signal for the vessel, whereas the original OCTA, the Mean-Subtraction algorithm and the Step-down Exponential filtering method have increased signal in comparison to the depth of 1.524 mm. At the depth of 1.661 mm, the vessel on the TAR-TES corrected OCTA has already vanished, whereas the original OCTA, the Mean-Subtraction algorithm and the Step-down Exponential filtering method still have clearly visible vessel.
Another comparison of different OCTA correction methods is presented in Figure 10, which provides A-scan profiles of the original and corrected OCTA images for the point defined by the cross-center point in Figure 4, Figure 5, Figure 6, Figure 7, Figure 8 and Figure 9. As we see, all correction methods decrease the OCTA signal and change the shape of the depth profile. The tail artifact is seen as the slowly falling curve after the peak, while the correct shape should be a symmetric peak, i.e., the profile of the vessel (Figure 3). When the signal is normalized to its maximal value, we can compare the peaks and tail artifacts for different tail-correction algorithms. The Mean-Subtraction algorithm basically cuts the lower part of the curve while the shape of the upper part of the curve remains unchanged. The Step-down Exponential filtering method changes the shape of the curve. It reduces the tail artifact, but less than the other two methods. The proposed TAR-TES method is the most effective in reducing the tail artifact, which can be seen from the depth profile. For example, the OCTA signal at 0.2 mm beneath the peak falls to 20% of its maximal value in the case of the TAR-TES method, whereas the Mean-Subtraction and the Step-down Exponential filtering method reduce the OCTA signal to about 40% of its maximal value at this point. The uncorrected OCTA signal is at 60% of its maximal value. Interestingly, the TAR-TES method also moves the peak forward for about 0.05 mm. It should be noted that the diameter of the vessel is about 0.1 mm.

4. Discussion

The presented manuscript introduces the TAR-TES algorithm for correcting tail artifacts in Optical Coherence Tomography Angiography (OCTA) images and evaluates its performance compared to alternative correction methods. The results demonstrate the effectiveness of the TAR-TES algorithm in addressing tail artifacts and improving the quality of OCTA images.
The TAR-TES algorithm, with an appropriately adjusted weight constant w1, showed promising results in artifact removal. Generally, the weight constant w1 determines the relation between the reflectivity of the tissue, light transmittance to that point and OCT signal. This relation also implicitly includes the efficiency of the detection system. Altogether, it would be very difficult to determine the weight constant w1 from the basic principles. It remains the free parameter of the correction algorithm, which should be adjusted for the maximal efficiency of the algorithm. When adjusted properly, the TAR-TES algorithm strikingly reduced the tail artifacts while preserving the vasculature image’s integrity. The balance achieved in artifact removal is essential to ensure that the corrected OCTA images remain true to the underlying vasculature morphology.
The adjustment of weight constant w1 can be viewed as a regularization problem, where the appropriate regularization aggressiveness is necessary to achieve optimal results. An investigation into different weight constants shed light on the significance of parameter selection in the TAR-TES algorithm. A weight constant of w1 = 0.01 displayed partial tail artifact removal, suggesting that a lower weight is insufficient for comprehensive correction. On the other hand, w1 = 0.05 exhibited complete tail artifact removal but led to signal loss and alterations in vessel shape. This underscores the importance of fine-tuning the weight constant based on the specific characteristics of OCTA images and the desired level of artifact reduction. It is essential to note that the parameters we used are not the only viable choices and could be optimized in different ways. Additionally, different samples may require different parameter settings. Different scattering properties imply different magnitude of the transmission to deeper layers and different magnitude of reflected light, which implies that different parameter w1 is optimal for each particular setting.
Comparative analysis with two alternative correction methods, namely the Mean-Subtraction algorithm and the Step-down Exponential filtering method, highlighted the advantages of the TAR-TES. The Mean-Subtraction algorithm, which subtracts a varying baseline value for each A-scan, reduces the OCTA signal after the vessels only because the baseline value is subtracted from the signal. The Step-down Exponential filtering method changes the shape of the A-scan profile of the OCTA signal. However, it did not completely eliminate tail artifacts. In contrast, the TAR-TES algorithm almost completely eliminated tail artifacts. On the other hand, the selected alternative methods may be less sensitive to the selection of the weight constant w2 and w3, with the exception of excessive correction weight in the case of the Mean-Subtraction algorithm that is easy to recognize. Due to this drastic reduction in the amount of vasculature in case of excessively aggressive correction with the Mean-Subtraction algorithm, we presented intermediate weights for the correction with this algorithm. In this case, the amount of vasculature is gradually decreasing, while the tail artifacts are also reduced. If we pick the specific vessel that we want to visualize, we can tune the Mean-Subtraction algorithm for this specific vessel so that only its tail artifact is removed while the vessel persists. However, such a procedure is almost the same as window setting in visualization. Therefore, its usefulness for the visualization is limited.
The A-scan profile analysis further clarifies how the TAR-TES algorithm reduces tail artifacts and enhances the depth-profile shape. The plot of the absolute signal values (Figure 10) shows that each correction method considerably reduces the OCTA signal strength. By normalizing the signal to its maximum value, it became evident that the TAR-TES algorithm is most successful in transforming the original OCTA signal into the expected symmetric peak.
The TAR-TES algorithm, while promising in addressing tail artifacts in Optical Coherence Tomography Angiography (OCTA) images, has several limitations. Its performance is sensitive to the selection of the weight constant w1, which can vary across datasets and imaging conditions. A high weight constant can lead to significant OCTA signal loss, potentially resulting in missing vasculature information. The algorithm may also alter the shape of vessels, which could impact the accuracy of quantitative vascular analysis. The performance of the algorithm can vary across different OCTA datasets due to the specific imaging conditions, patient populations, and imaging equipment used.
Depending on the clinical situation, motion artifacts can be a significant problem of the speckle-variance OCTA. Motion during the consecutive OCT scans also produce speckle patterns, similar as the blood flow, and therefore false flow image. One option to reduce this false flow signal is to perform complex-signal-based OCTA [7] or high-filtering of complex-signal-based OCTA [24]. Those techniques still suffer from the tail artifacts and the physics basis is similar as in the case of speckle-variance OCTA. Therefore, the ideas behind the TAR-TES algorithm may also be used for the tail artifact correction in these angiographic images, while the exact formulation may not be exactly as presented in this manuscript.
Possible future improvement in the TAR-TES algorithm for using it in the speckle-variance OCTA lies in the optimization of the parameter w1. An obvious improvement is the optimization of the algorithm on a broader range of datasets. Another possible direction is through the integration of artificial intelligence. The algorithm could be trained on simulated data or data obtained from phantom studies, enabling it to adapt and learn optimal parameter settings for the weight constant w1, tailored to specific imaging conditions and datasets. The incorporation of artificial intelligence would not only enhance the algorithm’s robustness and applicability across diverse scenarios but also make it more efficient and user-friendly. The third direction of future work is a modification of the algorithm by substituting the weight constant w1 with the function of depth. This improvement would correct the simplification that we do when going from Equation (4) to Equation (5), where we assume that the transmission factor Tn is equal to 1.
While these incremental improvements in the TAR-TES algorithm could provide superb results, compared to existing algorithms, the idea behind the TAR-TES algorithm (especially the physical model behind) may also serve as a ground for new, better correction algorithms. The problem of the TAR-TES and Mean Subtraction algorithm is that they can yield non-physical negative amplitudes, which are subsequently set to zero (hence grey areas for those two algorithms in Figure 5 and Figure 6). Interestingly, the Step-down Exponential filtering does not suffer from this problem, which is obvious from Equation (8)—correction is carried out using an always-positive scaling factor instead of subtraction. If the correction is small, the correction for Step-down Exponential filtering can be derived as
A n c o r r = A n × e w 3 × i = 1 n 1 A i c o r r A n 1 w 3 × i = 1 n 1 A i c o r r
While this derivation is not valid for our case as the argument in the exponent for the location of vessels is close to 1, the form of Equation (8) is very similar to the TAR-TES equation (Equation (6)). One difference is the summation of corrections for shallower layers in squares and the other difference is scaling the correction with structural image Sn. The latter is not a new idea as it is an essential part in the projection-resolved correction method [21].
The TAR-TES algorithm has the potential for widespread adoption and great impact on future clinical research and diagnostic applications of the OCTA. It is poised to enable true 3D vascular imaging by eliminating tail artifacts in OCTA images. This would enable true 3D reconstruction of vasculature topology and enhance the detection of deep-seated vessels beneath shallower vessels, providing a more accurate and comprehensive representation of complex vascular networks. This advancement holds significant promise in fields such as dermatology, ophthalmology, or oncology.

5. Conclusions

In conclusion, the TAR-TES algorithm emerges as a robust and valuable tool for correcting tail artifacts in OCTA images. Unlike the alternative correction methods, this algorithm is physics based. It provides enhanced image quality without sacrificing the accuracy of vessel morphology representation. This study underscores the importance of parameter optimization in the application of the TAR-TES algorithm and highlights its potential to improve the reliability of OCTA imaging in various clinical and research settings. Further research and validation across diverse datasets are warranted to fully harness the algorithm’s capabilities and ensure its broader applicability in the field of OCTA imaging.

Author Contributions

Conceptualization, U.S. and M.M.; Methodology, U.S.; Software, U.S.; Formal Analysis, U.S.; Investigation, U.S.; Resources, M.M. and U.S.; Data Curation, U.S.; Writing—Original Draft Preparation, M.M. and U.S.; Writing—Review and Editing, M.M. and U.S.; Visualization, U.S.; Supervision, M.M.; Project Administration, M.M.; Funding Acquisition, M.M. and U.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the Slovenian Research Agency through the research program P1-0389 and research projects J7-3150, J7-2600, and J3-3083.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board of the Faculty of mathematics and physics, University of Ljubljana. Given that OCTA imaging of hemangioma of healthy volunteer pose no risk and minimal discomfort to the volunteer, no separate ethics review was required, in accordance with institutional policies.

Informed Consent Statement

Informed consent was obtained from the subject involved in the study.

Data Availability Statement

The data that support the findings of this study are available upon reasonable request from the authors.

Acknowledgments

The authors acknowledge the support of Jost Stergar and Crt Keber in preparation of the OCTA image.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gerd, H.; Michael, W.L. “Coherence radar” and “Spectral radar”—new tools for dermatological diagnosis. J. Biomed. Opt. 1998, 3, 21–31. [Google Scholar]
  2. Huang, D.; Swanson, E.A.; Lin, C.P.; Schuman, J.S.; Stinson, W.G.; Chang, W.; Hee, M.R.; Flotte, T.; Gregory, K.; Puliafito, C.A.; et al. Optical Coherence Tomography. Science 1991, 254, 1178–1181. [Google Scholar] [CrossRef] [PubMed]
  3. Enfield, J.; Jonathan, E.; Leahy, M. In vivo imaging of the microcirculation of the volar forearm using correlation mapping optical coherence tomography (cmOCT). Biomed. Opt. Express 2011, 2, 1184–1193. [Google Scholar] [CrossRef] [PubMed]
  4. Jia, Y.; Tan, O.; Tokayer, J.; Potsaid, B.; Wang, Y.; Liu, J.J.; Kraus, M.F.; Subhash, H.; Fujimoto, J.G.; Hornegger, J.; et al. Split-spectrum amplitude-decorrelation angiography with optical coherence tomography. Opt. Express 2012, 20, 4710–4725. [Google Scholar] [CrossRef] [PubMed]
  5. Mariampillai, A.; Standish, B.A.; Moriyama, E.H.; Khurana, M.; Munce, N.R.; Leung, M.K.K.; Jiang, J.; Cable, A.; Wilson, B.C.; Vitkin, I.A.; et al. Speckle variance detection of microvasculature using swept-source optical coherence tomography. Opt. Lett. 2008, 33, 1530–1532. [Google Scholar] [CrossRef] [PubMed]
  6. Zhang, Q.; Lee, C.S.; Chao, J.; Chen, C.-L.; Zhang, T.; Sharma, U.; Zhang, A.; Liu, J.; Rezaei, K.; Pepple, K.L.; et al. Wide-field optical coherence tomography based microangiography for retinal imaging. Sci. Rep. 2016, 6, 22017. [Google Scholar] [CrossRef] [PubMed]
  7. Chen, C.-L.; Wang, R.K. Optical coherence tomography based angiography [Invited]. Biomed. Opt. Express 2017, 8, 1056–1082. [Google Scholar] [CrossRef]
  8. Chalam, K.V.; Sambhav, K. Optical Coherence Tomography Angiography in Retinal Diseases. J. Ophthalmic Vis. Res. 2016, 11, 84–92. [Google Scholar] [CrossRef]
  9. Meiburger, K.M.; Chen, Z.; Sinz, C.; Hoover, E.; Minneman, M.; Ensher, J.; Kittler, H.; Leitgeb, R.A.; Drexler, W.; Liu, M. Automatic skin lesion area determination of basal cell carcinoma using optical coherence tomography angiography and a skeletonization approach: Preliminary results. J. Biophotonics 2019, 12, e201900131. [Google Scholar] [CrossRef]
  10. Baran, U.; Choi, W.J.; Li, Y.; Wang, R.K. Tail artifact removal in OCT angiography images of rodent cortex. J. Biophotonics 2017, 10, 1421–1429. [Google Scholar] [CrossRef]
  11. Choi, W.J.; Li, Y.; Wang, R.K. Monitoring Acute Stroke Progression: Multi-Parametric OCT Imaging of Cortical Perfusion, Flow, and Tissue Scattering in a Mouse Model of Permanent Focal Ischemia. IEEE Trans. Med. Imaging 2019, 38, 1427–1437. [Google Scholar] [CrossRef] [PubMed]
  12. Batioglu, F.; Yanik, O.; Demirel, S.; Ozmert, E. Clinical Use of Optical Coherence Tomography Angiography in Retinal Diseases. Diagnostics 2023, 13, 1820. [Google Scholar] [CrossRef] [PubMed]
  13. Spaide, R.F.; Fujimoto, J.G.; Waheed, N.K. Image Artifacts in Optical Coherence Tomography Angiography. Retina 2015, 35, 2163–2180. [Google Scholar] [CrossRef] [PubMed]
  14. Zhang, W.; He, B.; Wu, Y.; Tao, Y.; Zhu, F.; Cai, W.; Liu, N.; Zhao, Q.; Xue, P. Tail artifacts removal of three-dimensional optical coherence tomography angiography with common parts extraction method. J. Biophotonics 2022, 15, e202200155. [Google Scholar] [CrossRef]
  15. Yuntao, L.; Jianbo, T. Blood vessel tail artifacts suppression in optical coherence tomography angiography. Neurophotonics 2022, 9, 021906. [Google Scholar]
  16. Choi, W.J.; Paulson, B.; Yu, S.; Wang, R.K.; Kim, J.K. Mean-Subtraction Method for De-shadowing of Tail Artifacts in Cerebral OCTA Images: A Proof of Concept. Materials 2020, 13, 2024. [Google Scholar] [CrossRef] [PubMed]
  17. Conor, L.; Harsha, R.; Marcel, B.; Vivek, J.S. Imaging and graphing of cortical vasculature using dynamically focused optical coherence microscopy angiography. J. Biomed. Opt. 2016, 21, 020502. [Google Scholar]
  18. Mahmud, M.S.; Cadotte, D.W.; Vuong, B.; Sun, C.; Luk, T.W.; Mariampillai, A.; Yang, V.X. Review of speckle and phase variance optical coherence tomography to visualize microvascular networks. J. Biomed. Opt. 2013, 18, 50901. [Google Scholar] [CrossRef] [PubMed]
  19. Liu, L.; Gao, S.S.; Bailey, S.T.; Huang, D.; Li, D.; Jia, Y. Automated choroidal neovascularization detection algorithm for optical coherence tomography angiography. Biomed. Opt. Express 2015, 6, 3564–3576. [Google Scholar] [CrossRef]
  20. Zhang, A.; Zhang, Q.; Wang, R.K. Minimizing projection artifacts for accurate presentation of choroidal neovascularization in OCT micro-angiography. Biomed. Opt. Express 2015, 6, 4130–4143. [Google Scholar] [CrossRef]
  21. Zhang, M.; Hwang, T.S.; Campbell, J.P.; Bailey, S.T.; Wilson, D.J.; Huang, D.; Jia, Y. Projection-resolved optical coherence tomographic angiography. Biomed. Opt. Express 2016, 7, 816–828. [Google Scholar] [CrossRef] [PubMed]
  22. Stefan, S.; Lee, J. Deep learning toolbox for automated enhancement, segmentation, and graphing of cortical optical coherence tomography microangiograms. Biomed. Opt. Express 2020, 11, 7325–7342. [Google Scholar] [CrossRef] [PubMed]
  23. Hemangioma. Available online: https://my.clevelandclinic.org/health/diseases/23365-hemangioma (accessed on 4 October 2023).
  24. Moiseev, A.; Ksenofontov, S.; Sirotkina, M.; Kiseleva, E.; Gorozhantseva, M.; Shakhova, N.; Matveev, L.; Zaitsev, V.; Matveyev, A.; Zagaynova, E.; et al. Optical coherence tomography-based angiography device with real-time angiography B-scans visualization and hand-held probe for everyday clinical use. J. Biophotonics 2018, 11, e201700292. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Picture of hemangioma that was imaged in this study. The red box is marking the OCT field of view.
Figure 1. Picture of hemangioma that was imaged in this study. The red box is marking the OCT field of view.
Sensors 23 09312 g001
Figure 2. Structure OCT images: en face (a), B-scan cross-sectional (b) and C-scan cross-sectional (c) view. Crosses on each view represent the location of the section for the other two views. All three sections define one point, which is in the vessel that is not visible in the OCT image. The red squares represent an area of 100 μm × 100 µm.
Figure 2. Structure OCT images: en face (a), B-scan cross-sectional (b) and C-scan cross-sectional (c) view. Crosses on each view represent the location of the section for the other two views. All three sections define one point, which is in the vessel that is not visible in the OCT image. The red squares represent an area of 100 μm × 100 µm.
Sensors 23 09312 g002
Figure 3. The OCTA tail artifact is visible as a slowly decreasing signal after the large vessel. The OCTA depth profile exhibits a sharp increase at the position of the vessel and a slow decrease in the signal beneath the vessel. The ideal profile would be a symmetric peak in place of the vessel (sketched in red on the depth profile).
Figure 3. The OCTA tail artifact is visible as a slowly decreasing signal after the large vessel. The OCTA depth profile exhibits a sharp increase at the position of the vessel and a slow decrease in the signal beneath the vessel. The ideal profile would be a symmetric peak in place of the vessel (sketched in red on the depth profile).
Sensors 23 09312 g003
Figure 4. Uncorrected and corrected OCTA images: en face (left), B-scan cross-sectional (middle) and C-scan cross-sectional (right) view. The upper row shows uncorrected OCTA with one example of the tail artifact marked with an arrows. The lower row shows corrected OCTA images with the TAR-TES method and w1 = 0.02. The removal of tail artifacts is balanced. (Note that the color scale is different for each row).
Figure 4. Uncorrected and corrected OCTA images: en face (left), B-scan cross-sectional (middle) and C-scan cross-sectional (right) view. The upper row shows uncorrected OCTA with one example of the tail artifact marked with an arrows. The lower row shows corrected OCTA images with the TAR-TES method and w1 = 0.02. The removal of tail artifacts is balanced. (Note that the color scale is different for each row).
Sensors 23 09312 g004
Figure 5. Corrected OCTA images with the TAR-TES method and w1 = 0.01 (upper row), and w1 = 0.05 (lower row). On the left is the en face view, in the middle is the B-scan cross-sectional view, and on the right is the C-scan cross-sectional view. (Note that the color scale is different for each row).
Figure 5. Corrected OCTA images with the TAR-TES method and w1 = 0.01 (upper row), and w1 = 0.05 (lower row). On the left is the en face view, in the middle is the B-scan cross-sectional view, and on the right is the C-scan cross-sectional view. (Note that the color scale is different for each row).
Sensors 23 09312 g005
Figure 6. Corrected OCTA images with the Mean-Subtraction algorithm for three correction weights. On the left is the en face view, in the middle is the B-scan cross-sectional view, and on the right is the C-scan cross-sectional view. The first row shows the corrected OCTA with the weight that we consider optimal, the second row shows the correction with two-times lower weight, and the third row shows the correction with a 2.5-times higher weight. (Note that the color scale is different for each row).
Figure 6. Corrected OCTA images with the Mean-Subtraction algorithm for three correction weights. On the left is the en face view, in the middle is the B-scan cross-sectional view, and on the right is the C-scan cross-sectional view. The first row shows the corrected OCTA with the weight that we consider optimal, the second row shows the correction with two-times lower weight, and the third row shows the correction with a 2.5-times higher weight. (Note that the color scale is different for each row).
Sensors 23 09312 g006
Figure 7. Corrected OCTA images with the Mean-Subtraction algorithm for three additional correction weights. (Note that the color scale is different for each row).
Figure 7. Corrected OCTA images with the Mean-Subtraction algorithm for three additional correction weights. (Note that the color scale is different for each row).
Sensors 23 09312 g007
Figure 8. Corrected OCTA images with the Step-down Exponential filtering method. On the left is the en face view, in the middle is the B-scan cross-sectional view, and on the right is the C-scan cross-sectional view. The first row shows the corrected OCTA with the weight that we consider optimal, the second row shows the correction with two-times lower weight, and the third row shows the correction with a 2.5-times higher weight. (Note that the color scale is different for each row).
Figure 8. Corrected OCTA images with the Step-down Exponential filtering method. On the left is the en face view, in the middle is the B-scan cross-sectional view, and on the right is the C-scan cross-sectional view. The first row shows the corrected OCTA with the weight that we consider optimal, the second row shows the correction with two-times lower weight, and the third row shows the correction with a 2.5-times higher weight. (Note that the color scale is different for each row).
Sensors 23 09312 g008
Figure 9. Original and corrected en face OCTA images at various depths. A correction was carried out with the TAR-TES method with w1 = 0.02, the Mean-Subtraction algorithm (MSA) with w2 = 3 and the Step-down Exponential filtering method (SDEF) with w3 = 0.02.
Figure 9. Original and corrected en face OCTA images at various depths. A correction was carried out with the TAR-TES method with w1 = 0.02, the Mean-Subtraction algorithm (MSA) with w2 = 3 and the Step-down Exponential filtering method (SDEF) with w3 = 0.02.
Sensors 23 09312 g009
Figure 10. The A-scan plot of original and corrected OCTA—original values (a) and scaled values to the maximum (b). A correction was carried out using the TAR-TES method with w1 = 0.02, the Mean-Subtraction algorithm (MSA) with w2 = 3 and the Step-down Exponential filtering method (SDEF) with w3 = 0.02.
Figure 10. The A-scan plot of original and corrected OCTA—original values (a) and scaled values to the maximum (b). A correction was carried out using the TAR-TES method with w1 = 0.02, the Mean-Subtraction algorithm (MSA) with w2 = 3 and the Step-down Exponential filtering method (SDEF) with w3 = 0.02.
Sensors 23 09312 g010
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Simoncic, U.; Milanic, M. Tail Artifact Removal via Transmittance Effect Subtraction in Optical Coherence Tail Artifact Images. Sensors 2023, 23, 9312. https://0-doi-org.brum.beds.ac.uk/10.3390/s23239312

AMA Style

Simoncic U, Milanic M. Tail Artifact Removal via Transmittance Effect Subtraction in Optical Coherence Tail Artifact Images. Sensors. 2023; 23(23):9312. https://0-doi-org.brum.beds.ac.uk/10.3390/s23239312

Chicago/Turabian Style

Simoncic, Urban, and Matija Milanic. 2023. "Tail Artifact Removal via Transmittance Effect Subtraction in Optical Coherence Tail Artifact Images" Sensors 23, no. 23: 9312. https://0-doi-org.brum.beds.ac.uk/10.3390/s23239312

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop