# Precise Detection of Wrist Pulse Using Digital Speckle Pattern Interferometry

Authors: Hengfei Zhang; Sijin Wu; Weixian Li; Yonghong Wang; Mingli Dong; Lianxiang Yang
Journal: Evidence-Based Complementary and Alternative Medicine
 (2018)

---

## Abstract

Pulse diagnosis is one of the four diagnostic methods of traditional Chinese medicine. However it suffers from the lack of objective and efficient detection method. We propose a noncontact optical method to detect human wrist pulse, aiming at the precise determination of the temporal and spatial distributions of pulse. The method uses the spatial-carrier digital speckle pattern interferometry (DSPI) to measure the micro/nanoscale skin displacement dynamically. Significant improvements in DSPI measurement have been made to allow the DSPI to detect the comprehensive information of the arterial pulsation at locations of Cun, Guan, and Chi. The experimental results prove that the spatiotemporal distributions of pulse can be obtained by the proposed method. The obtained data can be further used to describe most of the pulse parameters such as rate, rhythm, depth, length, width, and contour.


---

## Body

## 1. Introduction

Human arterial pulsation is considered as a vital sign. It includes much information of human health. In Chinese traditional medicine (TCM), the arterial pulsation at the locations of Cun, Guan, and Chi, which are near wrist joint, as shown in Figure1, is used for health evaluation and disease diagnosis. The method of wrist pulse investigation is called pulse diagnosis, one of the four diagnostic methods in TCM. The pulse diagnosis has played an important role during the past one thousand years in East Asia and also been a practical diagnostic method in modern society [1–3]. Many examples of health evaluation and disease diagnosis using the pulse diagnosis have been introduced [4–6]. The pulse diagnosis is also a practical method in other oriental medicines, such as Indian Ayurveda [7].Figure 1
Pulsation at the locations of Cun, Guan, and Chi.Traditionally, physicians of Chinese medicine touch patient’s wrist and feel the pulsations at the three locations: Cun, Guan, and Chi. The diagnostic test relies entirely on the doctor’s own experiential judgment. This subjective test may be affected by many factors. Moreover, the diagnostic results are hard to record and exchange. Thus, the quantification and standardization of the pulse diagnosis are very important. It requires the transformation of doctors’ subjective feelings into objective physical quantities. To achieve this goal, the precise determination of the temporal and spatial distributions of the wrist pulse is required.Owning to the development of modern transducers and computer technology, many physical methods have been presented to detect the wrist pulse during the last thirty years. These methods are mainly classified into pressure sensing and optical test methods. The pressure sensing methods are further divided into piezoresistive, piezoelectric, and piezocapacitive sensors. These sensors are attached to wrist skin tightly in order to sense the pulse pressure. For example, Sook Hyang Yoon et al. attached three semiconductor resistance strain gauge piezoresistive sensors at the three locations of Cun, Guan, and Chi to detect their pressure changes and pulse waves [8]. Kalange et al. used a single piezoelectric transducer to obtain the ideal pulse shape [7]. Hu developed a piezocapacitive sensor array that enabled multipoint detection of pulse [9]. Piezoelectric sensor based on polyvinylidene fluoride (PVDF) is a promising sensor for pulse detection due to its advantages of thin and soft structure, matching well with the acoustic impedance of human’s soft tissue, and ease to achieve multipoint measurement [10]. Another promising material for pressure sensor is newly developed extraordinarily stretchable all carbon collaborative nanoarchitectures which can be used as an epidermal sensor [11]. The pressure sensing method is the dominating method in human pulse detection because of its highly sensitive detection and ease of use. However, the insufficient information of its test results makes it an unideal pulse detection method. Only a single or several points of artery are tested. The important information of the spatial distribution such as pulse length and width is lost. Moreover, though it is sensitive to arterial pulsation, it is hard to convert the outputs of the pressure transducers to the depth of wrist pulse directly. Consequently, the pressure sensing method has not been an ideal method of pulse detection yet.Optical method can be divided into point measurement and full-field measurement methods. These methods detect the pulse by measuring the displacements of the human skin at selected locations or the full field. The optical point measurement, such as optical interferometry [12], laser triangulation [13], and optical fiber measurement [14], suffers from the same problem of the pressure sensing method that cannot get the spatial distribution of the wrist pulse.Full-field measurement using optical method is desired to allow the obtainment of the pulse spatial distributions. Efforts in this research area mainly focus on vision measurement and digital image correlation (DIC) techniques. For example, Zhang developed an image device to satisfy the full-field detection of pulse [15]. Shao and Dai proposed a DIC method to realize the real-time and full-field measurement of pulse [16]. Xue and Su subsequently improve the pulse detection method using DIC and exhibit much reasonable results [17]. Both vision measurement and DIC can get the temporal and spatial distributions of wrist pulse, allowing much comprehensive information of wrist pulse being obtained. However, the measuring resolution is limited to several micrometers or even worse by their measuring principle. This is not enough to precisely determine the pulse amplitude at any locations. For example, the maximum pulse amplitudes at some locations beside Cun, Guan, and Chi are always less than several micrometers. It means the pulses at Cun, Guan, and Chi can be detected by these methods, but the pulses at other locations may not be successfully detected. In addition, additional cooperation target such as grid line and water transfer printing speckle patterns on the skin surface is required when vision measurement or DIC is used, reducing the convenience of pulse detection.We introduce a precise pulse detection method based on digital speckle pattern interferometry (DSPI), allowing the highly sensitive detection of the temporal and spatial distributions of the pulse and the obtainment of more comprehensive and detailed pulse information. DSPI is a full-field optical technique which can measure the displacement/deformation distributions of object surfaces. Compared to vision measurement and DIC, DSPI enjoys significant advantages of much higher measuring accuracy and independence from cooperation target [18, 19]. As an optical interference technique, its resolution reaches 20 nm, which can trace back to the laser wavelength and is unsusceptible to device and external factors. Thus, slight change of the pulse can be detected by DSPI easily. We also improve the DSPI in data sampling, phase unwrapping, and data filtering, to achieve the in vivo detection of the human wrist pulse.


