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 x-axis coordinates of the radial artery were provided by the identification algorithm, and the y-axis 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.
National Major Science and Technology Projects of China2018ZX01031-2012018ZX09201011National Key Research and Development Program of China2018YFC2001100Chinese Academy of Sciences Key ProjectKFZD-SW-218Institute of MicroelectronicsChinese Academy of SciencesPCT/CN2017/092400PCT/CN2017/092388PCT/CN2017/0923991. 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 Methods2.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.
Photograph of the thermal imaging platform.
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:
A compact connection and a sparse connection are used to achieve the pixel combination principle of a corner.
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.
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.
Schematic illustration of locating GUAN. (a) Thermal image of the wrist, (b) edge image of the wrist, (c) the processed edge image of the wrist, (d) schematic illustration of locating x coordinate of GUAN, (e) the radial artery region segmented by the grayscale features of the edge point neighborhood, and (f) thermal image marked with GUAN.
Connecting operator and cutting operator. (a) Matrix diagram of the connecting operator. (b) Effect drawing of the connecting operator. (c) Matrix diagram of the cutting operator. (d) Effect drawing of the cutting operator.
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.
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 pwx is extracted by identifying the maximal connected component. Then, branches of pwx are cut by cutting operator along the long-axis direction of pwx. The cutting operator is chosen by a long-axis direction. When the long axis is along the x-axis, 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.
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 low-pass 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 ccx of the smoothed outline swx 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 ccx0 is the first local maximum of ccxand the slope of swx0is positive within the 100 pixels from the left of the salient, x0 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 Ri,j 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 RPi,j are selected and used to compose an image of the partial radial artery Ai,j. Then, the mean values of y coordinate values in each column of Ai,j (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 least-squares 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.
Algorithm 1: Algorithm of locating x coordinate of GUAN.
Extract the wx by identifying maximal connected component.
Smooth wx, swx=wx∗B.
Calculate ccx fromswx.ccx=sw''/1+sw'21.5.
Identify salient in swx and first local minimum of ccx near the salient.
Locate x coordinate of GUAN.
Algorithm 2: Algorithm of locating y coordinate of GUAN.
Establish Ri,j and calculate its mean value μR and standard deviation σR.
If the μR>t1&&σR<t2, where t1 is the mean gray value of border between the radial artery and other tissues and t2 is the standard deviation of radial artery region, t1=220,t2=25. RPi,j is used to compose the image of the partial radial artery Ai,j.
Obtain mean values of y values in each column of Ai,j and fit an equation.
Obtain y coordinate of GUAN.
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 three-dimensional 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 real-world 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 x-axis 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:
(1)Δθ=Xheight/2θ.
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:
(2)GPsinΔθ=hsin90−Δθ−α,(3)APsin90+α=GPsinΔθ.
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 x-axis and the distance O.
According to the coordinate system direction, we can see the x-axis coordinates of point P as follows:
(5)xP=−HP.
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) Geometric side view of thermal imaging platform (Case 1). (b) Geometric side view of thermal imaging platform (Case 2). (c) Top view of thermal imaging platform.
According to the coordinate system direction, we can see the x-axis coordinates of point P as follows:(7)xP=HP.
As shown in Figure 4(c), point Q is the mapping point of point P on the y-axis.
Then, the step angle of point Q can be expressed by the following formula:(8)Δβ=Ywidth/2β.
According to the sine theorem,(9)PQ=AP∗tanΔβ.
According to the coordinate system direction, we can see the y-axis coordinates of point P as follows:(10)yP=PQ∗signY.
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.
Experimental diagram of closing position.
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%.
Results of location of GUAN.
Subject
X distance (mm)
Y distance (mm)
Average
Standard deviation
RSD
Average
Standard deviation
RSD
1
19.55
1.41
0.072
33.14
0.25
0.008
2
14.44
0.41
0.028
31.93
0.19
0.006
3
15.81
1.09
0.069
32.62
0.50
0.015
4
19.53
1.65
0.084
36.90
0.46
0.012
5
26.03
1.56
0.060
30.52
1.47
0.048
6
16.91
1.35
0.080
26.54
0.63
0.024
7
20.31
1.27
0.063
35.46
1.07
0.030
8
24.58
1.45
0.059
32.42
0.45
0.014
5. Discussion5.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 trolley-type 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 high-precision 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 real-time 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 x-axis coordinates of the radial artery were provided by the identification algorithm, and the y-axis 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. 2018ZX01031-201 and 2018ZX09201011), National Key Research and Development Program of China (2018YFC2001100), and the Chinese Academy of Sciences Key Project (no. KFZD-SW-218). 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.
XuL.MengM. Q.-H.ShiC.WangK.LiN.Quantitative analyses of pulse images in Traditional Chinese Medicine200820317518910.1089/acu.2008.06322-s2.0-77953364571de MouraN. G. R.CordovilI.de Sá FerreiraA.Traditional Chinese medicine wrist pulse-taking is associated with pulse waveform analysis and hemodynamics in hypertension2016142100113LuoC.-H.ChungC.-Y.Non-invasive holistic health measurements using pulse diagnosis: II. Exploring TCM clinical holistic diagnosis using an ingestion test20168692693110.1016/j.eujim.2016.06.0162-s2.0-85007316403BaeJ.JeonY. J.KimJ. Y.KimJ. U.New assessment model of pulse depth based on sensor displacement in pulse diagnostic devices2013201310.1155/2013/9386412-s2.0-84885647630JeonS. H.KimK. K.LeeI. S.Pulse wave variation during the menstrual cycle in women with menstrual pain201620168108320810.1155/2016/10832082-s2.0-84984674356WangP.ZuoW.ZhangD.A compound pressure signal acquisition system for multichannel wrist pulse signal analysis201463615561565GongS.XuB.SunG.2012New York, NY, USAACMChungY.-F.HuC.-S.YehC.-C.LuoC.-H.How to standardize the pulse-taking method of traditional Chinese medicine pulse diagnosis201343434234910.1016/j.compbiomed.2012.12.0102-s2.0-84875377567WangD.ZhangD.LuG.A novel multichannel wrist pulse system with different sensor arrays20156472020203410.1109/tim.2015.24442582-s2.0-84959575127HuC.ChungY.YehC.LuoC.Temporal and spatial properties of arterial pulsation measurement using pressure sensor array2012201274512710.1155/2012/7451272-s2.0-80052730879ChungY.-F.HuC.-S.LuoC.-H.Possibility of quantifying TCM finger-reading sensations: II. An example of health standardization201243e263e27010.1016/j.eujim.2012.03.0042-s2.0-84867140039LuoC.-H.ChungY.-F.HuC.-S.Possibility of quantifying TCM finger-reading sensations: I. Bi-sensing pulse diagnosis instrument201243e255e26210.1016/j.eujim.2012.03.0032-s2.0-84867140070ClementeF.ArpaiaP.CimminoP.A piezo-film-based measurement system for global haemodynamic assessment201031569771410.1088/0967-3334/31/5/0072-s2.0-77951111441TyanC.-C.LiuS.-H.ChenJ.-Y.ChenJ.-J.LiangW.-M.A novel noninvasive measurement technique for analyzing the pressure pulse waveform of the radial artery200855128829710.1109/tbme.2007.9106812-s2.0-37349094157JunM.-H.KimY.-M.BaeJ.-H.JungC.ChoJ.-H.JeonY.Development of a tonometric sensor with a decoupled circular array for precisely measuring radial artery pulse201616676810.3390/s160607682-s2.0-84969922480AndersonR. R.ParrishJ. A.The optics of human skin1981771131910.1111/1523-1747.ep124791912-s2.0-0019813209NielsenK. P.ZhaoL.StamnesJ. J.StamnesK.MoanJ.The optics of human skin: aspects important for human health200813546PavlidisI.DowdallJ.SunN.PuriC.FeiJ.GarbeyM.Interacting with human physiology20071081-215017010.1016/j.cviu.2006.11.0182-s2.0-34548282008MehiganJ. T.McAlexanderR. A.Snuffbox arteriovenous fistula for hemodialysis1982143225225310.1016/0002-9610(82)90080-02-s2.0-0020080207SandersR. A.KeppelF. L.WaldropJ. I.External fixation of distal radial fractures: results and complications199116338539110.1016/0363-5023(91)90002-s2-s2.0-0025909814NakamuraT.CooneyW. P.LuiW.-H.Radial styloidectomy: a biomechanical study on stability of the wrist joint2001261859310.1053/jhsu.2001.209632-s2.0-0035129024LiuS.GengX.ZhangY.zhangS.An edge detection method based on wavelet transform at arbitrary angles2018E101–D92392240010.1587/transinf.2018edp71132-s2.0-85053851102