Next Article in Journal
Miniaturized Sensors for Detection of Ethanol in Water Based on Electrical Impedance Spectroscopy and Resonant Perturbation Method—A Comparative Study
Previous Article in Journal
Dilute Polymerization of Aniline on PDMS Substrate via Surface Modification Using (3-Aminopropyl)Triethoxysilane for Stretchable Strain Sensor
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quantification and Visualization of Reliable Hemodynamics Evaluation Based on Non-Contact Arteriovenous Fistula Measurement

1
Faculty of Clinical Engineering, Hyogo Medical University Hospital, Nishinomiya 663-8501, 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 City University, Osaka 558-8585, Japan
4
Department of Medicine, School of Medicine, Kurume University, Fukuoka 830-0011, Japan
5
Department of Dental Anesthesiology, Faculty of Dental Science, Kyushu University, Fukuoka 812-8582, Japan
6
Hemodialysis Center, Sugi Hospital, Fukuoka 837-0916, Japan
7
Hemodialysis Center, Tamaki-Aozora Hospital, Tokushima 779-3125, Japan
8
Department of Clinical Engineering, Faculty of Health and Welfare, Tokushima Bunri University, Kagawa 769-2193, Japan
9
Department of Clinical Laboratory and Laboratory Medicine, Kyushu University Hospital, Fukuoka 812-8582, Japan
10
Foundation for Biomedical Research and Innovation at Kobe, Kobe 650-0047, Japan
11
Hemodialysis Center, Wakaura Central Hospital, Wakayama 641-0054, Japan
12
Graduate School of Engineering, Nagoya Institute of Technology, Nagoya 466-8555, Japan
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(7), 2745; https://doi.org/10.3390/s22072745
Submission received: 2 February 2022 / Revised: 29 March 2022 / Accepted: 31 March 2022 / Published: 2 April 2022
(This article belongs to the Section Biomedical Sensors)

Abstract

:
The condition of arteriovenous fistula (AVF) blood flow is typically checked by using auscultation; however, auscultation should require a qualitative judgment dependent on the skills of doctors, and further attention to contact infection is required. For these reasons, this study developed a non-contact and non-invasive medical device to measure the pulse wave of AVFs by applying optical imaging technology. As a first step toward realization of the quantification judgment based on non-contact AVF measurement, we experimentally validated the developed system, whereby the hemodynamics of 168 subjects were visually and quantitatively evaluated based on clinical tests. Based on the evaluation results, the fundamental statistical characteristics of the non-contact measurement, including the average and median values, and distribution of measured signal-to-noise power ratio, were demonstrated. The clinical test results contributed to the future construction of quantified criteria for the AVF condition with the non-contact measurement.

1. Introduction

Currently, ~3.5 million patients are receiving hemodialysis therapy worldwide [1]. In Japan, ~340,000 dialysis patients were treated in 2019 [2], which implies that 1 in 380 people in Japan is receiving treatment. Hemodialysis therapy is one of the most common methods for blood purification, in which blood is pumped out of the body and passed through a purifier, known as a dialyzer, to remove urinary toxins, regulate electrolyte balance, and exclude excess fluid (Figure 1). The basic cycle of the treatment is 4 h per session, with 3 sessions per week. To perform the hemodialysis therapy, a unique blood vessel, known as the arteriovenous fistula (AVF), is used. This vessel is surgically created by anastomosing the flexural cutaneous vein and artery in distal and radial positions. As the intravascular pressure is higher in the radial artery than in the flexor cutaneous vein and the central resistance is lower than the peripheral vascular resistance at the fingertip, the arterial blood flows into the flexor cutaneous vein. These phenomena cause turgid veins and hyperreflux, which generate in an access line for blood delivery and withdrawal.
As the venous wall is thinner and softer than the arterial wall [3], the turbulent flow at the anastomotic site causes the venous walls to collide with each other. As a result, a fistula bruit sound known as a vascular murmur and a vibration known as a thrill [4] occur. The thrill can be measured by lightly touching the anastomosis with a finger. The AVF pulse wave W AVF ( t ) can be decomposed as
W AVF ( t ) = W h e a r t b e a t ( t ) + W t h r i l l ( t )
where W h e a r t b e a t ( t ) and W t h r i l l ( t ) are the pulse wave components of the heartbeat and thrill, respectively.
As the AVF often becomes stenotic [5], the AVF condition is checked by auscultation before and after thedialysis therapy, especially in Japan [6]. Auscultation is an excellent initial screening method owing to its simplicity and reliability. However, it can only help make qualitative decisions. In other words, the auscultation method confirms two kinds of AVF states: normal or abnormal condition. There are other issues as well. For exmple, auscultation results are affected by the pressure applied to the AVF vessel, it occasionally misses the sound of stenosis, the obtained data cannot be shared with co-workers, even if the co-workers are in the same hospital, and the decision can only be made by the staff who perform the measurement, and careful attention to contact infection is necessary because of the direct contact required with the patient.
To solve these problems, existing studies [7,8,9,10,11,12,13,14] developed a system with an electronic stethoscope to detect a fistula bruit (shunt sound) by using wavelets [7,8,9,10,11,12] and deep learning [13,14]. However, electronic stethoscope-based detection is sensitive to environmental noise [15]. In other related studies [15,16,17], the AVF pulse wave acquired by a photoplethysmogram was analyzed based on wavelet [15] and neural network analysis [16,17]. The photoplethysmogram still depends on the amount of pressure exerted by the sensor on the skin while in contact. Kamiyama et al. attempted to detect the stenosis site by capturing its image from above by using a charge-coupled device (CCD) camera [18]. Zhu et al. estimated the blood flow amount based on the weak AVF movement by capturing the image from a side view by using a digital camera [19].
Optical technology, e.g., photographic imaging based on a CCD camera, is efficient in realizing the non-contact and non-invasive measurement of the AVF condition. The advantage of the non-contact evaluation is to avoid concerns about infections because dialysis patients in particular are often infected with hepatitis viruses. In addition, contact AVF measurement should depend on the contact condition, so the non-contact measurement can be expected to solve the problem of the contact condition. However, in existing studies, the validation was limited because the phantom evaluation was conducted [18], and clinical tests with a sufficient number of subjects were not carried out [19]. Furthermore, the quantification of the AVF hemodynamic condition was rarely discussed, which is required to share the diagnostic results. Therefore, this study aims to develop a non-contact device to quantitatively and visually measure the AVF hemodynamic conditions by applying optical technology with color mapping. To validate the developed system, we conducted clinical experiments to clarify hemodynamics in 168 subjects. Then, a moving-average filter to extract AVF pulse waves with high signal-to-noise power ratios (SNRs) is required to quantify the AVF condition. In addition, this paper demonstrates the effectiveness of hemodynamic visualization and quantification based on the clinical test results, which implies that AVF quantification can solve the problems of the conventional qualitative diagnosis and significantly improve hemodialysis therapy.

