Next Article in Journal
Design and Enhancement of a Fog-Enabled Air Quality Monitoring and Prediction System: An Optimized Lightweight Deep Learning Model for a Smart Fog Environmental Gateway
Next Article in Special Issue
Disturbances in Electrodermal Activity Recordings Due to Different Noises in the Environment
Previous Article in Journal
Impact of Perception Errors in Vision-Based Detection and Tracking Pipelines on Pedestrian Trajectory Prediction in Autonomous Driving Systems
Previous Article in Special Issue
Smart Sensing Chairs for Sitting Posture Detection, Classification, and Monitoring: A Comprehensive Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reliable Stenosis Detection Based on Thrill Waveform Analysis Using Non-Contact Arteriovenous Fistula Imaging

by
Rumi Iwai
1,
Takunori Shimazaki
2,3,
Jaakko Hyry
4,
Yoshifumi Kawakubo
2,
Masashi Fukuhara
5,
Hiroki Aono
5,
Shingo Ata
3,
Takeshi Yokoyama
6 and
Daisuke Anzai
4,*
1
Department of Dental Anesthesiology, Kyushu University Hospital, Fukuoka 819-0395, Japan
2
Department of Clinical Engineering, Faculty of Health Care, Jikei University of Health Care Sciences, Osaka 532-0003, Japan
3
Graduate School of Engineering, Osaka Metropolitan University, Osaka 558-8585, Japan
4
Graduate School of Engineering, Nagoya Institute of Technology, Nagoya 466-8555, Japan
5
Department of Clinical Engineering, Shikoku Central Hospital of the Mutual Aid Association of Public School Teachers, Ehime 799-0193, Japan
6
Department of Dental Anesthesiology, Faculty of Dental Science, Kyushu University, Fukuoka 812-8582, Japan
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(15), 5068; https://doi.org/10.3390/s24155068
Submission received: 31 May 2024 / Revised: 1 August 2024 / Accepted: 2 August 2024 / Published: 5 August 2024
(This article belongs to the Special Issue Advanced Non-Invasive Sensors: Methods and Applications)

Abstract

:
Hemodialysis therapy is an extracorporeal circulation treatment that serves as a substitute for renal function. In Japan, patients receive this efficient four-hour treatment, three times per week, allowing them to maintain a social life nearly equivalent to that of healthy individuals. Before the treatment, two punctures are performed to establish extracorporeal circulation, and a high blood flow rate is essential to ensure efficient therapy. Specialized blood vessels created through arteriovenous fistula (AVF) surgery are utilized to achieve high blood flow rates. Although the AVF allows safe and efficient dialysis treatment, AVF stenosis leads to a serious problem in dialysis. To early detect this abnormal blood flow, auscultation and palpation methods are widely used in hospitals. However, these methods can only provide qualitative judgment of the AVF condition, so the results cannot be shared among other doctors and staff. Additionally, since the conventional methods require contact with the skin, some issues require consideration regarding infection and low reproducibility. In our previous study, we proposed an alternative method for auscultation using non-contact optical imaging technology. This study aims to construct a reliable AVF stenosis detection method using Thrill waveform analysis based on the developed non-contact device to solve the problem with the contact palpation method. This paper demonstrates the performance validation of the non-contact imaging in the normal AVF group (206 total data, 75 patients, mean age: 69.1 years) and in the treatable stenosis group (107 total data, 17 patients, mean age: 70.1 years). The experimental results of the Mann–Whitney U test showed a significant difference ( p = 0.0002 ) between the normal and abnormal groups, which indicated the effectiveness of the proposed method as a new possible alternative to palpation.

1. Introduction

