Research Article  Open Access
Xingguang Geng, Su Liu, Yitao Zhang, Jiena Hou, Shaolong Zhang, Jun Zhang, Haiying Zhang, "A Noncontact Method for Locating Radial Artery above Radial Styloid Process in Thermal Image", EvidenceBased Complementary and Alternative Medicine, vol. 2020, Article ID 4057154, 9 pages, 2020. https://doi.org/10.1155/2020/4057154
A Noncontact Method for Locating Radial Artery above Radial Styloid Process in Thermal Image
Abstract
A radial artery above the radial styloid process is called GUAN and is a critical position for collecting pulse wave in traditional Chinese medicine theory. Locating GUAN is a precondition for collecting radial pulse wave. However, existing methods for locating GUAN lead to large deviations. This paper proposes a novel nontouch method for locating GUAN based on thermal imaging and image processing. This method consists of three parts: the infrared thermal imaging location imaging platform, the wrist edge contour extraction algorithm based on arbitrary angle edge recognition, and radial protrusion recognition algorithm (x coordinate identification algorithm of GUAN) and radial artery fitting algorithm (y coordinate identification algorithm of GUAN). The infrared thermal imaging positioning imaging platform is used to ensure that the wrist of the subject enters the fixed imaging area in a fixed position during each measurement and transmits the thermal imaging images carrying the image information of radial processes and radial arteries to the upper computer. Arbitrary angle edge recognition algorithm is used to extract wrist contour and radial artery edge information. The xaxis coordinates of the radial artery were provided by the identification algorithm, and the yaxis coordinates of the radial artery were provided by the fitting algorithm. Finally, the x and y coordinates determine the GUAN position. The algorithm for locating GUAN could provide repeatable and reliable x and y coordinates. The proposed method shows that relative standard deviation (RSD) of x distance of GUAN is less than 9.0% and RSD of y distance of GUAN is less than 5.0%. The proposed method could provide valid GUAN coordinates and reduce deviations of locating GUAN.
1. Introduction
Radial pulse waves contain rich information about the physiological condition of a body. The analysis of the radial pulse waves is regarded as an important approach to evaluate patient status in traditional Chinese medicine (TCM) [1–3]. GUAN is a segment of the radial artery above the radial styloid process and a critical position for collecting pulse wave correctly in TCM diagnosis [4–8]. However, current methods for locating GUAN, including tactile sense and pressure sensor array, could not achieve precise positioning [5, 6, 9–13]. When sizes of those single sensors are almost less than 10 mm × 10 mm, it is difficult for people to place sensitive areas of a single sensor at the center of GUAN each time by tactile sense [4, 7, 13–15]. When sizes of those single sensors are large, those sensors tend to cover CUN and CHI, which are positions near GUAN [4, 7, 13–15]. Besides, locating GUAN is not achieved by the pressure sensor arrays with a small gap (less than 1 mm) [6, 10] because sensor units in the array could be susceptible to vibration interference among them through connection materials.
To improve position accuracy and repeatability, a nontouch method based on thermal imagery and image processing for locating GUAN is proposed in this paper. This method contains two parts: a specially made thermal imaging platform and an algorithm for recognizing GUAN.
Human skin has a large absorption coefficient in the visible and near UV spectrum. The images of the radial artery in deep forearm tissue are not obtained in this spectrum [16, 17]. A thermal camera provides vital information for the nontouch method. It could not only generate an image of a part of the radial artery but also reflect the wrist outline by detecting temperature differences among skin above the radial artery, skin above other tissues, and background [18].
The algorithm for locating GUAN mainly recognizes characteristics of the radial styloid process and radial artery [19–21], and coordinates of GUAN are generated according to these characteristics and TCM theory.
In this paper, we propose a novel method for locating GUAN based on thermal imaging and image processing. The method includes three steps: (1) getting the thermal image of a wrist through the proposed thermal imaging platform; (2) extracting continuous wrist edges by edge detection algorithm; (3) obtaining x and y coordinates of GUAN by recognizing features of the radial styloid process and radial artery. The proposed method could achieve better positioning accuracy and provide reliable coordinates of GUAN.
2. Materials and Methods
2.1. Thermal Imaging Platform
As shown in Figure 1, the proposed platform is applied to keep the wrist of the subject firmly and generate a thermal image. The platform mainly consists of the following parts: a thermal imager (MAG32, 8 to 14 μm, Magnity Electronics Co. Ltd.), wrist support, and a line laser. Before imaging, the wrist crease of the subject should be aligned with a line laser. It could ensure that the wrist is placed in roughly the same position and make sure that the wrist crease is at the edge of the image. Then, the subject holds a handle tightly on the wrist support to keep a unified posture and extrude wrist joint for imaging. Finally, the thermal imager collects image information of the radial artery and the radial styloid process.
2.2. Edge Detection Method
The quality of edge detection is related to the detection angle, scale, and threshold. Liu et al. propose a novel edge detection algorithm to detect edges at arbitrary angles, and the optimized number of angles in the algorithm is introduced [22]. This algorithm combines singularity detection with Gaussian wavelet transform and edge detection at arbitrary directions and contains three steps:(1)A compact connection and a sparse connection are used to achieve the pixel combination principle of a corner.(2)According to the above pixel line combination principle, the image is segmented into several pixel lines along the angle, and the local maximum value of the pixel line is identified and the nonmaximum value is suppressed. The maximum point of each pixel line is recombined into an image edge along a certain angle. According to the above steps, the image edges generated by the 0° to 45° detection angle are superimposed to generate multiangle image edges.(3)According to the angle expansion rule, the detection angle from 0° to 45° is extended to 0° to 360° to generate the final image edge.
2.3. Locating GUAN
When a subject stretches his/her wrist and holds a handle on the wrist support of the proposed platform, the radial styloid process and wrist joint appear in the image (Figure 2(a)). The edges of the image (Figure 2(b)) could be obtained using an edge detection algorithm. The radial artery and the radial styloid process could be recognized based on the edges Figure 3.
(a)
(b)
(c)
(d)
(e)
(f)
(a)
(b)
(c)
(d)
2.3.1. Locating X Coordinate of GUAN
The wrist edge could be extracted from other edges, and the radial styloid process could be recognized based on the wrist edge and its characteristics. Then, the algorithm of locating x coordinate of GUAN contains two parts: extracting wrist edge and locating x coordinate.(1)Extracting wrist edge: this preprocessing includes identifying the maximum connected domain of the wrist edge and cutting branches: If the maximum connected domain reaches the boundary of the left and right sides of the image, that is, there is no break point in the connected domain, the maximum connected domain can be considered as the edge of the arm wrist. When the maximum connected domain does not reach the left and right edges of the image, it is necessary to connect the breaking point of the arm edge to form an overall edge of the arm wrist through the left and right edges of the image. The left break point of the maximum connected domain is taken as the origin to search for the edge fragment within the range of 2 pixels in the five directions of the point: up, left up, left, left down, and down. If there are other connected domains in the search area, the two connected domains will be connected, and the intermediate breakpoint pixel will supplement pixels between the two fragments by interpolation and finally form a new connected domain. Moreover, the breakpoint on the left side of the new connected domain will be taken as the origin to further search for other edge fragments until it reaches the left edge of the image. When the left edge of the maximum connected domain of the wrist is connected, the new right break point of the connected domain is taken as the origin to search and connect the edge fragments within the range of 2 pixels in the five directions of the point, up, left up, left, left down, and down, and form a new connected domain until it reaches the right edge. Figure 3(a) shows the matrix diagram of the connecting operator. And Figure 3(b) shows the effect drawing of the connecting operator. A part of the wrist edge is extracted by identifying the maximal connected component. Then, branches of are cut by cutting operator along the longaxis direction of . The cutting operator is chosen by a longaxis direction. When the long axis is along the xaxis, the horizontal cutting operator and the vertical cutting operator will cut the branches on the top left, the bottom left, the top, the bottom, the top right, and the bottom right. Figure 3(c) shows a matrix diagram of cutting operators. And Figure 3(d) shows the effect drawing of the cutting operator. At last, the wrist edge can be extracted.(2)Locating x coordinate: the outline contains many steps which are caused by two row adjacent pixels and have an influence on the curvature of it. Therefore, it is smoothed through a Butterworth lowpass digital filter B where the sample frequency is 70 Hz, the passband corner frequency is 2 Hz, the stopband corner frequency is 3 Hz, passband ripple is 3 dB, and stopband attenuation is 30 dB. Then, the curvature curve of the smoothed outline is calculated (Figure 2(d)). Taking the salient as a landmark, the distance between recess and GUAN is less than the width of a finger (approximately 100 pixels in our platform) according to TCM theory. When is the first local maximum of and the slope of is positive within the 100 pixels from the left of the salient, is considered as x coordinate of GUAN. Finally, the x coordinate of GUAN is marked with a vertical line in the thermal image (Figure 2(d)). The main process mentioned above is described in Algorithm 1.
2.3.2. Locating Y Coordinate of GUAN
A 10× 10 pixel region is established at each edge point in the edge image. The mean value and standard deviation of these regions are compared with the threshold values which are derived from a gray histogram of the radial artery region. As shown in Figure 2(e), when the mean value and standard deviation of the regions satisfy the threshold conditions, the pixels in the regions are selected and used to compose an image of the partial radial artery . Then, the mean values of y coordinate values in each column of (blue labels in Figure 2(f)) are obtained and fit an equation (red slash line in Figure 2(f)) for radial artery pattern using a leastsquares method. The y coordinate of GUAN is obtained by the equation and the x coordinate of GUAN. Finally, the position of GUAN is shown at the thermal image (Figure 2(e)) and the main process mentioned above is described in Algorithm 2.