2. Principle of Non-Contact AVF Pulse Wave Measurement

Moving-Average Filter

Let us explain the principle of the quantification of the AVF pulse wave. We apply moving-average filtering to the measured image data. 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 ) as
g ( i , j ) = n = 1 N m = 1 M f ( i + m , j + n ) h ( m + 1 , n + 1 ) .
The purpose of moving-average filtering is to suppress scattered reflections and halos in the original image data. As the imaging data include pulse wave data for each pixel, a large amount of pulse wave data can be obtained at the same time. Therefore, an increase in the kernel exhibits an effect similar to that of additive averaging. The additive averaging method is effective in extracting weak biological signals from the superimposed noise. Similar techniques were applied to evoke potential measurements [20,21].

3. Methods

3.1. Experimental Protocol

In the clinical test, 168 subjects were rested in a seating position and imaged twice at 40 fps (frames per second) for 10 s at two points close to the apex of the anastomosis by using the developed non-contact imaging system (Figure 2). Table 1 and Table 2 show information and statics of the patients. The file format obtained from the AVF imaging was a raw image format with RGB information for each pixel. Figure 3 shows the setup and analysis flow of the non-contact AVF measurement system. In accordance with the Declaration of Helsinki, sufficient attention was paid to the protection of all subjects. The approval was obtained from the ethics committee of the hospital. The hardware and software specifications of the non-contact imaging system are listed in Table 3.

3.2. Hemodynamics Visualization Method Based on Color Mapping

The visualization of hemodynamics helps evaluate the AVF condition. In this study, color-mapping technology has been applied for visualization, wherein the blood flow direction, which is ejected from the anastomosis toward the trunk side, can be confirmed. Specifically, hemodynamics are visualized by using time-series data in which changes in luminance of the pixels of the raw image format are converted to a 256-level color tone (Figure 4).

3.3. Determination of the Optimal Kernel Size and Quantification

The AFV pulse wave waveform was obtained by increasing the kernel size of the moving average filter from 1 pixel (0.08 nm) to 201 pixels (16.08 nm) at every interval of 2 pixels from a 1 × 1 matrix to a 201 × 201 matrix. Figure 5 shows the relationship between the AVF imaging area and kernel size of the moving-average filter.
Next, let us explain the quantification of hemodynamics for AVF evaluation. The quantification based on the obtained image data is performed by using the following procedure. Because the typical heart rate is between 60 bpm and 100 bpm [22], the lower and higher cutoff frequencies of the band-pass filter (BPF) to extract the pulse wave component is set to 1 Hz and 5 Hz, respectively, based on our previous study [23]. Moreover, the lower and higher cutoff frequencies of the BPF to extract the noise component are set to 10 Hz and 19 Hz, respectively, so that the thrill wave frequency band can be ignored. A higher cutoff frequency of 19 Hz is set because of the Nyquist frequency. The obtained waveforms have been processed by using a fast Fourier transform with a sample size of 256 to realize each BPF to divide the waveform into pulse wave and noise components. Then, we have calculated the SNR as the ratio of the signal (pulse wave) and noise powers.
To calculate the SNR, the method to determine signal and noise components should be discussed. In this study, as noise may be mixed with the AVF component, the sum of the powers at the fundamental frequency and second harmonic components is used as the pulse wave domain band (Figure 6). In addition, the sum of the first- and second-largest signals at 10–19 Hz is considered as the noise domain. The optimal kernel size is determined to achieve the highest SNR; therefore, the quantification of the AVF condition can be defined as:
S N AVF = max ( AVF s / AVF n )
where AVF s and AVF n denote the powers of the pulse wave and noise components, respectively.