Approximately 4 million patients receive kidney replacement therapy (KRT) all over the world, of which 69% are hemodialysis patients [1]. Hemodialysis treatment removes waste products, such as the final breakdown products of dietary proteins and metabolic products from muscle metabolism. It also corrects electrolytes in the blood and removes excess water. The principle in hemodialysis is diffusion and its efficiency is important. A filter called a dialyzer is used, which consists of 10,000 hollow fibers, the sides having countless holes of a nanometer size. Diffusion occurs through the hollow fibers when blood is passed inside these fibers and passing dialysate occurs in the opposite direction on the outside of the fiber. A stable and high blood flow rate is essential for improving the diffusion efficiency.
Because hemodialysis requires extracorporeal circulation at high blood flow rates, an arteriovenous fistula (AVF) is mainly used in Japan, which is an anastomosis between the radial artery (hereafter referred to as the artery) and the cephalic vein (hereafter referred to as the vein) [2,3,4,5,6,7]. This allows blood to flow from arteries with high internal pressure to veins with low internal pressure, causing veins to become distended and allowing high flow rates to be easily obtained from shallow veins. This enables stable dialysis therapy; however, stenosis of the AVF remains a significant issue due to factors such as venous tortuosity, vascular wall thickening, the presence of venous valves, blood pressure drops during treatment and thrombosis [8,9,10,11]. When stenosis becomes challenging to treat, invasive procedures such as percutaneous transluminal angioplasty (PTA) or venous stenting are performed [12,13,14]. Stenosis not only reduces dialysis efficiency, but also imposes substantial physical and mental burdens on patients as it worsens, significantly decreasing their quality of life (QoL).
Because AVFs are easily stenotic vessels, blood flow is checked by auscultation and palpation before and after treatment. The auscultation is a method of listening to the pulsating and turbulent sounds produced at the anastomosis with a stethoscope, and palpation is a method of palpating the blood vessel by feeling the vibrations of the blood vessels (called Thrill) produced by the pulsating and turbulent flow at the anastomosis from the skin surface. The difference between a normal and stenotic AVF is that, when auscultating a normal AVF, there is a laminar blood flow, i.e., a constant or regular flow—you can hear the alternating strong and weak sounds, as well as during palpating, you can also feel the finer vibrations. However, when auscultating a stenotic AVF, this laminar flow changes into a turbulent one, where the high-frequency sounds of fistula bruit, i.e., a vascular murmur, increase, and during palpating the Thrill becomes smaller.
Both methods are simple and reliable screening methods, but they only allow qualitative judgment and data cannot be shared among others. In addition, contact infection must be taken into consideration because of skin contact, and reproducibility is low because the method varies greatly depending on the way the stethoscope or fingers are applied. Once stenosis is detected, the entire arm is coated with gel, and an ultrasound examination is conducted to identify the location of the stenosis [15,16,17]. This procedure requires skilled technicians and is time-consuming. Therefore, there is a pressing need to address the challenges of existing methods and to achieve early detection of stenosis.
In existing research, there are studies that quantify blood flow using wavelet transforms on sounds acquired from small microphones [18], determine stenosis [19], and assess severe stenosis using the Burg AR model [20]. Various analytical methods using multivariate Gaussian distribution (MGD) models have also been proposed [21]. Furthermore, analyses of blood flow conditions using photoplethysmography and neural networks [22], as well as the classification of stenosis levels, have been conducted [23]. As for non-contact methods, there are approaches for detecting stenosis sites using a CCD camera to capture images from above, and the blood flow is estimated by the slight movements of the AVF by using a digital camera that captures the images from the side.
In our previous research, we proposed a method using a light source and camera to address the limitations of existing techniques and quantify the AVF waveform amplitude, offering an alternative to the stethoscope [24]. This approach allows for the acquisition of comprehensive AVF information from the entire image at once, enabling the detection of stenosis without the need for probing. We also investigated how to determine if a whole arm is either normal or has stenosis (two-class classification) by using a machine learning algorithm, such as the Gradient-Boosting Decision Tree [25]. In this study, we further investigated the frequency bands of Thrill using the same method and validated the effectiveness of stenosis detection based on quantified values.

2. Principle Theory

2.1. Thrill Mechanism

Arteries are composed of a three-layered membrane, with a well-developed tunica media containing several smooth muscles and elastic fibers, making the vessel walls thick and highly elastic. Veins, on the other hand, have immature tunica media and, therefore, have thin and less elastic vessel walls [26]. As blood flows from the artery to the vein via the AVF, the arterial blood that spurts out of the narrow anastomosis hits violently against the thin venous wall, stretching it as it pulsates, as shown in Figure 1. The slight oscillations produced during this process are known as ‘Thrill’ [27]. The waveform of the AVF is observed as a composite of arterial and Thrill waveforms. The Thrill is also used to identify the site of stenosis in clinical practice. When the anastomosis narrows, the blood flow decreases, causing the Thrill to decrease and the pulsation to decrease. However, when the downstream region (body side) narrows, the Thrill decreases and the pulsation increases [28].

2.2. Non-Contact Imaging of the AVF

The arm is positioned so that the AVF is at the centre of the ring-shaped light source, as shown in Figure 2. The color of the light source is the same as in [24], which has shown the effectiveness of blue light. The slight displacement of the AVF due to pulsation and Thrill are captured using a CCD camera as changes in luminance, which are acquired at 40 fps. The saved file format is TIFF, a lossless compression method, and the resulting image file contains AVF waveform data for every pixel (0.08 mm × 0.08 mm). Thus, we obtain the original image data f ( i , j ) at a two-dimensional location [ i , j ] T .

2.3. Moving-Average Filter

As non-contact imaging contains a lot of noise from halation and diffuse reflection on the skin, a moving average filter [29] was used to remove this. As shown in Figure 3, this is a general method of smoothing high-frequency components by weighting an N×M matrix called a kernel. In the moving-average filter method, the kernel is defined as an N × M matrix whose component represents a weight coefficient h ( m , n ) . The original image data f ( i , j ) are convolutionally integrated with the kernel while moving one pixel at a time in a raster scan. Thus, we obtain the moving-averaged value g ( i , j ) , which is expressed with the following equation,
g ( i , j ) = n = 1 N m = 1 M f ( i + m , j + n ) h ( m + 1 , n + 1 ) .
Note that the kernel in this study performs N × N processing, and henceforth, N × N will be expressed as N. In addition, we assume N is an odd number. Since the moving-averaged output is measured in the time domain, the frequency spectrum of g ( i , j ) should be calculated with the Fast Fourier Transformation (FFT). Then, the signal and noise power components can be obtained in each specific frequency range. We define the ratio of signal power to noise power as SNR (signal-to-noise power ratio), which can be used for the detection of the AVF stenosis. It should be noted that the signal and noise frequency ranges are discussed in detail in Section 4.

3. Preliminary Experiment

To quantify the Thrill waveforms, the frequency bands were identified using the following procedure: The AVF waveform is a composite wave of the radial artery and Thrill, and the frequency characteristics of the radial artery waveform revealed the major harmonic components that form the wave. The frequency characteristics of the AVF waveform were used to determine the difference between the major frequency components of the radial artery and to identify the frequency band of Thrill.

3.1. Frequency Response of Arterial Waveforms Using Contact Method