## 2. Methods

### 2.1. Principle and Experimental Setup of Spatial-Carrier DSPI

DSPI is classified, based on their optical setups and phase determination methods, into phase-shifting DSPI (PS-DSPI) and spatial-carrier DSPI (SC-DSPI) [20]. The SC-DSPI is also known as digital holographic interferometry [21]. The SC-DSPI outperforms the PS- DSPI in measuring speed, so it is a more competitive candidate for pulse detection. The optical arrangement of the SC-DSPI is depicted by Figure 2. The laser beam is divided into two beams, the object beam and reference beam, by an adjustable splitter. The object beam illuminates the vibrating object surface via a beam expander along the illumination directionKi. The scattered object beam is collected by an imaging lens along the observation directionKo and then forms an image on the image sensor of the camera via an aperture. The reference beam is coupled into a fiber and then strikes the image sensor at a small angle between it and the optical axis. The object and reference beams encounter each other on the image sensor, yielding optical speckle interference.Figure 2
Optical setup for wrist pulse detection.The intensity distribution of the interferogram captured by the camera is given by(1)Ix,y=Ox,y2+Rx,y2+Ox,yRx,y∗+Rx,yOx,y∗,where O and R are the object and reference beams, respectively, (x,y) is the two-dimensional (2D) spatial coordinates on the image plane, and ∗ is the operator of complex conjugate.The inclination angle of the reference beam towards the optical axis determines the spatial frequency of the speckle interferogram. The maximal spatial frequency is calculated by(2)fmax=λ2sin⁡θmax2,where λ is the wavelength of the laser used in the experiment and θmax is the maximal inclination angle of the reference beam. According to the Nyquist-Shannon sampling theorem, the experimental setup needs to meet the condition 2fmax<1/p, wherep is the size of the camera pixel. This condition limits the range of the reference beam inclination angle. Details about the determination of the reference beam’s deflection can be found in [22].