4. Results

4.1. Preliminary Evaluation: Effect of Changes in Luminance, Hue, and Color Saturation on Pulse Waveform

Next, we investigated the optimal measure based on the AVF imaging data. Figure 7 shows a comparison among the observed pulse wave waveforms with respect to changes in luminance, hue, and saturation [24] by irradiating white light on the palms of a healthy subject. The evaluation measure of "luminance” is the luminous flux of a surface light source when viewed from an observation point that is expressed in cd / m 2 . The word "hue” defines the characterization of a color, and "saturation” is the degree of vividness of a color. The measure of "lightness” is an expression similar to luminance; however, the measure of lightness’ is excluded in the evaluation process because it is a psychophysical quantity related to vision that is recognized by reflected light.
Based on the comparison results, as shown in Figure 7, it is confirmed that the pulse pressure is calculated based on the luminance variation (measured by the peak-to-peak value of the waveform), which demonstrated the largest value among the three results, and the superimposed noise was relatively small. Hence, we decided to adopt the luminance variation of the AVF pulse wave measurement in the developed non-contact AVF measurement system.

4.2. Preliminary Evaluation: Optimal Wavelength of Imaging

In commercially available watch-type healthcare terminals, as well as in our previous studies, green light is an effective light source for contact measurement methods [25,26,27,28]. However, for non-contact AVF pulse wave imaging, the wavelength for reliable imaging has not yet been clarified. To determine the optimal wavelength for non-contact AVF imaging, it is necessary to conduct preliminary experiments by using several wavelengths, i.e., green light and blue light.
Figure 8 shows the measurement of the pulse wave waveform of the anastomotic area with blue light. In the experiments, effective waveforms were obtained with green land blue lights. However, as a result of three-dimensional image processing of the artery using a photometric method, it became difficult to observe three-dimensional imaging with green light. This is because the photometric method [29,30] is a technique to synthesize three-dimensional images from planar images captured from different directions. As the wavelength of the photometric method is sufficient accurate for small veins on the surface layer, it is difficult to create three-dimensional images owing to the longer wavelength compared to that of blue light. Because this phenomenon is unique to non-contact AVF imaging, we decided to use blue light as the source.

4.3. Preliminary Evaluation: Validation of Non-Contact AVF Imaging

Finally, to validate the pulse wave measured by the non-contact AVF imaging device, we confirmed how a large correlation between the radial artery waveforms and the developed non-contact imaging of a bidirectional Doppler flowmeter (ES-100V3, Hadeco, medical device certification number 21500BZZ00131000). Note that this Doppler flowmeter is based on contact measurement. We carried out the measurement with one healthy subject as shown in Figure 9. As shown in Figure 10, the developed non-contact system ensures the peak period corresponds to the commercial medical device with contact measurement. In other words, the blood flow meter matched the systolic cycle of the heart with regard to the measurement data. Based on the preliminary evaluations, we confirmed that the developed non-contact system has a fundamental function to measure the systolic cycle of the heart.

4.4. Color Mapping-Based Hemodynamic Visualization

Figure 11 shows the continuous AVF imaging results for one cycle (from 2.600 s to 3.360 s after the start of imaging) of a subject. At 2.850 s, the AVF was indicated in blue, and the anastomosis was tense owing to ejection. At 3.000 s, the blood moved from the anastomosis to the puncture site, and at 3.125 s, the puncture site was represented at its peak. At 3.150 s, the entire forearm turned blue. Finally, at 3.250 s, the intravascular pressure drastically dropped, and the anastomosis and puncture site changed in a reddish cycle, which was observed in other subjects as well.

4.5. Measurement Data with Moving-Average Filtering

Figure 12 shows the AVF pulse waves for the kernel size of the moving-average filter N = 1, 13, and 25. Each pulse wave waveform was normalized for comparison. It was confirmed that the AVF pulse wave was clearly obtained, suppressing the noise component, as the kernel size increased.

4.6. AVF Quantification with Optimal Kernel Size

As shown on the left side of Figure 13, the SNR curve is close to an exponential function as N increases. The highest S N AVF = 0.18 was obtained when the kernel size N = 113 was the optimal value. The normalized distribution of S N AVF with N = 113 is shown on the right side of Figure 13. The median S N AVF value of the normalized distribution was −0.37. Additionally, the validation results demonstrated that the first quartile (25 %tile), second (50 %tile), third (75 %tile), and the fourth quartiles (the maximum value) were −0.77, −0.60, 0.23, and 1.44, respectively.
Note that the quartile can be defined as a type of quartile that divides the number of measurement data into four parts with an equal size. For example, the first quartile (25 %tile) is calculated as the middle number between the minimum and the median of the measurement dataset. Additionally, to calculate the normalized distribution of S N AVF , we balanced the measurement data by the average value. Therefore, the average value (not median value) of normalized S N AVF should be of 0 and span from negative to positive value. We can see from Figure 13 that the median value (2nd quartile) is not different from the average value. Also, the variation of normalized S N AVF ranges from −0.77 to 1.44. It is difficult to construct certain criteria to judge the AVF condition based on only the results; however, the results can be used to realize the quantified evaluation of AVF measurement. For example, this study obtained the S N AVF distribution from the experiment with 168 normal dialysis patients in this study, so these data should be informative to analyze the experiment with non-normal patients in the future.
Furthermore, Figure 14 shows typical examples of the AVF pulse wave measurement results with optimal kernel N = 113 . As can be seen from these results, the AVF pulse wave waveforms were efficiently measured in the developed system. In addition, for providing the statics of the measured data, Figure 15 demonstrates the average and variance of the measured luminance for all patients.