A pulse wave sensor (BH1792GLC, manufactured by Rohm) was placed in contact with the radial artery and applied with a sampling frequency of 40 Hz for 10 s. The frequency response [30] processed by the band-pass filter (BPF) is shown in Figure 4. As indicated by the lines, the spectral peaks appear in integer multiples ( n 6 ) from the fundamental to the nth harmonic, and the peaks become smaller as n increases. It was also assumed to be noise as it was almost flat above the sixth harmonic. The cross-correlation coefficients between the original waveform and each of the waveforms obtained by low-pass filter (LPF) processing with the high cut-off frequency between each harmonic are shown in Figure 5. The cross-correlation coefficient tends to increase continuously and the rate of increase is saturated after 4.6 Hz. It is found that up to the third harmonic are greatly involved in the formation of the waveform, and the contribution of the fourth and subsequent harmonics is small.

3.2. Frequency Response of AVF Using Contact Method

Figure 6a shows the frequency response of the AVF waveform from the contact method filtered in the same way as Figure 4, with a higher peak than the third harmonic peak in the 5.0 Hz–7.5 Hz band. This phenomenon was not observed in the frequency analysis of the radial artery in either contact or non-contact, as shown in Figure 4 and Figure 6b.

3.3. Frequency Response of Arterial Waveforms Using Non-Contact Method

Figure 6b shows the frequency characteristics of the arterial waveforms obtained by the non-contact method through the same filter processing as in Figure 4. Because of the reduced photosensitivity of the non-contact method, the arterial waveforms from the contact method were observed up to the third harmonic, and the peaks after that were nearly of the same power spectrum.

3.4. Frequency Response of AVF Using Non-Contact Method

The spectrum of the radial artery decreased as the order n increased, but as shown in Figure 6c, the spectrum of the AVF was confirmed to increase in the band around 5–9 Hz. As this phenomenon was not observed in the radial artery, and since this band did not contribute to formation of the waveform, we derived this as the Thrill band, and for this study, we also defined 5–10 Hz as the Thrill frequency band.

4. Method

4.1. Experiment Procedure

As shown in Table 1, there were a total of 206 dialysis data (75 patients, mean age 69.1 years) with normal AVF and 107 dialysis data (17 patients, mean age 70.1 years) with treatable stenosis. In the experiment process the patients were placed in a resting position and their arm was inserted into the insertion port as shown in Figure 7, after which the images were taken for 10 s.
The setup and processing are shown in Figure 8, where the blue light with an adjusted brightness is captured by the camera as light reflecting from the arm. These received data also contain a lot of halation from common reflection components, which are generated from the midline and protrusions in the arm, as well as having diffuse refraction noise generated from the skin’s surface. Before the camera captures the images, we need to attenuate the halation with a polarizing filter. The data are then captured continuously and converted to TIFF image files. Subsequently, the noise generated by the diffuse reflection is mostly removed by using a moving average method, resulting in a clear AVF waveform. Then, the AVF waveform is converted from the time domain to the frequency domain by using FFT, and the peak values of each spectrum of the Thrill and noise are detected, so that we can then obtain the SNR.
The general criteria for determining AVF stenosis in ultrasound echography are the brachial artery flow volume (FV) ≤ 350 mL/min and the resistance index (RI) ≥ 0.6–0.7 [31,32,33,34], which are the standards of the participating hospitals in this study. So, FV ≤ 350 mL/min and RI ≥ 0.68 subjects who met these criteria were included and classified as the stenosis group in this study. The configuration of this software and hardware was the same [24] as in the previous study. Henceforth, the non-contact imaging method of this study will be referred to as the proposed method.

4.2. Quantification

As shown in Figure 9, 101 different AVF waves were obtained for one subject by moving the average processing with the kernel size N in an odd number ranging from 1 to 201. Next, after extracting the Thrill component with BPF with the low and high cut-off frequencies of 5Hz and 10Hz, respectively, 256 frame size points were processed with a Hanning window and the median of each peak in the power spectrum obtained by FFT calculated to P e a k t h r i l l ( N ) was obtained. Here, N is defined in Equation (1) as the kernel size of the moving-average filter. The P e a k n o i s e ( N ) is processed in the same way with a BPF (low cut-off frequency 15 Hz, high cut-off frequency 20 Hz) to obtain P e a k n o i s e ( N ) , using the following equation:
S N R t h r ( N ) = P e a k t h r i l l ( N ) P e a k n o i s e ( N )

4.3. Determination of S N R t h r ( N o p t ) and Statistical Determination of Normal and Stenotic AVFs

In S N R t h r ( N ) between the two groups, the normal group and the stenosis group, the optimal kernel N o p t is the N that is the largest difference between the two groups, and the optimal SNR is S N R t h r ( N o p t ) . After confirming the normality of S N R t h r ( N o p t ) between the normal group and the stenosis group using the Shapiro–Wilk test, the difference between the two groups was statistically evaluated.

5. Results

5.1. Effects of the Moving Average Process

As shown in Figure 10, at a kernel size of N = 1 , regular waveforms cannot be visually identified, and the AVF waveform is also indiscernible. High-frequency noise was also observed (indicated by the red triangle markers). At a kernel size of N = 101 , a slight periodicity became observable. The waveforms at the high-frequency noise did not diminish even with application of the moving average. At a kernel size of N = 201 , periodicity became clearly evident, and the amplitude increased significantly. The high-frequency noise also did not disappear. Thus, as the kernel size N increased, a reduction in high-frequency noise was observed, leading to an increase in the regularity of the AVF waveform. However, the high-frequency noise did not wholly disappear.

5.2. Determining N o p t

