Research Article  Open Access
Shuo Huang, Le Cheng, Bin Zhu, Ping Zhou, Yu Sun, Bing Zhang, Suiren Wan, "The TimeDomain Integration Method of Digital Subtraction Angiography Images", Computational and Mathematical Methods in Medicine, vol. 2018, Article ID 5284969, 5 pages, 2018. https://doi.org/10.1155/2018/5284969
The TimeDomain Integration Method of Digital Subtraction Angiography Images
Abstract
The clarity improvement and the noise suppression of digital subtraction angiography (DSA) images are very important. However, the common methods are very complicated. An image timedomain integration method is proposed in this study, which is based on the blood flow periodicity. In this method, the images of the first cardiac cycle after the injection of the contrast agent are integrated to obtain the timedomain integration image. This method can be used independently or as a postprocessing method of the denoising method on the signal image. The experimental results on DSA data from an aortic dissection patient show that the image timedomain integration method is efficient in image denoising and enhancement, which also has a good realtime performance. This method can also be used to improve the denoising and image enhancement effect of some common models.
1. Introduction
Circulatory system diseases, such as aortic dissection, have been the focus of medical research [1, 2] for their dangerousness and high incidence. Improving the clarity of the captured medical image helps us to diagnose more accurately, which is an important research field.
Digital subtraction angiography (DSA), as a realtime approach, is commonly employed in the clinical diagnosis of circulatory system disease [3, 4], especially in the realtime surgical monitoring and the medical examination among small branches of blood vessels which are difficult to be measured by other methods. In order to protect the patient, it is important to shorten the shooting time and to reduce the dosage of contrast media when capturing the images.
A lot of image denoising methods have been proposed in recent years, but they are all problematic when applied to the DSA images. For example, the image reconstruction method based on the level set theory [5], wavelet decomposition and reconstruction method [6, 7], Bayesian method [8], and image denoising method based on anisotropic diffusion [9] generally need a long operation time, which cannot meet the realtime requirements of DSA image processing. Moreover, images processed by these approaches are usually not clear enough to show the details such as edges and textures. In 2004, Candes et al. [10] proposed an image denoising method based on the sparse decomposition. On this basis, Needell and Vershynin [11] proposed the regularized orthogonal matching pursuit (ROMP) method; Scholefield and Dragotti [12] used a sparse quadtree decomposition representation to remove the noise in images; Adler et al. employed the shrinkage learning approach to acquire the highresolution reconstruction image [13–17]. However, the operation of these approaches is also very complicated. Therefore, it is necessary to find an image processing method which is more suitable for the realtime analysis of DSA images.
In this work, an image timedomain integration method based on blood flow periodicity has been proposed. In this algorithm, the DSA images of the first cardiac cycle after the injection of the contrast agent are extracted denoised by the wavelet reconstruction method firstly, and then these images are integrated to obtain the timedomain integration image, which is named after the TDI image in this paper. This method contributes to the diagnosis of circulatory system diseases.
2. Materials and Methods
2.1. The Noise Model
The theoretical grayscale at a certain pixel on the jth frame of DSA images can be obtained from the following equation by the Lambert–Beer law [18, 19]:where and are the Xray transmission amount before and after the addition of the contrast agent, respectively. is the volume of blood vessels at pixel . and are the number and amount of substance concentration of contrast agent particles at pixel , respectively. donates the absorption coefficient.
The image quality degrades in the original DSA image, as a result of the limitations on the imaging system’s resolution and the influence of additive noise such as Gaussian noise, which is donated by and can be expressed in as follows [6]:where and represent the point spread function and the additive noise, respectively. Operator “” is the convolution operator. Equation (1b) can be rewritten to matrix form using the block Toeplitz matrix , as shown in the following equation:
It is difficult to solve Equation (1c) when only is given. However, since the grayscale level of a certain pixel is proportional to the number of contrast agent particles in that pixel, the number of contrast agent particles follows the motion pattern of blood. And as for the blood motion pattern, on consideration of the periodicity of human heartbeat, the blood flow rate in human body is also cyclical, which can be expressed in the following equation:where is the velocity field of blood at time . denotes the average velocity field of blood at time , which is the mean flow velocity at the same time in multiple cardiac cycles. is the length of the cardiac cycle. characterizes the changes in flow velocity owing to factors such as the instability of human blood pressure. can be regarded as a zeromeanvalue distribution with a small variance, since patients are under the total anesthesia during the shooting process and their vital signs remain stable. According to the Wilke–Chang equation [20], the free diffusion rate of the contrast agent in the blood is much smaller than the blood flow rate, and thus the contrast agent obeys the same movement law as the blood. Therefore, the periodicity of the blood flow rate can be employed to improve the clarity of DSA images.
2.2. Image Integration
In order to decrease the shoot time and the contrast agent’s injection quantity, images in the first cardiac cycle (donated by the cardiac cycle S) after the injection of the contrast agent are analyzed in this study. Firstly, the time at which the cardiac cycle begins is set to be . Subsequently, the velocity of the ith contrast agent particle in this cardiac cycle is expressed as . Since the motion of the contrast agent particles is consistent with that of the blood, and on consideration of the velocity stability shown in Equation (2), the velocity of the particle at each position on its trajectory can be regarded as a sample of the blood flow field at that location. Therefore, once the substantial number of particles is extracted, the average velocity field of the contrast agent at the pixel in the entire cardiac cycle can be estimated by the mean velocity field value of particles which flowed through that pixel during the calculated cardiac cycle, as shown in the following equation:where represents the time when the ith particle approached pixel .
The total time length that the ith particle appears in pixel during one cardiac cycle satisfies Equation (3b), where represents the distance of the ith particle in the range of pixel and is set to be the magnitude of velocity component which is parallel to the image plane in pixel during that cardiac cycle since the photographing each DSA image can be regarded as a sample of each particle’s location:
A variate is set to represent the absorption capacity of light in the unit time of a single contrast agent particle. After that, the timeweighted grayscale value of the ith particle at pixel , , can be expressed by the following equation:
On combination of Equations (3a)–(3c), can be characterized by the sum of the time integral intensities of particles which have appeared in pixel during the cardiac cycle, as shown in the following equation:where is the average moving distance of the contrast agent particles within that pixel. Since the size of a pixel is small, is approximately equal to the length of each pixel, . According to Equations (1a) and (3d), when the frame rate tends to infinity, the following equation can be obtained:
Equation (3e) demonstrates that the overall timedomain integration value of pixel , , can be expressed as the integral of each picture’s gray value at that position in the entire cardiac cycle. On consideration that the time step is short in the actual case, Equation (3e) can be employed in the calculation of the captured DSA images. Therefore, the relationship shown in Equation (3f) can be established. And the image is named after the timedomain integration image or the TDI image:
Furthermore, to strength the denoising effect, the images are denoised by the median filter before the timedomain integration since Ling’s work [21] shows that the noiseless image is usually insensitive to a median filter. Table 1 shows the specific steps of our method.