4.7. Observation of the Reversed-Phase Pulse Wave

Although no reversed-phase signal component is observed in a normal pulse waveform, a reversed-phase waveform, as shown in Figure 16, was observed at the AVF periphery for all the subjects.

5. Discussion

5.1. Visualization by Color Mapping

Here, we discuss the reason for the reversed-phase waveform observed in the AVF pulse wave (Figure 16). To explain this, Figure 17 shows the color-mapping results of the AVF pulse waves of the parietal (01) and lateral (02) portions. When the arterial blood is ejected and the intra-anastomotic pressure is maximum, the CCD camera and AVF become closest to each other. Hence, the light received at the parietal part exhibits the largest amount. At point (01) on the parietal, the light-receiving volume reaches its maximum value at 0.4 s; however, the light-receiving volume decreases as the intravascular pressure decreases and reaches its minimum value at 0.9 s. If the principle of waveform acquisition is based on the same absorbance change as in general photoplethysmographs, the minimum value can be obtained at 0.4 s at the top of the head and the maximum value at 0.9 s. For this reason, it can be considered that the non-contact AVF pulse wave is a displacement change caused by the vibration of the CCD camera and AVF.
However, the side (02) is a reversed-phase waveform with a minimum value at 0.4 s because the shadow becomes larger as the AVF increases. The maximum value at 0.9 s is observed because the light-receiving amount increases again as the AVF decreases. This phenomenon can be explained by the alternation of the positive reflection of the parietal surface and the shadow of the lateral surface of the light caused by the vibration of the blood vessel. This is why the present method is based on the displacement change rather than the absorbance change. Because this phenomenon does not occur in the contact method, which is unique to non-contact, it should be considered while measuring large vessels that produce shadows in a non-contact manner.
The reason for the increase in luminance over the entire forearm at 3.150 s, as shown in Figure 11, can be explained as follows: as shown in Figure 18, the blood driven from the heart flows through the brachial artery and branches into the radial and ulnar arteries to the peripheral vessel. The blood returns via the cephalic vein to the heart. Thus, in addition to the normal pathway through the peripheral vessel, the bloodstream of the cephalic vein has two pathways that provide a short route toward the peripheral vessel.
As this phenomenon occurred during the anastomotic relaxation phase, normal blood flow was considered through the peripheral vessels. This allows us to identify the site of stenosis not only in the anastomotic area but also in the entire forearm. The visualization of the blood flow from the AVF back to the trunk side based on color mapping can be used to identify the site of stenosis. In this study, the quantization of the color mapping was 8 bits; however, we assume that higher resolution color mapping will be possible by improving quantization.

5.2. AVF Pulse Wave and Stereoscopic Image with Blue Light

To determine the wavelengths in the preliminary experiments, blue and green lights were considered in the AVF pulse wave measurement. In the preliminary experiment, the pulse waves were observed with almost no inferiority in the cases of green and blue lights. This is because while general photoplethysmography (PPG) uses the absorbance change associated with arterial pulsation, this method is based on the displacement change between the camera and AVF associated with the pulsation of the AVF. Furthermore, the reason why we could not create a three-dimensional image with green light is that the anastomosis is a large blood vessel; therefore, it absorbs a large amount of light, and the reflected light decreases. On the other hand, as most of the blue light is reflected by the epidermis, the reflected light increases, and a three-dimensional image can be created.

5.3. Effect of Moving-Average Filtering

Finally, this paper discusses the reason why a curve exhibits a trend that a sharp increase in SNR before the peak of 0.179 ( N = 113 ) and gradual decrease after the peak (Figure 13). The pulse wave component can be assumed to be a periodic signal obtained from each pixel in the same phase. Moreover, the noise component can be assumed to be random noise caused by scattered reflection obtained from each pixel, where we can confirm the effect of moving-average filtering [31]. Figure 19 shows the left figure of Figure 13 in an expanded view with N = 65–201. The SNR slowly decreases from N = 113 , and sharply drops from N = 187 . As N increases, the waveform at the side of the anastomosis is averaged, and the reverse phase waveform is added to the average. In this study, the bandwidth of the AVF observation was set to 1–5 Hz; however, a few thrill waveforms were observed in this band; therefore, it is suggested that thrill waveforms exist in the frequency band of 5–10 Hz, which was excluded in this study. Therefore, the frequency range of the thrill waveforms can be used as a substitute for palpation.