Figure 11 shows a graph representing the normalized S N R t h r for the normal AVF group ( S N R t h r n o r ), the normalized S N R t h r for the stenotic AVF group ( S N R t h r s t e ), and their residual values, plotted for each kernel N. S N R t h r n o r is represented by the blue line, while S N R t h r s t e is represented by the red line. The residuals of these SNRs are shown as dashed lines, and the black line represents the smoothed approximation curve of these residuals (6th order polynomial, and coefficient of determination R 2 = 0.93 ). Both S N R t h r n o r and S N R t h r s t e exhibited an increasing trend similar to a logarithmic function. The smoothed residual curve peaks at N o p t = 21 , which has been determined as the optimal kernel size for quantitatively distinguishing between normal and stenosis conditions.

5.3. Quantitative Assessment of Normal and Stenotic Groups

Figure 12a,b shows a comparison of S N R t h r n o r ( N o p t ) in the normal group and S N R t h r s t e ( N o p t ) in the stenosis group. Figure 12a shows the 75th percentile values for the normal and abnormal groups were 0.24 and 0.17, respectively, indicating a higher value for the normal group. The median values were 0.14 and 0.09 for the normal and abnormal groups, respectively, showing a slightly higher value for the normal group. The 25th percentile values were 0.07 and 0.05 for each group, respectively, indicating nearly identical values between the two groups. The Shapiro–Wilk normality test showed that p = 0.0001 , which did not follow a normal distribution, so a Mann–Whitney U test was performed, and a significant difference of p = 0.0002 was observed between the two groups.
The interquartile ranges for the normal and stenosis groups were 0.17 and 0.12, respectively. Figure 12b consists of a probability density, box and dot plots. The probability density shows the data distribution as a smooth shape, while the dot plot shows the plot from Figure 12a with an outlier added into it. The dot plot also shows the frequency distribution. When compared against the stenosis group, the distribution of the normal group is wider from the median and to the higher values, while the stenosis group is concentrated around the median.

6. Discussion

6.1. Regularity of the Thrill Waveforms

Improvement in waveforms was observed with the increase in kernels, confirming the effectiveness of the moving average method. Although Thrill is said to be random waveforms in clinical practice, the AVF wave showed an increase in the peak value shown in Figure 6c, which was not seen in the frequency response shown in Figure 4 and Figure 6b. This suggests that the Thrill waveforms do not only consist of irregular waves, but also contain regular waveforms.
The arterial wall is subjected to stretching normal to the blood flow and to wall shear stress (WSS) in the tangential direction. Stretching is caused by pulsation and WSS is caused by the spiral laminar flow (SLF) [35,36]. SLF is a widely observed phenomenon in human and animal blood vessels [37] and moderate WSS has been reported to contribute to the suppressive effect of atheroma [38,39,40] and renal damage due to renal artery stenosis [41].
In recent years, the presence of the SLF has also been confirmed in AVFs, and its involvement in vascular maturation has been reported [42,43]. In Marie et al. [44], SLF was observed in 80% of the group capable of sensing Thrill among the 56 subjects with AVF and it was found that Thrill was characterized by spirals rather than turbulence. Their discussion suggests that the frequency band of Thrill waveforms defined in this study may provide evidence of the periodicity created by SLF.
The experimental results regarding the amplitude components of AVF waves from previous studies [24] serve as a quantification method for the stretching of the vessels, which, if demonstrated for the quantification of the WSS component of the Thrill waveforms in this study, would enable the quantification of the two actions of the specific vessel, AVF. However, care should be taken as conditions such as the AVF and blood flow rate differ in Japan compared to other countries, especially the US [2]. It is essential to increase the number of subjects and investigate the regularity and the frequency band of Thrill waveforms defined in this study in the future.

6.2. Validity of the Proposed Method

A comparison of S N R t h r n o r ( N o p t ) in the normal group and S N R t h r s t e ( N o p t ) in the stenosis group, where the frequency band of the Thrill waveforms was 5–10 Hz, showed statistical differences, which were not evident in the box plots. The reason for this may be that the stenosis group in this experiment were subjects who could be treated normally, so the presence or absence of a small amount of stenosis could be shown in the statistics. This paper aimed to identify the frequency components of the Thrill, radial artery and noise. In our investigation, the Thrill component at 5–10 Hz was found to depend on the stenotic changes. Given that the Thrill waveform is sensitive to the blood flow, blood flow variation should have an effect on the stenotic changes, which results in a significant difference in the Thrill component. As shown in Figure 13, to further support the effectiveness of the proposed method, we examined the ultrasound tomographic images of the AVFs of two subjects from the stenosis group and found that subject A had severe stenosis in the brachial vein and subject B had severe stenosis in the radial cutaneous vein at the elbow, indicated by the arrows. When considering the influence of individual differences of the AVF condition on the quantification of stenosis, the experiment was conducted on 92 patients, so the current results include a variety of AVF conditions. Even in these results, the quantification of stenosis changes significantly, so we believe that the effectiveness of quantification of stenosis can be confirmed even with these differences. As future work, we would like to undertake a quantitative study on the quantification of stenosis due to these differences in each individual AVF condition.

6.3. Limitation of the Proposed Method