3. Results
Figure 1 shows a group of DSA images from a patient with aortic dissection. Figures 1(a1)–1(c1) are the DSA images at , , and , respectively. Figure 1(d1) is the TDI image of the dissecting aneurysm extracted by our algorithm. For the sake of comparison, the image of dissecting aneurysm region in Figures 1(a1)–1(c1) are extracted and then the normalized grayscale histograms of dissecting aneurysm regions among Figures 1(a1)–1(d1) are obtained, which are shown in Figures 1(a2)–1(d2), respectively. Table 2 shows the mean value, standard deviation, and coefficient of variation of DSA images of the dissecting aneurysm region in one cardiac cycle of the patient in Figure 1.

Table 2 shows that the coefficient of variation of the TDI image is higher than all the DSA images in that cardiac cycle, which means that our TDI image can enhance the details. Moreover, the gap between the peaks in Figure 1(d2) is clearer than those in Figures 1(a2)–1(c2), which means that our TDI image has higher resolution than the original images.
4. Discussion
Since the blood flow velocity has a certain degree of uncertainty in the actual situation, which directly influenced the grayscale value of the shot DSA image, Equation (1c) can be rewritten as follows:where donates the grayscale value after the removal of motion randomness.
Equations (3f) and (4) illustrate that the timedomain integration image can be obtained by the following equation:where matrix is defined to characterize the differences between and . A new symbol is defined here, which represents the noisefree timedomain integration image. . Then, the error analysis of Equation (5a) is implemented below.
According to Reference [22], when the cardiac cycle length is 0.8 seconds, the blood flow rate at the aortic inlet, denoted by , follows Equation (5b). The first part of Equation (5b) denotes the ejection phase, which is followed by a brief closure of blood after the closure of the aortic valve. The flow rate in the rest time is 0. Since the time varies when each particle enters the view field, the time when they arrive at the same position on the image also varies. Therefore, the number of images which satisfy is greater than one in most positions. Thus, Equation (5c) can be obtained.
According to Equation (5c), the signaltonoise ratio of image ( is defined as the noisefree image in this calculation) is higher than the mean signaltonoise ratio of original DSA images at each pixel. Therefore, the result of our method contains lower noise in the entire view field, and our method can suppress the noise.
Furthermore, Equation (5d) can be deduced from Equation (3e), where stands for the twonorm of matrix. According to Equation (5d), image is closer to the noisefree TDI image than all of the DSA images, which means that image has the highest clarity. Thus, the image can improve the clarity of the original DSA images.
5. Conclusion
In summary, this study presents a DSA image denoising and enhancement method based on the periodicity of blood flow. Firstly, the DSA images are reconstructed through the median filter, and then DSA images in a cardiac cycle are integrated and the overall timedomain integration image is obtained. According to the mathematical derivation as well as the verification using aortic dissecting aneurysm images, this study demonstrates that the TDI image of the contrast agent has a lower overall noise than original DSA images, and it is also clearer than the original image. This method can contribute to the feature location extraction and disease prophylaxis in circulatory disease, such as the first break’s position extraction of aortic dissecting aneurysms and the analysis of stress distribution on vessel wall [23], which are also our future work.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Authors’ Contributions
Shuo Huang and Le Cheng contributed equally to this work.
Acknowledgments
This work was supported by the National Key R&D Program of China under grant number 2017YFC0112801, the National Natural Science Foundation of China under grant numbers 61127002, 11572087, 61771130, and 31300780, the Research Center for Learning Science of Southeast University under grant number 3207038391, and the School of Biological Sciences and Medical Engineering under grant number 3207037434.
References
 R. Erbel, F. Alfonso, C. Boileau et al., “Diagnosis and management of aortic dissection: recommendations of the task force on aortic dissection, European society of cardiology,” European Heart Journal, vol. 22, no. 18, pp. 1642–1681, 2001. View at: Publisher Site  Google Scholar
 A. White, J. Broder, J. MandoVandrick, J. Wendell, and J. Crowe, “Acute aortic emergencies–part 2: aortic dissections,” Advanced Emergency Nursing Journal, vol. 35, no. 1, pp. 28–52, 2013. View at: Publisher Site  Google Scholar
 J. O. Balzer, M. Doss, A. Thalhammer, H.G. Fieguth, A. Moritz, and T. J. Vogl, “Urgent thoracic aortal dissection and aneurysm: treatment with stentgraft implantation in an angiographic suite,” European Radiology, vol. 13, no. 10, pp. 2249–2258, 2003. View at: Publisher Site  Google Scholar
 E. Sueyoshi, H. Nagayama, I. Sakamoto, and M. Uetani, “Carbon dioxide digital subtraction angiography as an option for detection of endoleaks in endovascular abdominal aortic aneurysm repair procedure,” Journal of Vascular Surgery, vol. 61, no. 2, pp. 298–303, 2015. View at: Publisher Site  Google Scholar
 S. Osher and J. A. Sethian, “Fronts propagating with curvaturedependent speed: algorithms based on HamiltonJacobi formulations,” Journal of Computational Physics, vol. 79, no. 1, pp. 12–49, 1988. View at: Publisher Site  Google Scholar
 S. Wan, B. I. Raju, and M. A. Srinivasan, “Robust deconvolution of highfrequency ultrasound images using higherorder spectral analysis and wavelets,” IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, vol. 50, no. 10, pp. 1286–1295, 2003. View at: Publisher Site  Google Scholar
 S. G. Mallat, “Multifrequency channel decompositions of images and wavelet models,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 37, no. 12, pp. 2091–2110, 1989. View at: Publisher Site  Google Scholar
 M. Lebrun, A. Buades, and J. M. Morel, “A nonlocal bayesian image denoising algorithm,” SIAM Journal on Imaging Sciences, vol. 6, no. 3, pp. 1665–1688, 2013. View at: Publisher Site  Google Scholar
 P. Perona and J. Malik, “Scalespace and edge detection using anisotropic diffusion,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 12, no. 7, pp. 629–639, 1990. View at: Publisher Site  Google Scholar
 E. J. Candes, J. Romberg, and T. Tao, “Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information,” IEEE Transactions on Information Theory, vol. 52, no. 2, pp. 489–509, 2004. View at: Publisher Site  Google Scholar
 D. Needell and R. Vershynin, “Uniform uncertainty principle and signal recovery via regularized orthogonal matching pursuit,” Foundations of Computational Mathematics, vol. 9, no. 3, pp. 317–334, 2009. View at: Publisher Site  Google Scholar
 A. Scholefield and P. L. Dragotti, “Image restoration using a sparse quadtree decomposition representation,” in Proceedings of 2009 IEEE International Conference on Image Processing, pp. 1473–1476, IEEE, Cairo, Egypt, November 2009. View at: Google Scholar
 A. Adler, Y. HelOr, and M. Elad, “A shrinkage learning approach for single image superresolution with overcomplete representations,” in Proceedings of 2010 European Conference on Computer Vision, pp. 622–635, Springer Berlin Heidelberg, Crete, Greece, September 2010. View at: Google Scholar
 M. Elad, M. A. T. Figueiredo, and Y. Ma, “On the role of sparse and redundant representations in image processing,” Proceedings of the IEEE, vol. 98, no. 6, pp. 972–982, 2010. View at: Publisher Site  Google Scholar
 D. Needell and R. Vershynin, “Signal recovery from incomplete and inaccurate measurements via regularized orthogonal matching pursuit,” IEEE Journal of Selected Topics in Signal Processing, vol. 4, no. 2, pp. 310–316, 2010. View at: Publisher Site  Google Scholar
 J. L. Starck, J. M. Fadili, and C. Chesneau, “Image deconvolution by stein block thresholding,” in Proceedings of 2009 IEEE International Conference on Image Processing, pp. 1329–1332, Cairo, Egypt, November 2009. View at: Google Scholar
 D. Tzikas, A. Likas, and N. Galatsanos, “Variational bayesian blind image deconvolution with studentT priors,” in Proceedings of 2007 IEEE International Conference on Image Processing, pp. I109–I112, San Antonio, TX, USA, September 2007. View at: Google Scholar
 M. Elizabeth, Concise Medical Dictionary, Oxford University Press, Oxford, UK, 2015.
 W. D. Jeans, “The development and use of digital subtraction angiography,” British Journal of Radiology, vol. 63, no. 747, pp. 161–168, 1990. View at: Publisher Site  Google Scholar
 C. R. Wilke and P. Chang, “Correlation of diffusion coefficients in dilute solutions,” AIChE Journal, vol. 1, no. 2, pp. 264–270, 1955. View at: Publisher Site  Google Scholar
 J. Ling and A. C. Bovik, “Smoothing lowSNR molecular images via anisotropic mediandiffusion,” IEEE Transactions on Medical Imaging, vol. 21, no. 4, pp. 377–384, 2002. View at: Publisher Site  Google Scholar
 C. Y. Wen, A. L. Yang, and J. W. Chai, “Investigation of pulsatile flow field in healthy thoracic aorta models,” Annals of Biomedical Engineering, vol. 38, no. 2, pp. 391–402, 2010. View at: Publisher Site  Google Scholar
 J. Charonko, S. Karri, J. Schmieg, S. Prabhu, and P. Vlachos, “In vitro comparison of the effect of stent configuration on wall shear stress using timeresolved particle image velocimetry,” Annals of Biomedical Engineering, vol. 38, no. 3, pp. 889–902, 2010. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Shuo Huang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.