### 2.2. Collection of Speckle Interferograms and Phase Reconstruction

The last two terms of (1) contain the information of the object beam’s phase. These terms are separated from the background in the frequency domain by performing a 2D fast Fourier Transform (FFT) on the recorded interferogram I(x,y). Subsequently the complex amplitude of the object beam’s wavefront is obtained after performing an inverse fast Fourier Transform (IFFT) on the selected frequency component. Due to the sensor discretization, the image captured by the camera composes m×n pixels whose sizes are Δx×Δy, where m and n are integer numbers. The speckle interferogram’s intensity can be rewritten as O(mΔx,nΔy). The phase of the wavefront can be calculated by (3)ΦwmΔx,nΔy=arctan⁡Im⁡OmΔx,nΔyRe⁡OmΔx,nΔy,where the subscript w denotes a wrapped phase and Im and Re indicate the imaginary and real parts, respectively.DSPI records a sequence of speckle interferograms and determines the phaseΦw(mΔx,nΔy,iΔτ) at the point (mΔx,nΔy) and the timeiΔτ, where Δτ is the time interval between two camera exposures. The phase difference between two adjacent phase distributions according to time (i-1) Δτ andiΔτ can be calculated by(4)ΔΦi,i-1=ΦwmΔx,nΔy,iΔτ-ΦwmΔx,nΔy,i-1Δτ,After the phase differences are obtained, a sine/cosine average filter is used to filter the noise at the obtained phase maps, yielding smoothed phase maps. The entire process of digital speckle interferogram collection and phase reconstruction is shown in Figure3.Figure 3
Process of the phase difference map calculation.The relationship between the phase difference and displacementd is given by(5)ΔΦ=2πλds,where s=Ki-Ko is the sensitivity vector which is determined by the geometric layout of the measuring setup.Usually the angle between the illumination and observation directions is very small. The absolute value of the sensitivity vectorsis almost 2. Consequently, the out-of-plane deformation is calculated by(6)wi,i-1=λ4πΔΦi,i-1,where wi,i-1 is the out-of-plane deformation during the time interval.


### 2.3. Phase Unwrapping

Phase unwrapping (PU) is a process to remove the 2π discontinuities within the phase image. It detects a 2π phase jump and adds or subtracts integer multiple of 2π to successive pixels following that phase jump based on a threshold mechanism. The existing algorithms can be classified into three categories, that is, temporal phase unwrapping (TPU), spatial phase unwrapping (SPU), and spatiotemporal three-dimensional phase unwrapping (STPU) [23].TPU is carried out along the time axis due to the independence of each camera pixel. The phase at each pixel is treated as a function of time and unwrapped by a one-dimensional unwrapping algorithm. SPU involves a spatial comparison of phase values at neighboring pixels, which can be further classified into path-dependent and path-independent methods. STPU unwraps the whole phase map along the 2D spatial domain and several selected points along the time axis. The combination of the spatial and temporal phase distributions yields the dynamic phase difference maps. The operating procedure of STPU is depicted by Figure4.Figure 4
Procedure of spatiotemporal three-dimensional phase unwrapping.


## 2.1. Principle and Experimental Setup of Spatial-Carrier DSPI