Although we have confirmed the effectiveness of the proposed detection, there are some limitations in the method. For example, stenosis ranges from mild to severe, while also requiring treatment such as PTA. In addition, even with the same degree of stenosis, the measurement results might differ based on the location, for example, if the stenosis is near an anastomosis or on the trunk side. This implies that they might be the reason why it was difficult to distinguish them visually on the graph. We believe that as the collection of stenosis data continues, that in the future the results will become clearer from the evaluation degree and the location of the stenosis. In addition, as this experiment was conducted on patients at a Japanese hospital, the experiment captured changes in luminance, so the results may differ even for the same person due to the degree of shading, sweat, skin color and condition. We plan to investigate how the results are affected in various environments in the future.

7. Conclusions

This paper developed reliable stenosis detection based on non-contact AVF imaging, which made use of optical technology. To evaluate the efficiency of the developed detection system, Thrill components were quantified with 92 patients and resulted in 313 data elements, using non-contact imaging in the normal group (206 total data, 75 patients) and in the treatable stenosis group (107 total data, 17 patients). A significant difference was observed in the Mann–Whitney U test between the normal group and the stenosis group, which demonstrates the effectiveness of the developed system.
The ultimate goal of this research is to perform non-contact screening with higher accuracy than auscultation and palpation methods, and to image the stenotic area. Additionally, the non-contact imaging method also allows us to identify the stenosis area and the degree of stenosis by looking at the color changes in the tone, similar to a heat map. This makes it possible to take action even at minor stenosis points. However, as a preliminary study, we cannot offer any additional physiological and pathological reasons for the AVF fluctuations in this study. As we believe that they are an important aspect to study, we will pursue these factors in future work.
Previous research confirmed the effectiveness of a non-contact quantitative method as an alternative to auscultation, and this study showed that it has the potential to be a new alternative method that can eliminate the disadvantages of the auscultation and palpation methods mentioned above. In the future, it would be necessary to further increase the number of subjects and carefully investigate the regularity and the frequency band of the Thrill waveforms defined in this study, taking into account the aforementioned factors.

Author Contributions

Conceptualization, R.I., Y.K., M.F., H.A., S.A., J.H. and T.S.; methodology, R.I., S.A., T.S. and M.F.; software, R.I., M.F. and H.A.; validation, R.I., T.S., T.Y. and D.A.; formal analysis, R.I.; investigation, R.I., H.A. and S.A.; resources, T.Y.; data curation, R.I., D.A. and T.S.; writing—original draft preparation, R.I., J.H. and T.S.; writing—review and editing, J.H., R.I. and T.S.; visualization, R.I. and J.H.; supervision, T.Y., D.A. and T.S.; project administration, T.Y. and D.A.; funding acquisition, D.A. and T.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by JSPS KAKENHI Grant Number 24K00885, 21H01325, and by JST Moonshot R&D Grant Number JPMJMS2214-06.

Institutional Review Board Statement

The Ethics Committee of Jikei University of Health Care Sciences approved this experiment (Ethics-ce-juhs No. 2022-1).

Informed Consent Statement

The hospital staff provided sufficient explanations based on the Declaration of Helsinki and informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are only available on request from the corresponding author but with certain limitations to maintain ethics, data privacy and anonymity for the test participants.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
AVFArteriovenous fistula
LPFLow-pass filter
KRTKidney replacement therapy
PTAPercutaneous transluminal angioplasty
QoLQuality of life
MGDMultivariate Gaussian distribution
BPFBand-pass filter
RIResistance index
FLFlow volume
WSSWall shear stress
SLFSpiral laminar flow
FFTFast Fourier transform
SNRSignal-to-noise power ratio