5.4. Limitation of Developed Non-Contact AVF Evaluation System

Here, let us discuss the limitation of the developed system. In the developed system, the test results were strongly affected by ambient light because the detection of the luminance should include not only the change of blood pulse but also the ambient light. Hence, the developed system employs a light-shielding box to prevent the ambient light from being superimposed on the measurement data.
It is noted that, as for the authorization for the present method, careful discussion is needed to obtain CE marking or any market; however, the developed system does not require strict selection criteria for the patient or any special skills of doctors, so that, it can be widely spread to the market.

6. Conclusions

This study has developed a non-contact AVF measurement system for the visualization and quantification of hemodynamics in order to solve the problems of contact auscultation. In the preliminary experiments, luminance has been found to be a better choice to acquire the AVF pulse wave waveform than hue and saturation. Furthermore, in existing studies of AVF measurement based on optical technology, validation was rarely conducted based on clinical tests with a sufficient number of subjects. Therefore, this study has conducted clinical tests with 168 normal dialysis patients to validate the quantification of non-contact AVF measurements. From the clinical test results, the median SNRAVF value of the normalized distribution was −0.37, and the first (25 %tile), second (50 %tile), third (75 %tile), and fourth quartiles (the maximum value) were −0.77, −0.60, 0.23, and 1.44, respectively, which should be informative to construct criteria for the quantified evaluation of AVF measurement.
Currently, we are carrying out another clinical test to make measurements before and after stenosis to determine the hemodynamics of the stenosis site to realize early detection of abnormalities. In the future, we will apply and use the developed system to the AVF condition evaluation and for patients with peripheral circulatory failure, respectively.

Author Contributions

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

Funding

This work was supported by JSPS KAKENHI Grant Number 21H01325.

Institutional Review Board Statement

All subjects gave their informed consent for inclusion before they participated in the study. The protocol was approved by the Ethics Committee of Tamaki–Aozora Hospital, Japan on 2nd August 2021 (R3–1).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
AVFArteriovenous Fistula
BPFBand-Pass Filter
CCDCharge-Coupled Device
FFTFast Fourier Transform
PPGPhotoplethysmography
SNRSignal-to-Noise power Ratio