3. Conversion between Pixel Distance and Actual Spatial Distance
The mapping of pixel distance to space distance plays an important role in the development of 3D mobile platform of desktop pulse diagnosis instrument. The camera is placed at an angle perpendicular to the ground. Due to the upward tilt of the elbow drag, the image obtained has trapezoid distortion, that is, the object above the image (near end) is wide; the object below the image (far end) is narrow. This trapezoid distortion will cause a large error in the subsequent calculation of threedimensional spatial coordinates, so it is unable to achieve the accurate positioning of the acquisition probe. According to the imaging principle of the camera's optical lens and other geometric relationships, the mapping relationship between the image coordinate system and the realworld coordinate system can be deduced. Only by knowing the camera's horizontal field angle, vertical field angle, the distance between the camera and the bracket, and the tilt angle of the elbow drag, the actual spatial coordinates of the guan position can be determined.
Assume that point P is the point location where the coordinate needs to be determined, and point O is the reference origin. Set the pixel position corresponding to P point as (x, y). Because of the upward inclination of the elbow drag, the calculation of the mapping value of the actual spatial coordinate of point P in the xaxis can be divided into two cases Figure 4. Case 1: as shown in Figure 4(a), point P is assumed to be on the left side of point G. Then, the step angle of point P can be expressed by the following formula: In formula (1), θ is half of the vertical field angle. Height is the height of the image. According to the sine theorem, the following formula can be obtained: In formula (2), is the tilt angle of the bracket, h is the distance between the camera and the bracket (perpendicular to the horizontal plane, and the intersection point with the horizontal plane is o point). Through the formula transformation, we can get the distance HP and the length of AP between the mapping point of point P in the xaxis and the distance O. According to the coordinate system direction, we can see the xaxis coordinates of point P as follows: Case 2: as shown in Figure 4(b), point P is assumed to be on the right side of point G. According to the sine theorem, the following formula can be obtained:
(a)
(b)
(c)
According to the coordinate system direction, we can see the xaxis coordinates of point P as follows:
As shown in Figure 4(c), point Q is the mapping point of point P on the yaxis.
Then, the step angle of point Q can be expressed by the following formula:
According to the sine theorem,
According to the coordinate system direction, we can see the yaxis coordinates of point P as follows:
In formula (10), sign is a sign function.
4. Experiments and Results
During the location process, each subject is required to hold the handle on the wrist support and removes the wrist after finishing image acquisition. An aluminum sheet with thermal insulation material (1 mm × 1 mm) is pasted on the wrist crease of a subject and is considered as a reference coordinate to eliminate position deviation (a thermal image with the aluminum sheet is shown in Figure 5). The process is repeated ten times and each subject is required to keep a unified posture. Eight volunteers (25.4 ± 2.1 years, 169.3 ± 6.4 cm, 59.4 ± 11.2 kg; 5 males and 3 females) from our lab participated in the experiment.
Table 1 shows x distance between GUAN and wrist crease (x coordinate of wrist crease − x coordinate of GUAN) and y distance of GUAN (y coordinate of GUAN). The standard deviations of x distance are all less than 1.7 mm. The standard deviations of y distance are almost less than 1 mm. The RSD of x distance is less than 9.0% and the RSD of y is less than 5.0%.