References

  1. Bello, A.K.; Okpechi, I.G.; Osman, M.A.; Cho, Y.; Htay, H.; Jha, V.; Wainstein, M.; Johnson, D.W. Epidemiology of haemodialysis outcomes. Nat. Rev. Nephrol. 2022, 18, 378–395. [Google Scholar] [CrossRef] [PubMed]
  2. Pisoni, R.L.; Zepel, L.; Fluck, R.; Lok, C.E.; Kawanishi, H.; Süleymanlar, G.; Wasse, H.; Tentori, F.; Zee, J.; Li, Y.; et al. International differences in the location and use of arteriovenous accesses created for hemodialysis: Results from the Dialysis Outcomes and Practice Patterns Study (DOPPS). Am. J. Kidney Dis. 2018, 71, 469–478. [Google Scholar] [CrossRef] [PubMed]
  3. Santoro, D.; Benedetto, F.; Mondello, P.; Pipitò, N.; Barillà, D.; Spinelli, F.; Ricciardi, C.A.; Cernaro, V.; Buemi, M. Vascular access for hemodialysis: Current perspectives. Int. J. Nephrol. Renov. Dis. 2014, 7, 281. [Google Scholar] [CrossRef] [PubMed]
  4. Dinwiddie, L.C. Overview of the role of a vascular access nurse coordinator in the optimization of access care for patients requiring hemodialysis. Hong Kong J. Nephrol. 2007, 9, 99–103. [Google Scholar] [CrossRef]
  5. Allon, M.; Robbin, M.L. Increasing arteriovenous fistulas in hemodialysis patients: Problems and solutions. Kidney Int. 2002, 62, 1109–1124. [Google Scholar] [CrossRef] [PubMed]
  6. Malovrh, M. Vascular access for hemodialysis: Arteriovenous fistula. Ther. Apher. Dial. 2005, 9, 214–217. [Google Scholar] [CrossRef] [PubMed]
  7. Konner, K.; Nonnast-Daniel, B.; Ritz, E. The arteriovenous fistula. J. Am. Soc. Nephrol. 2003, 14, 1669–1680. [Google Scholar] [CrossRef]
  8. Moist, L.M.; Lee, T.C.; Lok, C.E.; Al-Jaishi, A.; Xi, W.; Campbell, V.; Graham, J.; Wilson, B.; Vachharajani, T.J. Education in vascular access. Semin. Dial. 2013, 26, 148–153. [Google Scholar] [CrossRef] [PubMed]
  9. Cheung, A.K.; Imrey, P.B.; Alpers, C.E.; Robbin, M.L.; Radeva, M.; Larive, B.; Shiu, Y.T.; Allon, M.; Dember, L.M.; Greene, T.; et al. Intimal hyperplasia, stenosis, and arteriovenous fistula maturation failure in the hemodialysis fistula maturation study. J. Am. Soc. Nephrol. 2017, 28, 3005–3013. [Google Scholar] [CrossRef]
  10. MacRae, J.M.; Dipchand, C.; Oliver, M.; Moist, L.; Lok, C.; Clark, E.; Hiremath, S.; Kappel, J.; Kiaii, M.; Luscombe, R.; et al. Arteriovenous access failure, stenosis, and thrombosis. Can. J. Kidney Health Dis. 2016, 3, 2054358116669126. [Google Scholar] [CrossRef]
  11. Gao, M.; Wang, J. Risk factors of arteriovenous fistula stenosis of patients with maintenance hemodialysis. Evid.-Based Complement. Altern. Med. 2022, 2022, 2968122. [Google Scholar] [CrossRef] [PubMed]
  12. Tessitore, N.; Mansueto, G.; Bedogna, V.; Lipari, G.; Poli, A.; Gammaro, L.; Baggio, E.; Morana, G.; Loschiavo, C.; Laudon, A.; et al. A prospective controlled trial on effect of percutaneous transluminal angioplasty on functioning arteriovenous fistulae survival. J. Am. Soc. Nephrol. 2003, 14, 1623–1627. [Google Scholar] [CrossRef] [PubMed]
  13. Christidou, F.P.; Kalpakidis, V.I.; Iatrou, K.D.; Zervidis, I.A.; Bamichas, G.I.; Gionanlis, L.C.; Natse, T.A.; Sombolos, K.J. Percutaneous transluminal angioplasty (PTA) and venous stenting in hemodialysis patients with vascular access-related venous stenosis or occlusion. Radiography 2006, 12, 127–133. [Google Scholar] [CrossRef]
  14. Haage, P.; Vorwerk, D.; Piroth, W.; Schuermann, K.; Guenther, R.W. Treatment of hemodialysis-related central venous stenosis or occlusion: Results of primary Wallstent placement and follow-up in 50 patients. Radiology 1999, 212, 175–180. [Google Scholar] [CrossRef] [PubMed]
  15. Meola, M.; Marciello, A.; Di Salle, G.; Petrucci, I. Ultrasound evaluation of access complications: Thrombosis, aneurysms, pseudoaneurysms and infections. J. Vasc. Access 2021, 22, 71–83. [Google Scholar] [CrossRef] [PubMed]
  16. Lima, A.; Carrilho, P.; Germano, A. Clinical and ultrasound evaluation for hemodialysis access creation. Nefrologia 2022, 42, 1–7. [Google Scholar] [CrossRef]
  17. Zamboli, P.; Fiorini, F.; D’Amelio, A.; Fatuzzo, P.; Granata, A. Color Doppler ultrasound and arteriovenous fistulas for hemodialysis. J. Ultrasound 2014, 17, 253–263. [Google Scholar] [CrossRef]
  18. Sato, T.; Tsuji, K.; Kawashima, N.; Agishi, T.; Toma, H. Evaluation of blood access dysfunction based on a wavelet transform analysis of shunt murmurs. J. Artif. Organs 2006, 9, 97–104. [Google Scholar] [CrossRef]
  19. Wang, Y.N.; Chan, C.Y.; Chou, S.J. The detection of arteriovenous fistula stenosis for hemodialysis based on wavelet transform. Int. J. Adv. Comput. Sci. 2011, 1, 16–22. [Google Scholar]
  20. Chen, W.L.; Lin, C.H.; Chen, T.; Chen, P.J.; Kan, C.D. Stenosis detection using burg method with autoregressive model for hemodialysis patients. J. Med. Biomed. Eng. 2013, 33, 356–362. [Google Scholar]
  21. Sung, P.H.; Kan, C.D.; Chen, W.L.; Jang, L.S.; Wang, J.F. Hemodialysis vascular access stenosis detection using auditory spectro-temporal features of phonoangiography. Med. Biol. Eng. Comput. 2015, 53, 393–403. [Google Scholar] [CrossRef] [PubMed]
  22. Grochowina, M.; Leniowska, L.; Dulkiewicz, P. Application of Artificial Neural Networks for the Diagnosis of the Condition of the Arterio-venous Fistula on the Basis of Acoustic Signals. In Proceedings of the Brain Informatics and Health: International Conference, BIH 2014, Warsaw, Poland, 11–14 August 2014; Springer International Publishing: Cham, Switzerland, 2014; pp. 400–411. [Google Scholar]
  23. Du, Y.C.; Stephanus, A. Levenberg-Marquardt neural network algorithm for degree of arteriovenous fistula stenosis classification using a dual optical photoplethysmography sensor. Sensors 2018, 18, 2322. [Google Scholar] [CrossRef]
  24. Iwai, R.; Shimazaki, T.; Kawakubo, Y.; Fukami, K.; Ata, S.; Yokoyama, T.; Hitosugi, T.; Otsuka, A.; Hayashi, H.; Tsurumoto, M.; et al. Quantification and Visualization of Reliable Hemodynamics Evaluation Based on Non-Contact Arteriovenous Fistula Measurement. Sensors 2022, 22, 2745. [Google Scholar] [CrossRef] [PubMed]
  25. Shimazaki, T.; Kawakubo, Y.; Iwai, R.; Fukuhara, M.; Aono, H.; Mitsudo, J.; Hayashi, Y.; Ata, S.; Yokoyama, T.; Anzai, D. A Study on Stenosis Detection Based on Non-contact Thrill Wave Imaging and Gradient-Boosting Decision Tree. In Proceedings of the 2023 IEEE 17th International Symposium on Medical Information and Communication Technology (ISMICT), Lincoln, NE, USA, 10–12 May 2023; pp. 1–4. [Google Scholar] [CrossRef]
  26. Martini, F. Anatomy and Physiology, 2007 ed.; Rex Bookstore, Inc.: Quezon City, Philippines, 2006; pp. 533–534. [Google Scholar]
  27. Gilpin, V.; Nichols, W.K. Vascular access for hemodialysis: Thrills and thrombosis. J. Vasc. Nurs. 2010, 28, 78–83. [Google Scholar] [CrossRef] [PubMed]
  28. Asif, A.; Roy-Chaudhury, P.; Beathard, G.A. Early Arteriovenous Fistula Failure: A Logical Proposal for When and How to Intervene. Clin. J. Am. Soc. Nephrol. 2006, 1, 332–339. [Google Scholar] [CrossRef] [PubMed]
  29. Rafael, C.G.; Richard, E.W. Digital Image Processing; Prentice Hall: Hoboken, NJ, USA, 2002; pp. 116–124. [Google Scholar]
  30. Kyriacou, P.A.; Allen, J. Photoplethysmography: Technology, Signal Analysis, and Applications; Academic Press: Cambridge, MA, USA, 2021; p. 105. [Google Scholar]
  31. Ogawa, T.; Matsumura, O.; Matsuda, A.; Hasegawa, H.; Mitarai, T. Brachial artery blood flow measurement: A simple and noninvasive method to evaluate the need for arteriovenous fistula repair. Dial. Transplant. 2011, 40, 206–210. [Google Scholar] [CrossRef]
  32. Matsuura, K.; Gotoh, Y.; Sadaoka, S.; Takase, K.; Narimatsu, Y. Guidelines for basic techniques in vascular access intervention therapy (VAIVT). Interv. Radiol. 2018, 3, 28–43. [Google Scholar] [CrossRef]
  33. Malovrh, M. Native arteriovenous fistula: Preoperative evaluation. Am. J. Kidney Dis. 2002, 39, 1218–1225. [Google Scholar] [CrossRef] [PubMed]
  34. Kim, Y.O.; Yang, C.W.; Yoon, S.A.; Chun, K.A.; Kim, N.I.; Park, J.S.; Kim, B.S.; Kim, Y.S.; Chang, Y.S.; Bang, B.K. Access blood flow as a predictor of early failures of native arteriovenous fistulas in hemodialysis patients. Am. J. Nephrol. 2001, 21, 221–225. [Google Scholar] [CrossRef]
  35. Stonebridge, P.; Hoskins, P.; Allan, P.; Belch, J. Spiral laminar flow in vivo. Clin. Sci. 1996, 91, 17–21. [Google Scholar] [CrossRef]
  36. Lu, D.; Kassab, G.S. Role of shear stress and stretch in vascular mechanobiology. J. R. Soc. Interface 2011, 8, 1379–1385. [Google Scholar] [CrossRef] [PubMed]
  37. Ku, D.N. Blood flow in arteries. Annu. Rev. Fluid Mech. 1997, 29, 399–434. [Google Scholar] [CrossRef]
  38. Malek, A.M.; Alper, S.L.; Izumo, S. Hemodynamic shear stress and its role in atherosclerosis. JAMA 1999, 282, 2035–2042. [Google Scholar] [CrossRef] [PubMed]
  39. De Nisco, G.; Hoogendoorn, A.; Chiastra, C.; Gallo, D.; Kok, A.M.; Morbiducci, U.; Wentzel, J.J. The impact of helical flow on coronary atherosclerotic plaque development. Atherosclerosis 2020, 300, 39–46. [Google Scholar] [CrossRef] [PubMed]
  40. Ding, Z.; Fan, Y.; Deng, X.; Zhan, F.; Kang, H. Effect of swirling flow on the uptakes of native and oxidized LDLs in a straight segment of the rabbit thoracic aorta. Exp. Biol. Med. 2010, 235, 506–513. [Google Scholar] [CrossRef] [PubMed]
  41. Houston, J.G.; Gandy, S.J.; Milne, W.; Dick, J.B.; Belch, J.J.; Stonebridge, P.A. Spiral laminar flow in the abdominal aorta: A predictor of renal impairment deterioration in patients with renal artery stenosis? Nephrol. Dial. Transplant. 2004, 19, 1786–1791. [Google Scholar] [CrossRef] [PubMed]
  42. Srivastava, A.; Mittal, V.; Lal, H.; Javali, T.; Patidar, N.; Sureka, S.; Aggarwal, S. Spiral laminar flow, the earliest predictor for maturation of arteriovenous fistula for hemodialysis access. Indian J. Urol. IJU J. Urol. Soc. India 2015, 31, 240. [Google Scholar] [CrossRef] [PubMed]
  43. Cunnane, C.V.; Houston, J.G.; Moran, D.T.; Broderick, S.P.; Ross, R.A.; Walsh, M.T. Spiral Laminar Flow is Associated with a Reduction in Disturbed Shear in Patient-Specific Models of an Arteriovenous Fistula. Cardiovasc. Eng. Technol. 2022, 14, 152–165. [Google Scholar] [CrossRef]
  44. Marie, Y.; Guy, A.; Tullett, K.; Krishnan, H.; Jones, R.G.; Inston, N.G. Patterns of blood flow as a predictor of maturation of arteriovenous fistula for haemodialysis. J. Vasc. Access 2014, 15, 169–174. [Google Scholar] [CrossRef]