References

  1. Fresenius Medical Care Annual Report 2020. Available online: https://www.freseniusmedicalcare.com/fileadmin/data/com/pdf/Media_Center/Publications/Annual_Reports/FME_Annual_Report_2020_EN.pdf (accessed on 24 January 2022).
  2. Nitta, K.; Masakane, I.; Nakai, S.; Ogata, S.; Kimata, N.; Hanafusa, N.; Hamano, T. 2019 annual dialysis data report, JSDT renal data registry. J. Jpn. Soc. Dialysis Ther. 2020, 53, 579–632. [Google Scholar] [CrossRef]
  3. Corada, M.; Morini, M.F.; Dejana, E. Signaling pathways in the specification of arteries and veins. Arterioscler. Thrombosis Vasc. Biol. 2014, 34, 2372–2377. [Google Scholar] [CrossRef]
  4. Ene–Iordache, B.; Semperboni, C.; Dubini, G.; Remuzzi, A. Disturbed flow in a patient-specific arteriovenous fistula for hemodialysis: Multidirectional and reciprocating near–wall flow patterns. J. Biomech. 2015, 48, 2195–2200. [Google Scholar] [CrossRef]
  5. Remuzzi, A.; Bozzetto, M. Biological and physical actors involved in the maturation of arteriovenous fistula for hemodialysis. Cardiovasc. Eng. Technol. 2017, 8, 273–279. [Google Scholar] [CrossRef]
  6. Pablo, V. AV fistula diagnosis procedures. In Acoustical Signal Processing of Arterio Venous Fistula Bruits; Department of Electrical and Information Technology, Lund University: Lund, Sweden, 2015; pp. 10–11. [Google Scholar]
  7. 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]
  8. Panda, B.; Mandal, S.; Majerus, S. Vascular stenosis detection using temporal-spectral differences in correlated acoustic measurements. In Proceedings of the 2019 IEEE Signal Processing in Medicine and Biology Symposium, Philadelphia, PA, USA, 7 December 2019. [Google Scholar] [CrossRef]
  9. 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]
  10. Wang, H.; Wu, C.; Chen, C.; Lin, B. Novel noninvasive approach for detecting arteriovenous fistula stenosis. IEEE Trans. Biomed. Eng. 2014, 61, 1851–1857. [Google Scholar] [CrossRef]
  11. Chin, S.; Panda, B.; Damaser, M.S.; Majerus, S.J.A. Stenosis characterization and identification for dialysis vascular access. In Proceedings of the 2018 IEEE Signal Processing in Medicine and Biology Symposium, Philadelphia, PA, USA, 1 December 2018; pp. 1–5. [Google Scholar] [CrossRef]
  12. Gram, M.; Olesen, J.T.; Riis, H.C.; Selvaratnam, M.; Meyer-Hofmann, H.; Pedersen, B.B.; Schmidt, S.E. Stenosis detection algorithm for screening of arteriovenous fistulae. In Proceedings of the 15th Nordic-Baltic Conference on Biomedical Engineering and Medical Physics (NBC2011), Aalborg, Denmark, 14–17 June 2011; Springer: Berlin/Heidelberg, Fermany, 2011; pp. 241–244. [Google Scholar] [CrossRef]
  13. Ota, K.; Nishiura, Y.; Ishihara, S.; Adachi, H.; Yamamoto, Y.; Hamano, T. Evaluation of hemodialysis arteriovenous bruit by deep learning. Sensors 2020, 20, 4852. [Google Scholar] [CrossRef]
  14. Mansy, H.A.; Hoxie, S.J.; Patel, N.H.; Sandler, R.H. Computerised analysis of auscultatory sounds associated with vascular patency of haemodialysis access. Med. Biol. Eng. Comput. 2005, 43, 56–62. [Google Scholar] [CrossRef]
  15. 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]
  16. Yi–Chun, D.; 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]
  17. Chiang, P.Y.; Chao, P.C.P.; Tu, T.Y.; Kao, Y.H.; Yang, C.Y.; Tarng, D.C.; Wey, C.L. Machine learning classification for assessing the degree of stenosis and blood flow volume at arteriovenous fistulas of hemodialysis patients using a new photoplethysmography sensor device. Sensors 2019, 19, 3422. [Google Scholar] [CrossRef] [Green Version]
  18. Kamiyama, H.; Kitama, M.; Shimizu, H.O.; Yamashita, M.; Kojima, Y.; Shimizu, K. Fundamental study for optical transillumination imaging of arteriovenous fistula. Adv. Biomed. Eng. 2021, 10, 1–10. [Google Scholar] [CrossRef]
  19. Zhu, F.; Williams, S.; Putnam, H.; Campos, I.; Johnson, C.; Kappel, F.; Kotanko, P. Estimation of arterio-venous access blood flow in hemodialysis patients using video image processing technique. In Proceedings of the 38th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC2016), Orlando, FL, USA, 16–20 August 2016; pp. 207–210. [Google Scholar]
  20. Heckenlively, J.R.; Arden, G.B. Principles and Practice of Clinical Electrophysiology of Vision, 2nd ed.; MIT Press: Cambridge, MA, USA, 2006; p. 222. [Google Scholar]
  21. Kumar, A.; Bhattacharya, A.; Makhija, N. Evoked potential monitoring in anaesthesia and analgesia. Anaesthesia 2000, 55, 225–241. [Google Scholar] [CrossRef]
  22. Monroe, L. Physical Rehabilitation for the Physical Therapist Assistant, 1st ed.; Elsevier Health Sciences: Amsterdam, The Netherlands, 2014; p. 347. [Google Scholar]
  23. Shimazaki, T.; Hara, S.; Okuhata, H.; Nakamura, H.; Kawabata, T. Heart rate sensing during exercise by means of photoplethysmography. Trans. Jpn. Soc. Med. Biol. Eng. 2017, 54, 225–235. [Google Scholar]
  24. Logvinenko, A.D. The geometric structure of color. J. Vision Jan. 2015, 15, 16. [Google Scholar] [CrossRef] [PubMed]
  25. Spigulis, J.; Gailite, L.; Lihachev, A.; Erts, R. Simultaneous recording of skin blood pulsations at different vascular depths by multiwavelength Photoplethysmography. Appl. Opt. 2007, 46, 1754–1759. [Google Scholar] [CrossRef] [Green Version]
  26. Maeda, Y.; Sekine, M.; Tamura, T. Relationship between measurement site and motion artifacts in wearable reflected photoplethysmography. J. Med. Syst. 2011, 35, 969–976. [Google Scholar] [CrossRef]
  27. Shimazaki, T.; Hara, S.; Okuhata, H.; Nakamura, H.; Kawabata, T. Cancellation of motion artifact induced by exercise for PPG-based heart rate sensing. In Proceedings of the 36th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC2014), Chicago, IL, USA, 26–30 August 2014; pp. 3216–3219. [Google Scholar] [CrossRef]
  28. Simazaki, T.; Anzai, D.; Watanabe, K.; Nakajima, A.; Fukuda, M.; Ata, S. Heat stroke prevention in hot specific occupational environment enhanced by supervised machine learning with personalized vital signs. Sensors 2022, 22, 395. [Google Scholar] [CrossRef]
  29. Woodham, R.J. Photometric method for determining surface orientation from multiple images. Opt. Eng. 1980, 19, 139–144. [Google Scholar] [CrossRef]
  30. Santo, H.; Samejima, M.; Sugano, Y.; Shi, B.; Matsushita, Y. Deep photometric stereo networks for determining surface normal and reflectances. IEEE Trans. Pattern Anal. Mach. Intell. 2022, 44, 114–128. [Google Scholar] [CrossRef] [PubMed]
  31. Nakamura, M.; Nishida, S.; Shibasaki, H. Spectral properties of signal averaging and a novel technique for improving the signal-to-noise ratio. J. Biomed. Eng. 1989, 11, 72–78. [Google Scholar] [CrossRef]