DSPI is classified, based on their optical setups and phase determination methods, into phase-shifting DSPI (PS-DSPI) and spatial-carrier DSPI (SC-DSPI) [20]. The SC-DSPI is also known as digital holographic interferometry [21]. The SC-DSPI outperforms the PS- DSPI in measuring speed, so it is a more competitive candidate for pulse detection. The optical arrangement of the SC-DSPI is depicted by Figure 2. The laser beam is divided into two beams, the object beam and reference beam, by an adjustable splitter. The object beam illuminates the vibrating object surface via a beam expander along the illumination directionKi. The scattered object beam is collected by an imaging lens along the observation directionKo and then forms an image on the image sensor of the camera via an aperture. The reference beam is coupled into a fiber and then strikes the image sensor at a small angle between it and the optical axis. The object and reference beams encounter each other on the image sensor, yielding optical speckle interference.Figure 2
Optical setup for wrist pulse detection.The intensity distribution of the interferogram captured by the camera is given by(1)Ix,y=Ox,y2+Rx,y2+Ox,yRx,y∗+Rx,yOx,y∗,where O and R are the object and reference beams, respectively, (x,y) is the two-dimensional (2D) spatial coordinates on the image plane, and ∗ is the operator of complex conjugate.The inclination angle of the reference beam towards the optical axis determines the spatial frequency of the speckle interferogram. The maximal spatial frequency is calculated by(2)fmax=λ2sin⁡θmax2,where λ is the wavelength of the laser used in the experiment and θmax is the maximal inclination angle of the reference beam. According to the Nyquist-Shannon sampling theorem, the experimental setup needs to meet the condition 2fmax<1/p, wherep is the size of the camera pixel. This condition limits the range of the reference beam inclination angle. Details about the determination of the reference beam’s deflection can be found in [22].


## 2.2. Collection of Speckle Interferograms and Phase Reconstruction

The last two terms of (1) contain the information of the object beam’s phase. These terms are separated from the background in the frequency domain by performing a 2D fast Fourier Transform (FFT) on the recorded interferogram I(x,y). Subsequently the complex amplitude of the object beam’s wavefront is obtained after performing an inverse fast Fourier Transform (IFFT) on the selected frequency component. Due to the sensor discretization, the image captured by the camera composes m×n pixels whose sizes are Δx×Δy, where m and n are integer numbers. The speckle interferogram’s intensity can be rewritten as O(mΔx,nΔy). The phase of the wavefront can be calculated by (3)ΦwmΔx,nΔy=arctan⁡Im⁡OmΔx,nΔyRe⁡OmΔx,nΔy,where the subscript w denotes a wrapped phase and Im and Re indicate the imaginary and real parts, respectively.DSPI records a sequence of speckle interferograms and determines the phaseΦw(mΔx,nΔy,iΔτ) at the point (mΔx,nΔy) and the timeiΔτ, where Δτ is the time interval between two camera exposures. The phase difference between two adjacent phase distributions according to time (i-1) Δτ andiΔτ can be calculated by(4)ΔΦi,i-1=ΦwmΔx,nΔy,iΔτ-ΦwmΔx,nΔy,i-1Δτ,After the phase differences are obtained, a sine/cosine average filter is used to filter the noise at the obtained phase maps, yielding smoothed phase maps. The entire process of digital speckle interferogram collection and phase reconstruction is shown in Figure3.Figure 3
Process of the phase difference map calculation.The relationship between the phase difference and displacementd is given by(5)ΔΦ=2πλds,where s=Ki-Ko is the sensitivity vector which is determined by the geometric layout of the measuring setup.Usually the angle between the illumination and observation directions is very small. The absolute value of the sensitivity vectorsis almost 2. Consequently, the out-of-plane deformation is calculated by(6)wi,i-1=λ4πΔΦi,i-1,where wi,i-1 is the out-of-plane deformation during the time interval.


## 2.3. Phase Unwrapping