Figure 1. An illustration of the AVF waveforms formed from the pulse and Thrill.
Figure 1. An illustration of the AVF waveforms formed from the pulse and Thrill.
Sensors 24 05068 g001
Figure 2. Principle of non-contact imaging.
Figure 2. Principle of non-contact imaging.
Sensors 24 05068 g002
Figure 3. AVF imaging area and the kernel size ( N × M ) of the moving-average filter.
Figure 3. AVF imaging area and the kernel size ( N × M ) of the moving-average filter.
Sensors 24 05068 g003
Figure 4. Frequency response of the contact radial artery waveform.
Figure 4. Frequency response of the contact radial artery waveform.
Sensors 24 05068 g004
Figure 5. Determination of radial artery waveform components in the contact method by cross-validation.
Figure 5. Determination of radial artery waveform components in the contact method by cross-validation.
Sensors 24 05068 g005
Figure 6. (a) Frequency response of the contact AVF waveform, (b) frequency response of the non-contact radial artery waveform, and (c) frequency response of the non-contact AVF waveform. The blue arrows, the red frames, and the dashed lines indicate the peaks of Thrill components, the estimated frequency range of Thrill components, and the power spectrum level of the 3rd harmonic, respectively.
Figure 6. (a) Frequency response of the contact AVF waveform, (b) frequency response of the non-contact radial artery waveform, and (c) frequency response of the non-contact AVF waveform. The blue arrows, the red frames, and the dashed lines indicate the peaks of Thrill components, the estimated frequency range of Thrill components, and the power spectrum level of the 3rd harmonic, respectively.
Sensors 24 05068 g006
Figure 7. Experimental equipment setup.
Figure 7. Experimental equipment setup.
Sensors 24 05068 g007
Figure 8. Setup and analysis flow of the non-contact Thrill quantification system.
Figure 8. Setup and analysis flow of the non-contact Thrill quantification system.
Sensors 24 05068 g008
Figure 9. Procedure for S N R t h r ( N o p t ) .
Figure 9. Procedure for S N R t h r ( N o p t ) .
Sensors 24 05068 g009
Figure 10. Effects of moving average methods on different kernels.
Figure 10. Effects of moving average methods on different kernels.
Sensors 24 05068 g010
Figure 11. Determining N o p t .
Figure 11. Determining N o p t .
Sensors 24 05068 g011
Figure 12. Comparison of normal and stenosis groups by S N R t h r o p t represented by a box plot in (a). Comparison of the distribution of the same group by S N R t h r o p t represented by a rain cloud plot in (b).
Figure 12. Comparison of normal and stenosis groups by S N R t h r o p t represented by a box plot in (a). Comparison of the distribution of the same group by S N R t h r o p t represented by a rain cloud plot in (b).
Sensors 24 05068 g012
Figure 13. Ultrasound echocardiographic stenosis.
Figure 13. Ultrasound echocardiographic stenosis.
Sensors 24 05068 g013
Table 1. Experiment protocol.
Table 1. Experiment protocol.
Number of patients92 (normal AVF: 75, AVF with treatable stenosis: 17)
Number of data313 (normal AVF: 206, AVF with treatable stenosis: 107)
PositionSitting
Exposure time10 s
Light source colorBlue
Frame rate40 fps
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