5. Discussion
5.1. The Effect of Hand Posture and Position
In this paper, an imaging platform based on thermal imager is designed for imaging radial processes and radial arteries. In this paper, the closed position imaging platform is used to maintain the wrist state of the subject and generate thermal imaging images. The imaging platform is mainly composed of the thermal imager, wrist bracket, a line laser, and a host computer. A line laser can ensure that the wrist is placed in roughly the same position and make sure that the wrist crease is at the edge of the image. Wrist branches can highlight the radial artery of the wrist. However, when placing the arm, the crease lines on the wrist should be included in the image. Meanwhile, the height of the arm should be slightly lower than that of the wrist. This can further ensure that the high posterior metacarpal bone presents a complete envelope.
5.2. The Effect of Location of GUAN
The x distance between wrist crease and GUAN is introduced for evaluating the accuracy of the x coordinate of GUAN. We set a reference point (the aluminum sheet) in order to eliminate position deviation and investigate deviations caused by other factors because position deviation has little influence on the guiding sensor.
Table 1 shows that x distance of GUAN is fairly reproducible (RSD< 9.0% and standard deviation < 1.7 mm). The deviation could be mainly due to wrist shape changes and could not be eliminated. However, it could be accepted because it could be covered by the sensitive area of a single sensor and it is also smaller than the spatial resolution of a sensor array. Y distance of GUAN is also fairly reproducible (RSD< 5.0% and standard deviation < 1.5 mm). The standard deviation of y distance is effective in locating y coordinate of GUAN for a single sensor or a sensor array. Furthermore, most of the deviation of y distance is smaller than the radial artery diameter. Therefore, the deviation of y distance provides valid values for guiding sensor. Conclusively, the proposed location algorithm is effective and necessary for locating GUAN.
5.3. The Limitation of This Study
In this paper, only eight young people are used to verify the accuracy of the system. With the growth of age and the loss of human collagen, blood vessels are more convex, so the GUAN position of the wrist of the elderly is more easily detected, so it is not considered. With high BMI subject, it is usually more difficult to discern the radial artery pulsation and styloid process. The detection accuracy of the system is affected by the accuracy of the infrared camera itself. However, the higher the accuracy, the higher the price of the infrared camera itself. For this defect, we will make the results in the manufacture of the trolleytype pulse diagnosis instrument. Meanwhile, in future work, we will improve the proposed location algorithm to decrease the standard deviations of x and y coordinates of GUAN.
6. Further Work
In future work, we will improve the proposed location algorithm and buy highprecision infrared thermal imaging equipment to decrease the standard deviations of x and y coordinates of GUAN. In addition, the wrist support, the thermal imager, and pressure sensors will be integrated together to achieve the realtime location of GUAN and the collection of the radial pulse waves.
7. Conclusions
In this paper, a novel method for locating GUAN based on thermal imaging and image processing is proposed. The infrared thermal imaging positioning imaging platform is used to ensure that the wrist of the subject enters the fixed imaging area in a fixed position during each measurement and transmits the thermal imaging images carrying the image information of radial processes and radial arteries to the upper computer. Arbitrary angle edge recognition algorithm is used to extract wrist contour and radial artery edge information. The xaxis coordinates of the radial artery were provided by the identification algorithm, and the yaxis coordinates of the radial artery were provided by the fitting algorithm. Finally, the x and y coordinates determine the GUAN position. The results of locating GUAN show that the proposed location algorithm not only has better repeatability but also provides a valid reference for collecting pulse waves of GUAN. Finally, the novel method could provide valid GUAN coordinates.
Data Availability
The raw/processed data required to reproduce these findings cannot be shared at this time as the data also form part of an ongoing study. But in the near future, with the permission of the institution, data sharing will be considered.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This study was supported by funds from the National Major Science and Technology Projects of China (no. 2018ZX01031201 and 2018ZX09201011), National Key Research and Development Program of China (2018YFC2001100), and the Chinese Academy of Sciences Key Project (no. KFZDSW218). The works presented in this paper are subject to pending China and international patents filed by Institute of Microelectronics of Chinese Academy of Sciences (IMECAS) in China (PCT/CN2017/092400, PCT/CN2017/092388, and PCT/CN2017/092399). The second author of the article, Su Liu, is currently working in NavInfo Co., Ltd. Although he has changed his job, the research results should belong to the original unit, so his article was signed as the original unit.
References
 L. Xu, M. Q.H. Meng, C. Shi, K. Wang, and N. Li, “Quantitative analyses of pulse images in Traditional Chinese Medicine,” Medical Acupuncture, vol. 20, no. 3, pp. 175–189, 2008. View at: Publisher Site  Google Scholar
 N. G. R. de Moura, I. Cordovil, and A. de Sá Ferreira, “Traditional Chinese medicine wrist pulsetaking is associated with pulse waveform analysis and hemodynamics in hypertension,” Journal of Integrative Medicine, vol. 14, no. 2, pp. 100–113, 2016. View at: Google Scholar
 C.H. Luo and C.Y. Chung, “Noninvasive holistic health measurements using pulse diagnosis: II. Exploring TCM clinical holistic diagnosis using an ingestion test,” European Journal of Integrative Medicine, vol. 8, no. 6, pp. 926–931, 2016. View at: Publisher Site  Google Scholar
 J. Bae, Y. J. Jeon, J. Y. Kim, and J. U. Kim, “New assessment model of pulse depth based on sensor displacement in pulse diagnostic devices,” Evidencebased Complementary and Alternative Medicine, vol. 2013, 2013. View at: Publisher Site  Google Scholar
 S. H. Jeon, K. K. Kim, I. S. Lee et al., “Pulse wave variation during the menstrual cycle in women with menstrual pain,” Biomed Research International, vol. 2016, Article ID 1083208, 8 pages, 2016. View at: Publisher Site  Google Scholar
 P. Wang, W. Zuo, and D. Zhang, “A compound pressure signal acquisition system for multichannel wrist pulse signal analysis,” IEEE Transactions on Instrumentation and Measurement, vol. 63, no. 6, pp. 1556–1565, 2014. View at: Google Scholar
 S. Gong, B. Xu, G. Sun et al., Accurate Cirrhosis Identification with WristPulse Data for Mobile Healthcare, ACM, New York, NY, USA, 2012.
 Y.F. Chung, C.S. Hu, C.C. Yeh, and C.H. Luo, “How to standardize the pulsetaking method of traditional Chinese medicine pulse diagnosis,” Computers in Biology and Medicine, vol. 43, no. 4, pp. 342–349, 2013. View at: Publisher Site  Google Scholar
 D. Wang, D. Zhang, and G. Lu, “A novel multichannel wrist pulse system with different sensor arrays,” IEEE Transactions on Instrumentation and Measurement, vol. 64, no. 7, pp. 2020–2034, 2015. View at: Publisher Site  Google Scholar
 C. Hu, Y. Chung, C. Yeh, and C. Luo, “Temporal and spatial properties of arterial pulsation measurement using pressure sensor array,” Evidencebased Complementary and Alternative Medicine, vol. 2012, Article ID 745127, 2012. View at: Publisher Site  Google Scholar
 Y.F. Chung, C.S. Hu, C.H. Luo et al., “Possibility of quantifying TCM fingerreading sensations: II. An example of health standardization,” European Journal of Integrative Medicine, vol. 4, no. 3, pp. e263–e270, 2012. View at: Publisher Site  Google Scholar
 C.H. Luo, Y.F. Chung, C.S. Hu et al., “Possibility of quantifying TCM fingerreading sensations: I. Bisensing pulse diagnosis instrument,” European Journal of Integrative Medicine, vol. 4, no. 3, pp. e255–e262, 2012. View at: Publisher Site  Google Scholar
 F. Clemente, P. Arpaia, and P. Cimmino, “A piezofilmbased measurement system for global haemodynamic assessment,” Physiological Measurement, vol. 31, no. 5, pp. 697–714, 2010. View at: Publisher Site  Google Scholar
 C.C. Tyan, S.H. Liu, J.Y. Chen, J.J. Chen, and W.M. Liang, “A novel noninvasive measurement technique for analyzing the pressure pulse waveform of the radial artery,” IEEE Transactions on Biomedical Engineering, vol. 55, no. 1, pp. 288–297, 2008. View at: Publisher Site  Google Scholar
 M.H. Jun, Y.M. Kim, J.H. Bae, C. Jung, J.H. Cho, and Y. Jeon, “Development of a tonometric sensor with a decoupled circular array for precisely measuring radial artery pulse,” Sensors, vol. 16, no. 6, p. 768, 2016. View at: Publisher Site  Google Scholar
 R. R. Anderson and J. A. Parrish, “The optics of human skin,” Journal of Investigative Dermatology, vol. 77, no. 1, pp. 13–19, 1981. View at: Publisher Site  Google Scholar
 K. P. Nielsen, L. Zhao, J. J. Stamnes, K. Stamnes, and J. Moan, “The optics of human skin: aspects important for human health,” Solar Radiation and Human Health, vol. 1, pp. 35–46, 2008. View at: Google Scholar
 I. Pavlidis, J. Dowdall, N. Sun, C. Puri, J. Fei, and M. Garbey, “Interacting with human physiology,” Computer Vision and Image Understanding, vol. 108, no. 12, pp. 150–170, 2007. View at: Publisher Site  Google Scholar
 J. T. Mehigan and R. A. McAlexander, “Snuffbox arteriovenous fistula for hemodialysis,” The American Journal of Surgery, vol. 143, no. 2, pp. 252253, 1982. View at: Publisher Site  Google Scholar
 R. A. Sanders, F. L. Keppel, and J. I. Waldrop, “External fixation of distal radial fractures: results and complications,” The Journal of Hand Surgery, vol. 16, no. 3, pp. 385–391, 1991. View at: Publisher Site  Google Scholar
 T. Nakamura, W. P. Cooney, W.H. Lui et al., “Radial styloidectomy: a biomechanical study on stability of the wrist joint,” The Journal of Hand Surgery, vol. 26, no. 1, pp. 85–93, 2001. View at: Publisher Site  Google Scholar
 S. Liu, X. Geng, Y. Zhang, S. zhang et al., “An edge detection method based on wavelet transform at arbitrary angles,” IEICE Transactions on Information and Systems, vol. E101–D, no. 9, pp. 2392–2400, 2018. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Xingguang Geng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.