Phase unwrapping (PU) is a process to remove the 2π discontinuities within the phase image. It detects a 2π phase jump and adds or subtracts integer multiple of 2π to successive pixels following that phase jump based on a threshold mechanism. The existing algorithms can be classified into three categories, that is, temporal phase unwrapping (TPU), spatial phase unwrapping (SPU), and spatiotemporal three-dimensional phase unwrapping (STPU) [23].TPU is carried out along the time axis due to the independence of each camera pixel. The phase at each pixel is treated as a function of time and unwrapped by a one-dimensional unwrapping algorithm. SPU involves a spatial comparison of phase values at neighboring pixels, which can be further classified into path-dependent and path-independent methods. STPU unwraps the whole phase map along the 2D spatial domain and several selected points along the time axis. The combination of the spatial and temporal phase distributions yields the dynamic phase difference maps. The operating procedure of STPU is depicted by Figure4.Figure 4
Procedure of spatiotemporal three-dimensional phase unwrapping.


## 3. Experiment and Result Discussion

A 5 W diode pumped solid state laser with a central wavelength of 532nm (Coherent Verdi G5) and CMOS camera (CatchBEST MU3S230 (SGYYO)) with a frame rate of 80 frames per second (fps) at 2.3 megapixels were used in the experiment. The frame rate can be enhanced when the pixel resolution is reduced. The camera was triggered by an external clock signal which guaranteed the uniformly sampling of the wrist pulse. A nonspherical lens whose focal length is 100 mm was used for imaging.Reference [16] shows that the maximum pulse displacement velocity is about 10 μm/s. This velocity is transferred to the phase change rate of about 37.6 cycles per second according to (6). In order to satisfy the Nyquist-Shannon sampling theorem, at least 2 points per phase cycle should be sampled. Actually, at least 4 sample points per phase cycles should be provided to avoid the influence of phase noise. Thus, a sample rate of 150.4 Hz is necessary to satisfy the wrist pulse detection. In the experiment, the frame rate was set as 160 fps which is faster than the minimum requirement of sample rate, guaranteeing the effective sampling. This frame rate was realized by reducing the image region of interesting to 1500×500 pixels, corresponding to a measurement area of 75 mm × 25 mm.The right wrist of a volunteer was fixed by two lashing bands on a backplane. In addition, the output of the laser was adjusted to 0.25 W to avoid human’s discomfort. In the experiment, totally 28800 digital speckle interferograms were captured for further processing. The vibration of pulse during the test time of 3 minutes was then obtained. The experimental setup is shown in Figure5.Figure 5
Experimental setup for wrist pulse detection.Results within a period of 30 s were selected for exhibition. Figure6(a) illustrates the absolute displacement variation at the location of pixel (290, 480). The complex and irregular data indicate the influence of a composite factor that includes the rigid body motion of human wrist and random displacement fluctuation of human skin. The quality of the result was much improved by applying a 0.7 Hz to 2 Hz band pass filtering towards the amplitude spectrum which is obtained by performing a FFT on the raw data. The amplitude spectrum and denoised data are shown in Figures 6(b) and 6(c), respectively.Figure 6
Absolute displacement of pulse at the location of pixel (290, 480): (a) raw data obtained by performing the temporal phase unwrapping; (b) amplitude spectrum; (c) raw and denoised data.
(a) (b) (c)The displacement detail and frequency spectrum of the pulse vibration are depicted by Figure7. Figure 7(a) shows the peak to peak value of the pulse amplitude is approximately 9 μm at the selected location. Figure 7(b) shows that the pulse frequency is 0.933 Hz which is equivalent to 60 pulses per minute.Figure 7
Detail of the pulse amplitude and frequency: (a) pulse displacement; (b) pulse amplitude spectrum.
(a) (b)The spatial distributions of pulse vibration are given by Figure8. These displacement maps correspond to different time of the pulse vibration.Figure 8
Spatial displacement distributions of the pulse.
(a) (b) (c) (d) (e) (f) (g) (h) (i) (j)


## 4. Conclusions

The paper presents a human wrist arterial pulsation detection method using SC-DSPI. The feasibility of pulse detection using the presented method is illustrated by the theoretical analysis as well as the experiment. The advantage of noninvasive, full-field, and high-precision detection has also been demonstrated.



---

*Source: 4187349-2018-06-07.xml*