Iwai, R.; Shimazaki, T.; Hyry, J.; Kawakubo, Y.; Fukuhara, M.; Aono, H.; Ata, S.; Yokoyama, T.; Anzai, D. Reliable Stenosis Detection Based on Thrill Waveform Analysis Using Non-Contact Arteriovenous Fistula Imaging. Sensors 2024, 24, 5068. https://doi.org/10.3390/s24155068

AMA Style

Iwai R, Shimazaki T, Hyry J, Kawakubo Y, Fukuhara M, Aono H, Ata S, Yokoyama T, Anzai D. Reliable Stenosis Detection Based on Thrill Waveform Analysis Using Non-Contact Arteriovenous Fistula Imaging. Sensors. 2024; 24(15):5068. https://doi.org/10.3390/s24155068

Chicago/Turabian Style

Iwai, Rumi, Takunori Shimazaki, Jaakko Hyry, Yoshifumi Kawakubo, Masashi Fukuhara, Hiroki Aono, Shingo Ata, Takeshi Yokoyama, and Daisuke Anzai. 2024. "Reliable Stenosis Detection Based on Thrill Waveform Analysis Using Non-Contact Arteriovenous Fistula Imaging" Sensors 24, no. 15: 5068. https://doi.org/10.3390/s24155068

APA Style

Iwai, R., Shimazaki, T., Hyry, J., Kawakubo, Y., Fukuhara, M., Aono, H., Ata, S., Yokoyama, T., & Anzai, D. (2024). Reliable Stenosis Detection Based on Thrill Waveform Analysis Using Non-Contact Arteriovenous Fistula Imaging. Sensors, 24(15), 5068. https://doi.org/10.3390/s24155068

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