Figure 1. Overview of hemodialsis therapy.
Figure 1. Overview of hemodialsis therapy.
Sensors 22 02745 g001
Figure 2. Experimental equipment. (a) External view of the experimental equipment. (b) Internal view of the light-shielding box.
Figure 2. Experimental equipment. (a) External view of the experimental equipment. (b) Internal view of the light-shielding box.
Sensors 22 02745 g002
Figure 3. Setup and analysis flow of the non-contact AVF measurement system.
Figure 3. Setup and analysis flow of the non-contact AVF measurement system.
Sensors 22 02745 g003
Figure 4. Correspondence between brightness and color tone for the color mapping.
Figure 4. Correspondence between brightness and color tone for the color mapping.
Sensors 22 02745 g004
Figure 5. AVF imaging area and the kernel size ( N × N ) of the moving-average filter.
Figure 5. AVF imaging area and the kernel size ( N × N ) of the moving-average filter.
Sensors 22 02745 g005
Figure 6. Frequency spectrum of AVF pulse wave waveform.
Figure 6. Frequency spectrum of AVF pulse wave waveform.
Sensors 22 02745 g006
Figure 7. Pulse wave measurement based on changes in (a) luminance, (b) hue, and (c) saturation. The measurement was conducted on the palm of a healthy person’s hand. The horizontal axis shows the data sample index with time resolution of 25 ms (the sampling frequency of 40 Hz). It was confirmed that the best way to measure the AVF pulse wave was achieved based on the change with the luminance.
Figure 7. Pulse wave measurement based on changes in (a) luminance, (b) hue, and (c) saturation. The measurement was conducted on the palm of a healthy person’s hand. The horizontal axis shows the data sample index with time resolution of 25 ms (the sampling frequency of 40 Hz). It was confirmed that the best way to measure the AVF pulse wave was achieved based on the change with the luminance.
Sensors 22 02745 g007
Figure 8. AVF pulse wave measured with blue light. (a) Photograph of a left wrist. (b) AVF pulse wave measured at the two points close to the anastomosis. The low- and high-cutoff frequencies of band pass filter were 1 Hz and 10 Hz, respectively. Because the similar waveforms are observed at different points, the measured data is not sensitive to observation points. (c) Three-dimensional image by the photometric method with blue light. As compared with the photo shown in (a), good agreement is confirmed, which means that the three-dimensional image based on the luminance is efficient to represent the AVF form.
Figure 8. AVF pulse wave measured with blue light. (a) Photograph of a left wrist. (b) AVF pulse wave measured at the two points close to the anastomosis. The low- and high-cutoff frequencies of band pass filter were 1 Hz and 10 Hz, respectively. Because the similar waveforms are observed at different points, the measured data is not sensitive to observation points. (c) Three-dimensional image by the photometric method with blue light. As compared with the photo shown in (a), good agreement is confirmed, which means that the three-dimensional image based on the luminance is efficient to represent the AVF form.
Sensors 22 02745 g008
Figure 9. Setup for comparison between contact and non-contact measurements.
Figure 9. Setup for comparison between contact and non-contact measurements.
Sensors 22 02745 g009
Figure 10. Comparison with reference data. (a) Reference data obtained by bidirectional Doppler flowmeter. (b) Measurement data by the developed system with blue light source at radial artery. Since the observed pulse wave was inverted in the developed system, both periods were compared in the minimum values. It was confirmed that the period of the peak value of the blood flow meter (reference) and the peak value of the radial artery were in agreement.
Figure 10. Comparison with reference data. (a) Reference data obtained by bidirectional Doppler flowmeter. (b) Measurement data by the developed system with blue light source at radial artery. Since the observed pulse wave was inverted in the developed system, both periods were compared in the minimum values. It was confirmed that the period of the peak value of the blood flow meter (reference) and the peak value of the radial artery were in agreement.
Sensors 22 02745 g010
Figure 11. Continuous color -mapping image of AVF pulse wave for one cycle (2.600 to 3.360 s, the band pass filter was applied). (a) Anastomotic tense phase (2.820 s). (b) Transition from anastomosis to puncture site (3.000 s). (c) Peak puncture site phase (3.125 s). (d) Anastomotic relaxation phase (3.150 s). (e) End of anastomotic relaxation phase (3.325 s).
Figure 11. Continuous color -mapping image of AVF pulse wave for one cycle (2.600 to 3.360 s, the band pass filter was applied). (a) Anastomotic tense phase (2.820 s). (b) Transition from anastomosis to puncture site (3.000 s). (c) Peak puncture site phase (3.125 s). (d) Anastomotic relaxation phase (3.150 s). (e) End of anastomotic relaxation phase (3.325 s).
Sensors 22 02745 g011
Figure 12. Effect of the kernel size N on noise reduction in the AVF pulse wave measurement. The random noise on the AVF pulse wave was suppressed by the moving-average filtering. Larger kernel size made the observed AVF pulse wave clearer.
Figure 12. Effect of the kernel size N on noise reduction in the AVF pulse wave measurement. The random noise on the AVF pulse wave was suppressed by the moving-average filtering. Larger kernel size made the observed AVF pulse wave clearer.
Sensors 22 02745 g012
Figure 13. Determination of optimal kernel size N. S N A V F = 0.179 was achieved when the optimal kernel N = 113 .
Figure 13. Determination of optimal kernel size N. S N A V F = 0.179 was achieved when the optimal kernel N = 113 .
Sensors 22 02745 g013
Figure 14. Typical examples of AVF pulse wave measurement with optimal kernel N = 113 .
Figure 14. Typical examples of AVF pulse wave measurement with optimal kernel N = 113 .
Sensors 22 02745 g014
Figure 15. Average and variance of measured luminance for all patiants with optimal kernel N = 113 .
Figure 15. Average and variance of measured luminance for all patiants with optimal kernel N = 113 .
Sensors 22 02745 g015
Figure 16. Measurement example of reversed-phase AVF pulse wave. (a) Photograph of measurement points on AVF. (b) Measured waveforms obtained at the two measurement points.
Figure 16. Measurement example of reversed-phase AVF pulse wave. (a) Photograph of measurement points on AVF. (b) Measured waveforms obtained at the two measurement points.
Sensors 22 02745 g016
Figure 17. Detailed demonstration of the reversed-phase AVF pulse waveform. The luminance of the (01) measurement point at 0.4 s increased, whereas that of the (02) point decreased. On the other hand, at 0.9 s, the luminance of the (02) point increased, and that of the (01) point inversely decreased. It was suggested that the shadows caused by the masses may have an effect on the phenomenon of the reversed-phase AVF pulse waveform.
Figure 17. Detailed demonstration of the reversed-phase AVF pulse waveform. The luminance of the (01) measurement point at 0.4 s increased, whereas that of the (02) point decreased. On the other hand, at 0.9 s, the luminance of the (02) point increased, and that of the (01) point inversely decreased. It was suggested that the shadows caused by the masses may have an effect on the phenomenon of the reversed-phase AVF pulse waveform.
Sensors 22 02745 g017
Figure 18. Mechanism of blood flow in cephalic vein.
Figure 18. Mechanism of blood flow in cephalic vein.
Sensors 22 02745 g018
Figure 19. Expanded view of the SNR variation with the kernel size N = 65–201. The reason for the slow decrease in SNR from S N AVF = 0.179 ( N = 113 ) can be explained by the fact that the reversed-phase waveform at the side of the anastomosis is also averaged as the kernel size N increases.
Figure 19. Expanded view of the SNR variation with the kernel size N = 65–201. The reason for the slow decrease in SNR from S N AVF = 0.179 ( N = 113 ) can be explained by the fact that the reversed-phase waveform at the side of the anastomosis is also averaged as the kernel size N increases.
Sensors 22 02745 g019
Table 1. Information and protocol of clinic test.
Table 1. Information and protocol of clinic test.
DateFrom 12 October to 5 November 2021
PlaceTamaki-Aozora Hospital, Tokushima, Japan
Number of patients168
Patient selection criteriaPatient who has a normal shunt condition with a blood flow of more than 200 mL/min
PreparationPurpose and rule of the test were sufficiently explained to all patients beforehand
Rules for rejection of testsThey can refuse to participate at any time during the test
Measurement protocolMeasurement in resting sitting condition for 10 s
Table 2. Statics of patients.
Table 2. Statics of patients.
Average age69.1 years old
Average dialysis years8.2 years
Average blood flow226.8 mL/min
Standard deviation of blood flow28.0 mL/m
Table 3. Hardware specifications.
Table 3. Hardware specifications.
PartsModel NumberVendersFeatures
Image processing softwareHalcon Rev.18MVTec
CMOS cameraa2A1920–160 μmBASLER1920 × 1200 pixels, 160 fps
CCTV lensFA0802DCHIOPTF#1.4–16
Polar screen#52–556Edmund Optics
IlluminatorIMAR–130DB–8chLEIMACCenter wavelength: 465 nm
Polarizer#45–204LEIMACMounted in the front of a illuminator and CCTV Lens
ControllerIDGB–30M8PG–TPLEIMACAC 100–240 V, DC 12 V, 30 W
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

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. https://doi.org/10.3390/s22072745

AMA Style

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(7):2745. https://doi.org/10.3390/s22072745

Chicago/Turabian Style

Iwai, Rumi, Takunori Shimazaki, Yoshifumi Kawakubo, Kei Fukami, Shingo Ata, Takeshi Yokoyama, Takashi Hitosugi, Aki Otsuka, Hiroyuki Hayashi, Masanobu Tsurumoto, and et al. 2022. "Quantification and Visualization of Reliable Hemodynamics Evaluation Based on Non-Contact Arteriovenous Fistula Measurement" Sensors 22, no. 7: 2745. https://doi.org/10.3390/s22072745

APA Style

Iwai, R., Shimazaki, T., Kawakubo, Y., Fukami, K., Ata, S., Yokoyama, T., Hitosugi, T., Otsuka, A., Hayashi, H., Tsurumoto, M., Yokoyama, R., Yoshida, T., Hirono, S., & Anzai, D. (2022). Quantification and Visualization of Reliable Hemodynamics Evaluation Based on Non-Contact Arteriovenous Fistula Measurement. Sensors, 22(7), 2745. https://doi.org/10.3390/s22072745

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