CN115868944A - Method for recognizing slow pulse characteristics by pulse envelope detection of cunguanchi three-part pressure sensing array - Google Patents

Method for recognizing slow pulse characteristics by pulse envelope detection of cunguanchi three-part pressure sensing array Download PDF

Info

Publication number
CN115868944A
CN115868944A CN202211165274.7A CN202211165274A CN115868944A CN 115868944 A CN115868944 A CN 115868944A CN 202211165274 A CN202211165274 A CN 202211165274A CN 115868944 A CN115868944 A CN 115868944A
Authority
CN
China
Prior art keywords
pulse
characteristic
target
cun
guan
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202211165274.7A
Other languages
Chinese (zh)
Other versions
CN115868944B (en
Inventor
孙珂
李昕欣
孙毅
郑熙坤
杨恒
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangdong Xinhuangpu Joint Innovation Institute Of Traditional Chinese Medicine
Original Assignee
Guangdong Xinhuangpu Joint Innovation Institute Of Traditional Chinese Medicine
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Guangdong Xinhuangpu Joint Innovation Institute Of Traditional Chinese Medicine filed Critical Guangdong Xinhuangpu Joint Innovation Institute Of Traditional Chinese Medicine
Priority to CN202211165274.7A priority Critical patent/CN115868944B/en
Publication of CN115868944A publication Critical patent/CN115868944A/en
Application granted granted Critical
Publication of CN115868944B publication Critical patent/CN115868944B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

The embodiment of the invention provides a method for recognizing a slow pulse characteristic by detecting pulse envelopes through a cunguan-chi three-part pressure sensing array, and relates to the technical field of sensing. Pulse wave information of cun position, guan position and chi position of a tested person under M kinds of external pressure is obtained firstly. Then based on the three pulse wave information, three corresponding pulse form change information are obtained. Then obtaining three pulse moistening coefficients corresponding to cun position, guan position and chi position according to the three pulse shape change information; the pulse moistening coefficient reflects the shape change trend of the cross section of the pulse of the testee along with the increase of the pressure value of the external pressure. If the pulse body moistening coefficients corresponding to the cun position, the guan position and the chi position are all larger than or equal to the preset threshold value, determining that the pulse condition characteristics of the tested person accord with the slow pulse characteristics. Therefore, the pulse wave information of cun position, guan position and chi position of the tested person is analyzed to identify the pulse condition characteristics of the slow pulse, and the accuracy of pulse condition identification is ensured.

Description

Method for recognizing slow pulse characteristics by pulse envelope detection of cunguanchi three-part pressure sensing array
Technical Field
The invention relates to the technical field of sensing, in particular to a method for recognizing a slow pulse characteristic by detecting pulse envelopes through cunguan three pressure sensing arrays.
Background
The pulse condition is the image and the representation of the beating rhythm of the radial artery, such as the speed, the size, the intensity, the floating and sinking, the deficiency and excess, etc. The 28 pulse systems of TCM include 28 types of pulse conditions. Each pulse has its own characteristics, while the slow pulse, one of the 28 pulse conditions, is characterized by its large and soft appearance and its soft size.
In the pulse diagnosis method in the prior art, the pulse condition is generally analyzed and judged by adopting pulse wave signals acquired by a sensor. This way, there is no problem in differentiating between several simple pulse conditions, such as slippery pulse, astringent pulse, slow pulse, and rapid pulse. But the complex pulse condition characteristics of the slow pulse, which also relate to the pulse condition changes at the cun, guan and chi positions, cannot be identified.
Disclosure of Invention
The invention aims to provide a method for recognizing a slow pulse characteristic by detecting a pulse envelope by a cunguan-chi three-part pressure sensing array so as to solve the problems in the prior art.
Embodiments of the invention may be implemented as follows:
in a first aspect, the invention provides a method for recognizing a slow pulse characteristic by detecting a pulse envelope with a cunguan-chi three-part pressure sensing array, which comprises the following steps:
acquiring pulse wave information of cun position, guan position and chi position of a tested person under M kinds of external pressure respectively;
obtaining pulse body form change information corresponding to the cun position, the guan position and the chi position respectively based on pulse wave information corresponding to the cun position, the guan position and the chi position respectively;
obtaining pulse body moistening coefficients corresponding to the cun position, the guan position and the ulnar position according to pulse body form change information corresponding to the cun position, the guan position and the ulnar position respectively;
wherein, the pressure values of the M kinds of the external pressure are gradually increased; the pulse body moistening coefficient reflects the form change trend of the cross section of the pulse body of the testee along with the increase of the pressure value;
and if the pulse body moistening coefficients corresponding to the cun position, the guan position and the chi position are all larger than or equal to a preset threshold value, determining that the pulse condition characteristics of the measured person accord with the slow pulse characteristics.
In a second aspect, the present invention provides a device for recognizing a slow pulse characteristic by detecting a pulse envelope with a cunguan-chi three-part pressure sensing array, comprising:
the acquisition module is used for acquiring pulse wave information of the cun position, the guan position and the chi position of the measured person under M kinds of external pressure respectively;
a processing module to:
obtaining pulse body form change information corresponding to the cun position, the guan position and the chi position respectively based on pulse wave information corresponding to the cun position, the guan position and the chi position respectively;
obtaining pulse body moistening coefficients corresponding to the cun position, the guan position and the ulnar position according to pulse body form change information corresponding to the cun position, the guan position and the ulnar position respectively; wherein the pressure value of the external pressure of the M kinds of the external pressures is gradually increased; the pulse body moistening coefficient reflects the form change trend of the cross section of the pulse body of the testee along with the increase of the pressure value;
and if the pulse body moistening coefficients corresponding to the cun position, the guan position and the chi position are all larger than or equal to a preset threshold value, determining that the pulse condition characteristics of the measured person accord with the slow pulse characteristics.
In a third aspect, the present invention provides an electronic device comprising: a memory and a processor, wherein the memory stores machine readable instructions executable by the processor, and when the electronic device runs, the processor executes the machine readable instructions to implement the method for detecting the pulse envelope and identifying the bradycardia characteristic by the three-dimensional pressure sensing array according to any one of the foregoing embodiments.
In a fourth aspect, the present invention provides a computer-readable storage medium storing a computer program, which is executed by a processor to implement the method for detecting pulse envelope and identifying bradycardia features of a pulse sensing array of the cun-guan-scale three-part pressure sensing according to any of the foregoing embodiments.
Compared with the prior art, the embodiment of the invention provides a method for recognizing the slow pulse characteristics by detecting the pulse envelope by the cun-guan-chi three-part pressure sensing array, which is characterized in that the cun position, the guan position and the chi position of a detected person are respectively pulse wave information under M kinds of external pressure. And then obtaining pulse body form change information corresponding to the cun position, the guan position and the chi position respectively based on the pulse wave information corresponding to the cun position, the guan position and the chi position respectively. Then obtaining the pulse body moistening coefficients corresponding to the cun position, the guan position and the ulna position according to the pulse body form change information corresponding to the cun position, the guan position and the ulna position respectively; the pulse soft-superficial coefficient reflects the trend of the shape change of the cross section of the pulse of the subject along with the increase of the pressure value of the external pressure. If the pulse body moistening coefficients corresponding to the cun position, the guan position and the chi position are all larger than or equal to the preset threshold value, determining that the pulse condition characteristics of the tested person accord with the slow pulse characteristics. Therefore, the pulse wave information of the cun position, the guan position and the ulnar position of the tested person is analyzed to identify the pulse condition characteristics of the slow pulse, and the accuracy of pulse condition identification is ensured.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the embodiments will be briefly described below, it should be understood that the following drawings only illustrate some embodiments of the present invention and therefore should not be considered as limiting the scope, and for those skilled in the art, other related drawings can be obtained according to the drawings without inventive efforts.
Fig. 1 is a schematic view of an application scenario of a sensor unit and a schematic view of a structure of a sensor array according to an embodiment of the present invention.
Fig. 2 is a first flowchart of a method for recognizing a slow pulse feature by detecting a pulse envelope with a cun-guan-chi pressure sensor array according to an embodiment of the present invention.
Fig. 3 is a schematic diagram illustrating analysis of a slow pulse feature recognition principle according to an embodiment of the present invention.
Fig. 4 is a second flowchart of a method for recognizing a slow pulse feature by detecting a pulse envelope with a cun-guan-chi three-part pressure sensing array according to an embodiment of the present invention.
Fig. 5 is a first schematic diagram of a lateral envelope curve according to an embodiment of the present invention.
FIG. 6 is a schematic diagram of processing a point to be smoothed using the Catmull-Rom Splines algorithm.
Fig. 7 is a second schematic diagram of a lateral envelope curve according to an embodiment of the present invention.
Fig. 8 is a third schematic diagram of a lateral envelope curve according to an embodiment of the present invention.
FIG. 9 shows an embodiment of the present invention providing a cun-position external pressure F 1 Corresponding 6 transverse envelope curves.
FIG. 10 shows an embodiment of the present invention providing a cun-position external pressure F 2 Corresponding 6 transverse envelope curves.
FIG. 11 shows an external pressure F for the cun position provided by an embodiment of the present invention 3 Corresponding 5 transverse envelope curves.
Fig. 12 is a schematic view of a linear fit straight line corresponding to the cunning position according to an embodiment of the present invention.
FIG. 13 shows an embodiment of the present invention providing a shut-off external pressure F 1 Corresponding 6 transverse envelope curves.
FIG. 14 shows an embodiment of the present invention providing a shut-off external pressure F 2 Corresponding 6 transverse envelope curves.
FIG. 15 shows an embodiment of the present invention providing a shut-off external pressure F 3 Corresponding 6 transverse envelope curves.
Fig. 16 is a schematic view of a linear fit straight line corresponding to the off-position according to an embodiment of the present invention.
FIG. 17 shows an embodiment of the present invention providing a gauge external pressure F 1 Corresponding 5 transverse envelope curves.
FIG. 18 shows the external pressure F at the scale position provided by the embodiment of the invention 2 Corresponding 6 transverse envelope curves.
FIG. 19 shows an embodiment of the present invention providing a gauge external pressure F 3 Corresponding 6 transverse envelope curves.
Fig. 20 is a schematic view of a linear fit straight line corresponding to the scale according to the embodiment of the invention.
Fig. 21 is a schematic structural diagram of a device for recognizing a slow pulse feature by detecting a pulse envelope with a cunguan-chi three-part pressure sensing array according to an embodiment of the present invention.
Fig. 22 is a schematic structural diagram of an electronic device according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention.
Thus, the following detailed description of the embodiments of the present invention, presented in the figures, is not intended to limit the scope of the invention, as claimed, but is merely representative of selected embodiments of the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
It should be noted that: like reference numbers and letters refer to like items in the following figures, and thus, once an item is defined in one figure, it need not be further defined and explained in subsequent figures. Furthermore, the appearances of the terms "first," "second," and the like, if any, are only used to distinguish one description from another and are not to be construed as indicating or implying relative importance. It should be noted that the features of the embodiments of the present invention may be combined with each other without conflict.
Cun-kou pulse method is to identify pulse condition by measuring the floating, sinking, slow, rapid, strong and weak pulse signals of radial artery at cunkou position of wrist with fingers. The currently accepted pulse conditions can be divided into 28 types.
In the prior art, an instrument is used for identifying pulse conditions, but because the size of the adopted pressure sensor or strain sensor is larger, in the existing pulse diagnosis method, the number of the adopted sensors is generally not more than 3, the sensor array is generally used for measuring pulse waves, the sensor array is pressed at the radial artery of the cunkou through external loading, a pulse wave time sequence curve changing along with time is obtained, and then the pulse wave time sequence curve is further analyzed to identify the pulse conditions. However, this method can only identify a few simple pulse conditions, and cannot accurately identify the slow pulse.
Based on the above technical problems, the inventors have made creative efforts to propose the following technical solutions to solve or improve the above problems. It should be noted that the above prior art solutions have shortcomings which are the results of practical and careful study of the inventor, therefore, the discovery process of the above problems and the solutions proposed by the embodiments of the present application in the following description should be the contribution of the inventor to the present application in the course of the invention creation process, and should not be understood as technical contents known by those skilled in the art.
For the complex pulse condition of the slow pulse, the ancient book describes the characteristics of the complex pulse condition, including: "it means gentle, even going, equal up and down, floating big and soft", "slow, not tight, and even going, i.e. the slow pulse has the feeling of even and soft under the fingers and longitudinal and slow pulse. In general, the core characteristic of slow pulse is "big and soft in shape, soft in size".
Therefore, the slow vessels have the characteristics of soft and soft, so that the cross-sectional shape of the radial artery of the wrist tends to be widened along with the increase of finger pressure. The inventor finds out through long-term observation and investigation that the pulse wave pressure signals acquired by a plurality of parallel sensors at the same time can form a cross-section transverse envelope curve. Therefore, pulse wave pressure signals of the cun-guan-chi parts of the wrists under different pressures are collected, transverse envelope curves of the cun-guan-chi parts under different pressures can be obtained, the curves are further analyzed, morphological change trends of the pulse body cross sections corresponding to the cun-guan-chi parts of the testee can be obtained, and therefore the pulse delaying characteristics can be identified by judging the soft-superficial degree of the pulse body.
Therefore, the method for recognizing the slow pulse features by detecting the pulse envelopes by the cun-guan-chi pressure sensing array provided by the embodiment of the invention can obtain the pulse body form change information of the cun position, the guan position and the chi position respectively by utilizing the pulse wave information of the cun position, the guan position and the chi position under M kinds of external pressures, and further recognize the slow pulse features. The following detailed description is made by way of examples, with reference to the accompanying drawings.
Referring to fig. 1, fig. 1 is a schematic view of an application scenario of a sensor unit according to an embodiment of the present invention. The radial artery at the cun-kou is divided into cun-guan-chi region, also called cun-position, guan-position and chi-position. Wherein, the sensor unit covers the radial artery at the cunkou of the wrist of the tested person.
Alternatively, the sensor unit may be composed of three identical sensor arrays 100, each sensor array 100 covers the radial artery vessel at the middle, and the sensor arrays 100 are arranged in a direction approximately perpendicular to the flow direction of the radial artery.
Referring to fig. 1, a sensor array 100 may be composed of a plurality of micro sensors 110 uniformly arranged, and the plurality of micro sensors 110 are packaged on a flexible substrate 120. The micro sensor 110 may be a micro pressure sensor. Where K may represent a center-to-center distance between two adjacent micro sensors 110, and K represents a center-to-center distance between two outermost micro sensors 110.
Due to the width limitation at the radial artery of the wrist of the human body, there is also a certain requirement for the size of the sensor array 100. In an alternative example, k may satisfy not more than 1.25mm and k may satisfy 6.25mm or more. The number of micro sensors 110 in each sensor array 100 is not limited to that shown in fig. 2, and may be not less than 5.
When data acquisition is performed on a measured person, the sensor array 100 in the sensor unit can acquire a pulse wave pressure signal of the measured person at a preset sampling frequency not lower than 16Hz, and the signal sampling time difference of any two micro sensors 110 can be controlled between 0 and 0.1 second.
Referring to fig. 1, the pulse wave pressure signals collected at the same time by the sensor array 100 corresponding to the cun position, the guan position, and the chi position, respectively, may be used to generate the transversal envelope curves S1, S2, and S3, respectively.
It should be noted that the structure of the sensor unit is not limited to the above example, and in another example, the sensor unit may only include one sensor array, in which three rows of sequentially arranged micro sensors are packaged on the flexible substrate.
The electronic device in the embodiment of the present invention may be, but is not limited to, a computer, a personal computer, a smart phone, a server, or a wearable device (in which the sensor unit is integrated), and the embodiment of the present invention is not limited in this respect.
Referring to fig. 2, fig. 2 is a schematic flow chart of a method for recognizing a slow pulse feature by detecting a pulse envelope with a cunguan-chi pressure sensing array according to an embodiment of the present invention. The method is applied to the electronic equipment and can comprise the following steps S100-S400:
s100, obtaining pulse wave information of the cun position, the guan position and the chi position of the measured person under M kinds of external pressure.
When the sensor unit is used for acquiring pulse wave pressure signals of a detected person, pressure values with different sizes are applied to the radial artery at the cunkou, and the detected pulse wave pressure signals are different in size. When the data acquisition is finished, the sampling data acquired by the sensor unit can be obtained, and the sampling data can comprise pulse wave information of cun position, guan position and chi position under M kinds of external pressure.
And S200, obtaining pulse body form change information corresponding to the cun position, the guan position and the chi position respectively based on the pulse wave information corresponding to the cun position, the guan position and the chi position respectively.
It can be understood that, for any one of the cun position, the guan position and the chi position, the pulse shape change information of the pulse taking part can reflect the deformation trend of the pulse taking part under different external pressures.
S300, obtaining the pulse moistening coefficients corresponding to the cun position, the guan position and the chi position according to the pulse form change information corresponding to the cun position, the guan position and the chi position.
Wherein, the pressure values of the M kinds of external pressure are gradually increased. The pulse soft-superficial coefficient can reflect the shape change trend of the cross section of the pulse of the testee along with the increase of the pressure value of the external pressure.
That is, for any one of cun, guan and chi pulse-taking parts, the pulse-taking part form-changing tendency with the increase of the pressure value of the external pressure can be reflected by the pulse-taking part form-changing coefficient, so as to reflect the pulse-taking part form-changing degree. The degree of moistening reflects the degree of deformation of the pulse body under external pressure. The higher the degree of moistening, the more obvious the deformation of the pulse body under pressure; the lower the degree of moistening, the lower the degree of deformation of the pulse body when it is stressed.
S400, if the pulse moistening coefficients corresponding to the cun position, the guan position and the chi position are all larger than or equal to a preset threshold value, determining that the pulse condition characteristics of the measured person accord with the slow pulse characteristics.
When the pulse-moistening coefficients corresponding to cun, guan and chi positions are all greater than or equal to the predetermined threshold, it is determined that the pulse-moistening degrees of cun, guan and chi positions are all high, and the core features of large and soft pulse-relaxing are satisfied.
The method for recognizing the slow pulse characteristic by measuring the pulse envelope by the cun-guan-chi three-part pressure sensing array provided by the embodiment of the invention obtains three corresponding pulse body form change information by utilizing the pulse wave information of the cun position, the guan position and the chi position of a measured person under M external pressures, and then obtains pulse body moistening coefficients corresponding to the cun position, the guan position and the chi position respectively based on the three pulse body form change information, wherein the pulse body moistening coefficients reflect the form change trend of the cross section of the measured person along with the increase of the pressure value of the external pressure. When the three pulse body moistening coefficients are all larger than or equal to the preset threshold value, determining that the pulse condition characteristics of the detected person accord with the slow pulse characteristics. Therefore, the pulse wave information of cun position, guan position and chi position is utilized to process and analyze, and the accuracy of pulse condition identification is ensured.
The principle of recognizing the characteristics of the slow pulse is first introduced here.
Referring to fig. 3, the external pressure 1, the external pressure 2 and the external pressure 3 are sequentially increased, and a transverse envelope curve of each external pressure can be formed by acquiring pulse wave pressure signals of the three cun-guan-chi parts as shown in fig. 3.
Observing the transverse envelope curve in fig. 3, it can be seen that there is a phenomenon that the curve shape is similar to being squeezed to two sides as the pressure value of the external pressure increases, i.e. as the pressure value of the external pressure increases, d1, d2, d3 marked in the figure become larger in sequence. This corresponds to the characteristics of the slow pulse, and therefore, whether the pulse condition characteristics correspond to the characteristics of the slow pulse can be determined based on analyzing the trend of the change of the form of the transverse envelope curve along with the increase of the pressure value of the external pressure.
In an alternative example, the sampled data is available using a sensor unit, which may include three sensor arrays. The sampling data comprises pulse wave information corresponding to cun position, guan position and chi position respectively, the pulse wave information corresponding to cun position, guan position and chi position is collected by three sensor arrays corresponding to cun position, guan position and chi position respectively, and each sensor array comprises a plurality of micro sensors.
Each type of pulse wave information may include all pulse wave pressure signals acquired by each miniature sensor over N pulse cycles at each external pressure. The N pulse periods may be continuous, partially continuous, or dispersed, and are not limited herein.
It should be noted that M and N are positive integers, and M and N may satisfy N ≧ 3 and M ≧ 3 in order to ensure the accuracy of pulse condition identification.
The following describes a process of obtaining pulse shape change information corresponding to the target pulse-taking part based on pulse wave information corresponding to the target pulse-taking part, taking the target pulse-taking part as an example.
Wherein, the target pulse-taking part can be any one of cun position, guan position and ulna position.
On the basis of fig. 3, please refer to fig. 4, the sub-steps of the step S200 may include S210 to S250.
S210, pulse wave information corresponding to the target pulse taking part is obtained.
Pulse wave information corresponding to the target pulse taking part can be extracted from the sampling data. The pulse wave information may include M kinds of pulse wave data corresponding to the M kinds of external pressures.
Each type of pulse wave data may include pulse wave pressure signals acquired by each micro-sensor over N pulse cycles at a corresponding external pressure.
In the processing, the same processing is performed for each type of pulse wave data corresponding to the external pressure to obtain the characteristic distance average value of the target pulse-taking portion corresponding to the external pressure, and steps S220 to S240 will be described below in terms of processing one type of pulse wave data corresponding to the external pressure.
S220, aiming at any external pressure, generating a transverse envelope curve corresponding to the target pulse taking position in each pulse period under the external pressure based on pulse wave information corresponding to the target pulse taking position, and obtaining N transverse envelope curves.
It is understood that, for any kind of pulse wave data corresponding to external pressure, a transverse envelope curve corresponding to the target pulse taking position in each pulse period under the external pressure can be generated based on the pulse wave data. For example, when N is 6, 6 transverse envelope curves corresponding to the target pulse-taking portion within 6 pulse periods under the external pressure are obtained.
The pulse beat frequency and the heart beat frequency of a human body are generally consistent, and are normally 60 to 100 times per minute. The pulse cycle is counted as the passing of one pulse cycle every time the pulse beats. The sensor array may perform multiple data acquisitions during each pulse cycle.
The preset sampling frequency of the sensor unit determines the number of sampling time points in each pulse period. In each pulse period, the number s of sampling time points depends on the size T of the pulse period and the size f of the preset sampling frequency: s = T × f. At each sampling time point, each micro sensor of the sensor unit will acquire a pulse wave pressure signal, which has a unit of kPa.
The process of generating the transversal envelope curve corresponding to the pulse period in step S220 is briefly described below by taking one pulse period as an example.
The substep of step S220 may include steps S221 to S223 for each pulse period under the external pressure using the pulse wave data corresponding to the external pressure.
S221, comparing the pulse wave pressure signals collected by the target micro sensor at each sampling time point in the pulse period, and taking the sampling time point corresponding to the maximum pulse wave pressure signal as the target sampling time point.
It is understood that the location of the target micro sensor may correspond to the center of the radial artery of the subject. In a pulse period, the sampling time point corresponding to the maximum pulse wave pressure signal acquired by the target micro sensor can be used as the target sampling time point.
S222, screening out pulse wave signals collected by each micro sensor at a target sampling time point in a pulse period from the pulse wave data corresponding to the external pressure.
And S223, obtaining a transverse envelope curve by using the pulse wave pressure signals acquired by each miniature sensor at the target sampling time point.
A coordinate system can be established by taking the sensor number of the micro-sensor as the horizontal axis and the signal intensity of the pulse wave pressure signal as the vertical axis. The pulse wave pressure signal acquired by each micro sensor at the target sampling time point can be marked as an independent point in a coordinate system, and the independent points corresponding to all the micro sensors are connected in sequence to obtain a transverse envelope curve.
In conjunction with fig. 2 and 5, assuming that the sensor array includes 18 micro sensors as an example, the center-to-center distance K between two adjacent micro sensors 110 is 0.65mm, and the center-to-center distance K between two outermost micro sensors 110 is 11.05mm. Referring to fig. 5, for a certain pulse period, 18 discrete independent points are obtained by using 18 pulse wave pressure signals collected at a target sampling time point in the pulse period, and the 18 independent points are smoothed, so that a transverse envelope curve can be obtained as shown in fig. 5. In fig. 5, the horizontal axis indicates the sensor number chn (0 to 17) of the micro sensor, and the vertical axis indicates the pressure value of the pulse wave pressure signal from which the dc component is removed, in kPa.
It should be noted that, in order to achieve a better resolution effect, the sampling time point corresponding to the maximum pulse wave signal in the pulse period is selected for the target sampling time period. The target sampling point time point may also be any sampling time point in the pulse cycle, and the selection of the specific target sampling point time point may be chosen based on the accuracy level of the data acquired by the sensor array.
In fig. 5, a smoothing algorithm is required to obtain the transverse envelope curve from 18 discrete points. In an alternative example, the smoothing algorithm may be a Catmull-Rom Splines algorithm, and the smoothing principle is briefly introduced as follows:
as in FIG. 6, assume that it is necessary to construct a point P to be smoothed 0 To P 1 Curve P (t), P between -1 and P2 Is a control point which determines the trend and curvature of the connecting line between two points to be smoothed.
P (t) is a cubic curve that can be expressed as to be solved:
P(t)=at 3 +bt 2 +ct+d
after solving, we get:
Figure BDA0003861084900000071
it is understood that the above examples are only examples, and other algorithms for smooth connection of discrete points may be used in practical applications, and are not limited herein.
And S230, calculating the characteristic distance corresponding to each transverse envelope curve to obtain N characteristic distances.
With reference to fig. 5, on the transverse envelope curve, the distance between two peaks can reflect the degree of deformation of the pulse-taking portion when pressed. The distance between two peaks can be laterally represented by a characteristic distance. That is, the characteristic distance may be indicative of a degree of deformation of the target site under the external pressure.
Optionally, the lateral envelope curve comprises peaks and valleys. The following describes calculating the characteristic distance corresponding to a transverse envelope curve by taking the transverse envelope curve as an example.
The substeps of the above step S230 may include S231 to S236:
and S231, determining a plurality of corresponding target sensor numbers according to any one transverse envelope curve.
In step S231, for a certain transverse envelope curve, the process of determining the corresponding plurality of target sensor numbers may include the following steps S2311 to S2313:
and S2311, pulse wave information acquired by the sensor array corresponding to the transverse envelope curve is acquired, and a time sequence oscillogram of each miniature sensor on a time sequence is generated.
For a transverse envelope curve corresponding to a certain pulse period at a certain external pressure at the target pulse taking part:
because the target pulse taking part is under the external pressure, the corresponding pulse wave data comprises: the target pulse taking part corresponds to pulse wave pressure signals collected by each micro sensor of the sensor array in N pulse periods. Therefore, the time-series waveform diagram corresponding to each micro-sensor can be generated by using the pulse wave data corresponding to the external pressure at the target pulse-taking part. The time-series waveform diagram can show the change of the pulse wave pressure signal at a certain point of the section of the target pulse-taking part along with the time.
And S2312, calculating curve similarity between each timing waveform diagram and a target timing waveform diagram corresponding to the target micro sensor.
As described above, in the sensor array corresponding to the target pulse taking position, the position of the target micro sensor may correspond to the center of the radial artery of the subject. The target time sequence oscillogram is a time sequence oscillogram corresponding to the target micro sensor.
It is understood that the way of calculating the curve similarity may be any existing method for calculating the curve similarity or a certain curve similarity calculation method modified by the existing method.
And S2313, all the sensor numbers of the micro sensors corresponding to the time sequence oscillograms with the curve similarity exceeding the set threshold are used as target sensor numbers.
It is understood that when the similarity of a curve exceeds a predetermined threshold, the micro-sensor corresponding to the similarity of the curve corresponding to the time-series waveform diagram can be regarded as an effective channel, and the sensor number of the micro-sensor can be regarded as the target sensor number.
The effective channel can be understood that a time sequence waveform diagram formed by the data collected by the micro sensor is consistent with a time sequence waveform diagram of a standard pulse condition, and the data collected by the micro sensor has referential property and analyzability. It should be noted that the size of the threshold may be determined according to practical applications, and is not limited herein.
And S232, determining a corresponding characteristic curve segment on the transverse envelope curve based on the plurality of target sensor numbers.
Normally, the plurality of target sensor numbers may be consecutive. Referring to fig. 7, assuming that a plurality of target sensor numbers 4 to 13 are obtained, the corresponding characteristic curve segment on the transverse envelope curve is the curve between points K1 to K2.
And S233, if only one convex peak exists on the characteristic curve segment, marking the characteristic distance corresponding to the transverse envelope curve as zero.
In a possible example, it may occur that only one upwardly convex peak is present on the characteristic curve segment. With reference to fig. 8, assuming that the number of the target sensor is 7 to 14, only one convex peak exists on the characteristic curve segment between the points K1 to K2, and at this time, the characteristic distance corresponding to the transverse envelope curve is directly marked as zero. The Span btw Gs in FIG. 8 represents the feature distance.
And S234, if at least two convex peaks exist on the characteristic curve segment, determining a first characteristic point and a second characteristic point corresponding to the characteristic curve segment.
It will be appreciated that the first characteristic point, the second characteristic point may each be located in a region between the characteristic curve segment and the horizontal axis of the transverse envelope curve.
One possible way to determine the first and second feature points is given below.
Optionally, the process of determining the first characteristic point and the second characteristic point corresponding to the characteristic curve segment in step S234 may include steps S2341 to S2343:
s2341, a horizontal straight line parallel to a horizontal axis of the transverse envelope curve is used for obtaining a surrounding area formed by the horizontal straight line and the characteristic curve section.
It should be noted that the ordinate of the horizontal straight line can be a set multiple of the maximum ordinate of the characteristic curve segment. The maximum ordinate of the characteristic curve segment is the ordinate corresponding to the highest point on the characteristic curve segment, and is also the height of the characteristic curve segment. The set multiple may be 1/3 to 1/2.
S2342, the enclosing region is divided into a first half section and a second half section with a longitudinal straight line passing through the intermediate division point and being parallel to the longitudinal axis of the lateral envelope curve as a boundary.
It will be appreciated that the intermediate division point may be located between two convex peaks on the characteristic curve segment. The highest upward convex peak on the characteristic curve segment may be referred to as the main peak, and the upward convex peak adjacent to the main peak on the left side (or adjacent to the right side) may be referred to as the neighboring peak. The intermediate division point may be located between the main peak and the neighboring peak on the characteristic curve segment. Specifically, the intermediate dividing point may be a trough between the main peak and the adjacent peak, or the abscissa of the intermediate dividing point may be the mean of the abscissas between the main peak and the adjacent peak (in this case, the intermediate dividing point is the same distance from the horizontal axes of the main peak and the adjacent peak).
In an alternative example, in connection with fig. 7, it is assumed that the ordinate of the horizontal straight line L1 is 1/3 times the maximum ordinate of the characteristic curve segment, and the intermediate division point through which the longitudinal straight line L2 passes is the trough K between the main peak and the adjacent peak. The hatched portion in fig. 7 is the enclosing region formed by the horizontal straight line L1 and the characteristic curve segment. The longitudinal straight line L2 divides the hatched portion into a first half section and a second half section.
S2343, determining a first feature point corresponding to the first half section and a second feature point corresponding to the second half section according to the preset height coefficient and the preset weight coefficient corresponding to each target sensor number.
Optionally, the coordinates of the feature points are (P) x ,P y ):
Figure BDA0003861084900000091
Figure BDA0003861084900000092
Wherein h is a preset height coefficient, and h belongs to [0,1].
When the characteristic point is a first characteristic point, vld _ chns is the number of target sensor numbers corresponding to the first half section, and x i Number of ith target sensor corresponding to first half section, y i The signal intensity of the pulse wave pressure signal corresponding to the ith target sensor number corresponding to the first half section is obtained; w is a i A preset weight coefficient corresponding to the number of the ith target sensor corresponding to the first half section;
when the characteristic point is a second characteristic point, vld _ chns is the number of target sensor numbers corresponding to the second half section, and x i Number of ith target sensor corresponding to second half section, y i The signal intensity of the pulse wave pressure signal corresponding to the ith target sensor number corresponding to the second half section is obtained; w is a i And the preset weight coefficient corresponds to the ith target sensor number corresponding to the second half section.
It should be noted that h and w may be adjusted for convenience of calculation i To obtain first and second feature points of different geometric definitions. The following lists 3 different first and second feature points:
(1) When h =0.5, w corresponding to each target sensor number i When the center of gravity of the first half section is equal to the center of gravity of the second half section, the first characteristic point and the second characteristic point respectively represent the center of gravity of the first half section and the center of gravity of the second half section, and the corresponding characteristic distance is the horizontal distance between the centers of gravity of the two half sections;
(2) When h =1, w corresponding to the maximum target sensor number and the minimum target sensor number i All are 1 and the rest are numbered with the corresponding w i When the first characteristic point and the second characteristic point are both 0, the first characteristic point and the second characteristic point respectively represent the maximum target sensor number and the minimum target sensor number which respectively correspond to the head end point and the tail end point on the characteristic curve segment (under special conditions, the longitudinal coordinates of the head end point and the tail end point are the same, and the corresponding characteristic distances are the sameI.e. the effective width of the surrounding area).
(3) When h =1, w corresponding to target sensor numbers of main peak and adjacent peak i All are 1 and the rest are numbered with the corresponding w i When the first characteristic point and the second characteristic point are both 0, the first characteristic point and the second characteristic point respectively represent peak points, namely a main peak and an adjacent peak, and the corresponding characteristic distance is the horizontal distance between the two peak points.
The example is only a special example, and the preset height coefficient and the preset weight coefficient corresponding to each target sensor number are based on the actual application situation, and are not limited to the above example.
And S235, calculating the horizontal distance between the first characteristic point and the second characteristic point to obtain the characteristic distance corresponding to the transverse envelope curve.
In connection with the example of fig. 7, when h =0.5, each target sensor is numbered with a corresponding w i When both are 1, the first feature point G1 and the second feature point G2 represent the gravity center points of the first half section and the second half section, respectively. The characteristic distance d is the horizontal distance between the two gravity points G1 and G2. The illustration in fig. 7 is only one possible example and is not intended as a limitation of the present embodiment.
It should be noted that, for a certain transverse envelope curve, if only one convex peak exists on the obtained characteristic curve segment, the characteristic distance corresponding to the transverse envelope curve is calculated, and steps S231 to S233 are executed; and if at least two convex peaks exist on the obtained characteristic curve segment, calculating the characteristic distance corresponding to the transverse envelope curve, and executing steps S231-S232 and S234-S235.
And S236, traversing each transverse envelope curve to obtain N characteristic distances.
And traversing each transverse envelope curve to obtain N characteristic distances corresponding to the N transverse envelope curves.
S240, calculating a characteristic distance mean value corresponding to the target pulse taking part under the external pressure based on the N characteristic distances.
It can be understood that N characteristic distances can be obtained by using pulse wave data corresponding to any external pressure, and the average value of the characteristic distances corresponding to the external pressure can be obtained by averaging the N characteristic distances. Likewise, the characteristic distance mean value can represent the deformation degree of the target vein-taking part under the external pressure.
And S250, traversing each external pressure to obtain a characteristic distance average value corresponding to the target pulse taking position under each external pressure.
The pulse shape change information of the target pulse taking part can comprise: and taking the characteristic distance mean value corresponding to the pulse part of the target under each external pressure. That is, the pulse shape change information of the target pulse taking part includes: m characteristic distance means corresponding to the M external pressures.
For each pulse taking part in cun position, guan position and ulna position, the above steps S210 to S250 and their respective substeps are performed to obtain the pulse body form change information of each pulse taking part: m characteristic distance means corresponding to the M external pressures.
Hereinafter, a method of calculating the soft pulse rate will be described by taking the target pulse-taking portion as an example.
The substeps of the above step S300 may include:
s310, performing linear fitting on M characteristic distance mean values corresponding to the target vein-taking part under M external pressures to obtain a linear fitting straight line corresponding to the target vein-taking part.
S320, taking the slope of the linear fitting straight line as a pulse body moistening coefficient corresponding to the target pulse part.
When the pressure value of the external pressure increases, the characteristic distance average value also increases, which indicates that the deformation degree of the target vein-taking part tends to increase, and when the pulse body moistening coefficient is greater than or equal to a preset threshold value, the pulse body moistening degree of the target vein-taking part is greater.
It should be noted that this embodiment is only one possible embodiment of the pulse-moistening coefficient.
In other possible embodiments, the pulse volume saturation coefficient of the target pulse-taking portion may be, by using the M characteristic distance mean values corresponding to M kinds of external pressures: and carrying out quadratic fitting on the M characteristic distance mean values under the M types of external pressure to obtain the slope of a quadratic fitting straight line, the correlation coefficient between the characteristic distance mean value and the external pressure and the like.
Alternatively, in step S400, the size of the preset threshold may be obtained by comparing data of the test group and the control group. A plurality of slow pulse testees are used as an experimental group, a plurality of flat pulse testees (normal pulse conditions) are used as a comparison group, pulse wave information of the three parts of cun, guan and chi of the experimental group is utilized to analyze to obtain pulse body form change information, pulse wave information of the three parts of cun, guan and chi of the comparison group is utilized to analyze to obtain pulse body form change information, and two groups of analysis results are compared to determine a preset threshold value.
It should be noted that, the execution sequence of each step in the foregoing method embodiments is not limited to that shown in the drawings, and the execution sequence of each step is subject to the practical application.
Assume that M is 3, N is 6, and the predetermined threshold is 0.04. The pulse wave information corresponding to cun position, guan position and chi position is used, and the example of determining the pulse condition of the person to be measured by using the corresponding transverse enveloping curves of cun position, guan position and chi position in each pulse period is given below.
M is 3, i.e. three external pressures F applied to the radial artery at the cuncial opening 1 、F 2 、F 3 Wherein the magnitude of the external pressure is F 1 <F 2 <F 3
In each of the transverse envelope curves in fig. 9 to 11, 13 to 15, and 17 to 19 of the drawings, the range of the corresponding target sensor number and the boundary between the first half section and the second half section are indicated by a dotted line, G1 and G2 respectively indicate the gravity center points of the first half section and the second half section, and Span btw Gs indicates the characteristic distance corresponding to the calculated transverse envelope curve.
First, the case where the cun position is under three kinds of external pressures is explained:
referring to FIG. 9, views (a) to (F) in FIG. 9 show the inch outside pressure F 1 6 transverse envelope curves corresponding to the next 6 pulse periods.
Referring to FIG. 10, the drawings (a) to (f) in FIG. 10 show the cun position outsidePressure F 2 6 transverse envelope curves corresponding to the next 6 pulse periods.
Referring to FIG. 11, FIGS. 11 (a) - (e) show the inch outside pressure F 3 5 transverse envelope curves corresponding to the next 5 pulse periods.
Statistical analysis was performed on the characteristic distances corresponding to each of the lateral envelope curves shown in fig. 9 to 11, as shown in table 1 below:
TABLE 1
Figure BDA0003861084900000111
Three kinds of external pressure F 1 、F 2 、F 3 The corresponding feature distance mean values are 2.15, 2.41 and 3.39 respectively, and a linear fitting straight line as shown in fig. 12 is obtained after linear fitting. The linear fitting straight line embodies the dynamic change trend that two half sections on the transverse envelope curve corresponding to cun position are gradually extruded to two sides along with the increase of the pressure value of external pressure, and conforms to the characteristic of pulse moistening.
In FIG. 12, the slope of the linear fit line is 0.62 and the goodness of fit R 2 =0.9. At this time, the slope was 0.62>0.04, the degree of pulse moistening in cun position was determined.
The goodness of fit refers to the degree of fit of a regression line (i.e., a linear fit line) to an observed value (i.e., a feature distance mean). Goodness of fit R 2 The closer to 1, the better the fitting degree of the regression straight line to the observed value is; otherwise, R 2 The smaller the regression line, the worse the degree of fitting to the observed values.
Next, the case where the off position is under three kinds of external pressures will be explained:
referring to fig. 13, in fig. 13, the diagrams (a) to (F) show the off-position external pressure F 1 6 transverse envelope curves corresponding to the next 6 pulse periods.
Referring to fig. 14, in fig. 14, the diagrams (a) to (F) show the off-position external pressure F 2 6 transverse envelope curves corresponding to the next 6 pulse periods.
Referring to fig. 15, in fig. 15, the diagrams (a) to (F) show the off-position external pressure F 3 6 transverse envelope curves corresponding to the next 6 pulse periods.
The characteristic distances corresponding to each of the lateral envelope curves shown in fig. 13 to 15 were statistically analyzed as shown in table 2 below:
TABLE 2
Figure BDA0003861084900000121
Three kinds of external pressure F 1 、F 2 、F 3 The corresponding feature distance mean values are 0.00, 1.87, and 2.38, respectively, and a linear fit straight line as shown in fig. 16 is obtained after linear fitting. This linear fitting straight line has embodied two half sections on the corresponding horizontal envelope curve of closed position and has been crowded to the dynamic variation trend of both sides along with the pressure value of external pressure increases gradually, accords with the characteristic that the pulse body is soft-shelled.
In FIG. 16, the slope of the linear fit line is 1.19 and the goodness of fit R 2 =0.9. At this time, the slope was 1.79>0.04, the degree of soft-superficial pulse at the point of closure is determined to be strong.
Finally, the description is given for the case of three external pressures for the scale:
referring to fig. 17, views (a) to (e) in fig. 17 show the scale external pressure F 1 The next 5 pulse periods correspond to 5 transverse envelope curves.
Referring to fig. 18, views (a) to (F) in fig. 18 show the scale external pressure F 2 6 transverse envelope curves corresponding to the next 6 pulse periods.
Referring to FIG. 19, the drawings (a) to (F) in FIG. 19 show the scale outside pressure F 3 6 transverse envelope curves corresponding to the next 6 pulse periods.
The feature distances corresponding to each of the lateral envelope curves shown in fig. 17 to 19 are counted, and the following table 3 is obtained:
TABLE 3
Figure BDA0003861084900000122
Three kinds of external pressure F 1 、F 2 、F 3 The corresponding feature distance mean values are 2.48, 2.90, and 3.15, respectively, and a linear fit straight line as shown in fig. 20 is obtained after linear fitting. The linear fitting straight line embodies the dynamic change trend that two half sections on the transverse envelope curve corresponding to the position of the ruler are gradually extruded to two sides along with the increase of the pressure value of the external pressure, and conforms to the characteristic of pulse moistening.
In FIG. 20, the slope of the linear fit straight line is 0.33, and the goodness of fit R 2 =0.98. At this time, the slope was 0.33>0.04, the degree of the pulse in the ulnar region was determined to be large.
The pulse condition characteristics of the tested person are judged to accord with the slow pulse characteristics by integrating the pulse softness conclusions of the cun position, the guan position and the ulna position.
It should be noted that the above example is only an exemplary illustration, and the specific way of determining the slow pulse is not limited thereto.
Compared with the prior art, the embodiment of the invention has the following beneficial effects:
the limitation that the pulse width information cannot be sensitively sensed by using a single-point sensor or a few isolated point sensors in the prior art is overcome;
the pulse wave information acquired by the sensor unit can be used for drawing a transverse envelope curve, the deformation degree of the pulse taking part under external pressure is reflected by the characteristic distance corresponding to the transverse envelope curve, and the limitation that the data acquired by the sensor can only be used for drawing a longitudinal time sequence oscillogram in the prior art is overcome;
analyzing the waveform changes of the three transverse envelope curves of cun position, guan position and chi position under different external pressures to indirectly reflect the change trend of the deformation degree of the pulse taking part, thereby quantifying the moistening degree of the pulse body, extracting more comprehensive and complete pulse condition characteristics and ensuring the accuracy of pulse condition characteristic identification.
In order to execute the corresponding steps in the above method embodiments and various possible embodiments, an implementation of the apparatus for recognizing a slow pulse feature by detecting a pulse envelope with a cun-guan-chi three-part pressure sensing array is given below.
Referring to fig. 21, fig. 21 is a schematic structural diagram of a device for recognizing a slow pulse characteristic by detecting a pulse envelope with a cunguan scale and a pressure sensing array according to an embodiment of the present invention. The apparatus 200 comprises: an acquisition module 210 and a processing module 220.
The obtaining module 210 is configured to obtain pulse wave information of the cun position, the guan position, and the chi position of the subject under the M kinds of external pressures.
A processing module 220 for:
obtaining pulse body form change information corresponding to cun position, guan position and chi position respectively based on pulse wave information corresponding to cun position, guan position and chi position respectively;
obtaining the pulse body moistening coefficients corresponding to the cun position, the guan position and the ulna position according to the pulse body form change information corresponding to the cun position, the guan position and the ulna position respectively; wherein, the pressure value of the external pressure in the M external pressures is gradually increased; the pulse soft-boiling coefficient reflects the form change trend of the cross section of the pulse of the testee along with the increase of the pressure value;
if the pulse body moistening coefficients corresponding to the cun position, the guan position and the chi position are all larger than or equal to the preset threshold value, determining that the pulse condition characteristics of the tested person accord with the slow pulse characteristics.
In an optional embodiment, the pulse wave information corresponding to the cun position, the off position and the chi position is respectively acquired by three sensor arrays corresponding to the cun position, the off position and the chi position, each sensor array comprises a plurality of micro sensors, and the pulse wave information comprises pulse wave pressure signals acquired by each micro sensor in N pulse periods. The processing module 220 may be specifically configured to:
acquiring pulse wave information corresponding to a target pulse taking part, wherein the target pulse taking part is any one of cun position, guan position and chi position; aiming at any external pressure, generating a transverse envelope curve corresponding to the target pulse taking position in each pulse period under the external pressure based on pulse wave information corresponding to the target pulse taking position to obtain N transverse envelope curves; calculating the characteristic distance corresponding to each transverse envelope curve to obtain N characteristic distances; the characteristic distance represents the deformation degree of the target pulse taking part under external pressure; calculating a characteristic distance mean value corresponding to the target pulse taking part under the external pressure based on the N characteristic distances; traversing each external pressure to obtain a characteristic distance average value corresponding to the target pulse taking part under each external pressure; wherein, the pulse body form change information of the target pulse taking part comprises: and taking the characteristic distance mean value corresponding to the pulse part of the target under each external pressure.
In an alternative embodiment, the horizontal axis of the transverse envelope curve corresponds to the number of sensors in the sensor array, wherein the sensors are sequentially arranged, and the transverse envelope curve comprises peaks and valleys. The processing module 220 may be specifically configured to:
determining a plurality of corresponding target sensor numbers according to any one transverse envelope curve; determining a corresponding characteristic curve segment on the transverse envelope curve based on the plurality of target sensor numbers; if only one convex peak exists on the characteristic curve segment, marking the characteristic distance corresponding to the transverse envelope curve as zero; if at least two convex peaks exist on the characteristic curve segment, determining a first characteristic point and a second characteristic point corresponding to the characteristic curve segment; the first characteristic point and the second characteristic point are both positioned in the area between the characteristic curve segment and the transverse axis of the transverse envelope curve; calculating the horizontal distance between the first characteristic point and the second characteristic point to obtain the characteristic distance corresponding to the transverse envelope curve; and traversing each transverse envelope curve to obtain N characteristic distances.
In an alternative embodiment, the processing module 220 may specifically be configured to: acquiring pulse wave information acquired by a sensor array corresponding to the transverse envelope curve, and generating a time sequence oscillogram of each miniature sensor on a time sequence; calculating the curve similarity between each time sequence oscillogram and a target time sequence oscillogram corresponding to the target micro sensor; the position of the target micro sensor corresponds to the center of the radial artery of the measured person; and taking the sensor numbers of the micro sensors corresponding to the time sequence oscillograms with the curve similarity exceeding the set threshold value as the target sensor numbers.
In an alternative embodiment, the vertical axis of the transverse envelope curve represents the signal strength of the pulse wave pressure signal. The processing module 220 may be specifically configured to: obtaining an enclosing area formed by the horizontal straight line and the characteristic curve section by using the horizontal straight line parallel to the horizontal axis of the transverse enveloping curve; the vertical coordinate of the horizontal straight line is a set multiple of the maximum vertical coordinate of the characteristic curve segment; dividing the surrounding area into a first half section and a second half section by taking a longitudinal straight line which passes through the middle division point and is parallel to the longitudinal axis of the transverse envelope curve as a boundary; the middle division point is positioned between the two convex wave crests on the characteristic curve segment; and determining a first characteristic point corresponding to the first half section and a second characteristic point corresponding to the second half section according to the preset height coefficient and the preset weight coefficient corresponding to each target sensor number.
In an optional embodiment, the processing module 220 may be further specifically configured to: performing linear fitting on M characteristic distance mean values corresponding to the target pulse taking part under M external pressures to obtain a linear fitting straight line corresponding to the target pulse taking part; the slope of the linear fitting line is used as the pulse body moistening coefficient corresponding to the target pulse part.
It is clear to those skilled in the art that, for convenience and brevity of description, the specific working process of the apparatus 200 described above may refer to the corresponding process in the foregoing method embodiment, and is not described herein again.
Referring to fig. 22, fig. 22 is a schematic structural diagram of an electronic device according to an embodiment of the present invention. The electronic device 300 includes a processor 310, a memory 320, and a bus 330, the processor 310 being coupled to the memory 320 via the bus 330.
The memory 320 can be used for storing a software program, such as a device for recognizing the slow pulse characteristic by detecting the pulse envelope of the cunguan-chi three-part pressure sensing array shown in fig. 21. The Memory 320 may be, but is not limited to, a Random Access Memory (RAM), a Read Only Memory (ROM), a Flash Memory (Flash), a Programmable Read-Only Memory (PROM), an Erasable Read-Only Memory (EPROM), an electrically Erasable Read-Only Memory (EEPROM), and the like.
The processor 310 may be an integrated circuit chip having signal processing capabilities. The Processor 310 may be a general-purpose Processor including a Central Processing Unit (CPU), a Network Processor (NP), and the like; but also Digital Signal Processors (DSPs), application Specific Integrated Circuits (ASICs), field Programmable Gate Arrays (FPGAs) or other Programmable logic devices, discrete Gate or transistor logic devices, discrete hardware components. Memory 320 stores machine-readable instructions executable by processor 310. The processor 310, when executing the machine-readable instructions, implements the method for recognizing the characteristics of the slow pulse by detecting the pulse envelope of the pulse sensing array of the cunguanchi type disclosed in the above embodiments.
It will be appreciated that the configuration shown in fig. 22 is merely illustrative and that electronic device 300 may include more or fewer components than shown in fig. 22 or have a different configuration than shown in fig. 22. The components shown in fig. 22 may be implemented in hardware, software, or a combination thereof.
The embodiment of the invention also provides a computer readable storage medium, wherein a computer program is stored on the computer readable storage medium, and when the computer program is executed by a processor, the method for detecting the pulse envelope and identifying the bradycardia characteristic by the inch-gate-size three-part pressure sensing array disclosed by the embodiment is realized. The computer readable storage medium may be, but is not limited to: u disk, removable hard disk, ROM, RAM, PROM, EPROM, EEPROM, FLASH disk or optical disk, etc.
To sum up, the embodiment of the invention provides a method for recognizing a slow pulse characteristic by detecting pulse envelopes through a cun-guan-chi three-part pressure sensing array. And then obtaining pulse body form change information corresponding to the cun position, the guan position and the chi position respectively based on the pulse wave information corresponding to the cun position, the guan position and the chi position respectively. Then obtaining the pulse body moistening coefficients corresponding to the cun position, the guan position and the ulna position according to the pulse body form change information corresponding to the cun position, the guan position and the ulna position respectively; the pulse soft-superficial coefficient reflects the trend of the shape change of the cross section of the pulse of the subject along with the increase of the pressure value of the external pressure. If the pulse body moistening coefficients corresponding to the cun position, the guan position and the chi position are all larger than or equal to the preset threshold value, determining that the pulse condition characteristics of the tested person accord with the slow pulse characteristics. Therefore, the pulse wave information of the cun position, the guan position and the ulnar position of the tested person is analyzed to identify the pulse condition characteristics of the slow pulse, and the accuracy of pulse condition identification is ensured.
The above description is only for the specific embodiments of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims (7)

1. A method for recognizing a slow pulse characteristic by detecting pulse envelopes through a cunguan three-part pressure sensing array is characterized by comprising the following steps:
acquiring pulse wave information of cun position, guan position and chi position of a tested person under M kinds of external pressure respectively;
obtaining pulse body form change information corresponding to the cun position, the guan position and the chi position respectively based on pulse wave information corresponding to the cun position, the guan position and the chi position respectively;
obtaining the pulse body moistening coefficients corresponding to the cun position, the guan position and the ulna position according to the pulse body form change information corresponding to the cun position, the guan position and the ulna position respectively;
wherein the pressure values of the M kinds of the external pressures are gradually increased; the pulse body moistening coefficient reflects the form change trend of the cross section of the pulse body of the testee along with the increase of the pressure value;
and if the pulse body moistening coefficients corresponding to the cun position, the guan position and the chi position are all larger than or equal to a preset threshold value, determining that the pulse condition characteristics of the measured person accord with the slow pulse characteristics.
2. The method according to claim 1, wherein the pulse wave information corresponding to the cun position, the off position, and the chi position is respectively collected by three sensor arrays corresponding to the cun position, the off position, and the chi position, each sensor array including a plurality of micro sensors, the pulse wave information including pulse wave pressure signals collected by each micro sensor in N pulse periods;
the step of obtaining the pulse body form change information corresponding to the cun position, the guan position and the chi position respectively based on the pulse wave information corresponding to the cun position, the guan position and the chi position respectively comprises the following steps:
acquiring the pulse wave information corresponding to the target pulse taking part; the target pulse taking part is any one of the cun position, the guan position and the chi position;
aiming at any external pressure, generating a transverse envelope curve corresponding to the target pulse taking position in each pulse period under the external pressure based on pulse wave information corresponding to the target pulse taking position to obtain N transverse envelope curves;
calculating a characteristic distance corresponding to each transverse envelope curve to obtain N characteristic distances; the characteristic distance represents the deformation degree of the target pulse taking part under the external pressure;
calculating a characteristic distance mean value corresponding to the target pulse taking position under the external pressure based on the N characteristic distances;
traversing each external pressure to obtain a characteristic distance average value corresponding to the target pulse taking position under each external pressure;
wherein, the pulse body form change information of the target pulse taking part comprises: and under each external pressure, the characteristic distance mean value corresponding to the pulse taking part of the target is obtained.
3. The method of claim 2, wherein the horizontal axis of the lateral envelope curve corresponds to the number of sensors in the array of sensors in which the microsensors are sequentially arranged;
the step of calculating the characteristic distance corresponding to each of the transverse envelope curves to obtain N characteristic distances includes:
aiming at any one transverse envelope curve, determining a plurality of corresponding target sensor numbers;
determining a corresponding characteristic curve segment on the transverse envelope curve based on a plurality of the target sensor numbers;
if only one convex peak exists on the characteristic curve segment, marking the characteristic distance corresponding to the transverse envelope curve as zero;
if at least two convex peaks exist on the characteristic curve segment, determining a first characteristic point and a second characteristic point corresponding to the characteristic curve segment; the first characteristic point and the second characteristic point are both located in a region between the characteristic curve segment and a transverse axis of the transverse envelope curve;
calculating the horizontal distance between the first characteristic point and the second characteristic point to obtain the characteristic distance corresponding to the transverse envelope curve;
and traversing each transverse envelope curve to obtain N characteristic distances.
4. The method of claim 3, wherein the step of determining a corresponding plurality of target sensor numbers comprises:
acquiring pulse wave information acquired by a sensor array corresponding to the transverse envelope curve, and generating a time sequence oscillogram of each miniature sensor on a time sequence;
calculating the curve similarity between each time sequence oscillogram and a target time sequence oscillogram corresponding to the target micro sensor; the position of the target micro sensor is close to the center of the radial artery of the testee;
and taking the sensor numbers of the micro sensors corresponding to the time sequence oscillogram with the curve similarity exceeding a set threshold value as target sensor numbers.
5. The method according to claim 3, wherein the vertical axis of the transverse envelope curve represents the signal strength of the pulse wave pressure signal; the step of determining the first characteristic point and the second characteristic point corresponding to the characteristic curve segment includes:
obtaining an enclosed area formed by the horizontal straight line and the characteristic curve section by using the horizontal straight line parallel to the horizontal axis of the transverse envelope curve; the vertical coordinate of the horizontal straight line is a set multiple of the maximum vertical coordinate of the characteristic curve segment;
dividing the enclosed region into a first half-section and a second half-section by taking a longitudinal straight line which passes through a middle division point and is parallel to a longitudinal axis of the transverse envelope curve as a boundary; the middle division point is positioned between the two convex peaks on the characteristic curve segment;
and determining a first characteristic point corresponding to the first half section and a second characteristic point corresponding to the second half section according to a preset height coefficient and a preset weight coefficient corresponding to each target sensor number.
6. The method of claim 5, wherein the coordinates of the feature points are (P) x ,P y ),
Figure FDA0003861084890000031
Figure FDA0003861084890000032
Wherein h is the preset height coefficient, and h belongs to [0,1]](ii) a When the characteristic point is the first characteristic point, vld _ chns is the number of the target sensor numbers corresponding to the first half section, and x is i The number y of the ith target sensor corresponding to the first half section i Compiling the ith target sensor corresponding to the first half sectionThe signal intensity of the pulse wave pressure signal corresponding to the signal; w is a i A preset weight coefficient corresponding to the ith target sensor number corresponding to the first half section;
when the characteristic point is the second characteristic point, vld _ chns is the number of the target sensor numbers corresponding to the second half section, and x is i The number y of the ith target sensor corresponding to the second half section i The signal intensity of the pulse wave pressure signal corresponding to the ith target sensor number corresponding to the second half section is obtained; w is a i And a preset weight coefficient corresponding to the ith target sensor number corresponding to the second half section.
7. The method according to claim 2, wherein the step of obtaining the pulse volume moistening coefficients corresponding to the cun position, the guan position and the ulna position according to the pulse volume morphological change information corresponding to the cun position, the guan position and the ulna position comprises:
performing linear fitting on M characteristic distance mean values corresponding to the target vein-taking part under M external pressures to obtain a linear fitting straight line corresponding to the target vein-taking part;
and taking the slope of the linear fitting straight line as the pulse body moistening coefficient corresponding to the target pulse taking part.
CN202211165274.7A 2022-09-23 2022-09-23 Method for measuring pulse envelope and identifying pulse buffering characteristics by three pressure sensing arrays of cunguan ruler Active CN115868944B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211165274.7A CN115868944B (en) 2022-09-23 2022-09-23 Method for measuring pulse envelope and identifying pulse buffering characteristics by three pressure sensing arrays of cunguan ruler

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211165274.7A CN115868944B (en) 2022-09-23 2022-09-23 Method for measuring pulse envelope and identifying pulse buffering characteristics by three pressure sensing arrays of cunguan ruler

Publications (2)

Publication Number Publication Date
CN115868944A true CN115868944A (en) 2023-03-31
CN115868944B CN115868944B (en) 2023-09-08

Family

ID=85769990

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211165274.7A Active CN115868944B (en) 2022-09-23 2022-09-23 Method for measuring pulse envelope and identifying pulse buffering characteristics by three pressure sensing arrays of cunguan ruler

Country Status (1)

Country Link
CN (1) CN115868944B (en)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101049247A (en) * 2006-04-06 2007-10-10 香港理工大学 Three portions and nine pulse-takings pulse condition detector of pulse condition sensor of herbalist doctor, and pulse condition detection method
KR20080084007A (en) * 2007-03-14 2008-09-19 한국 한의학 연구원 Diagnosis system of deficient and forceful pulse
CN102068241A (en) * 2005-04-24 2011-05-25 西安信唯信息科技有限公司 Pulse condition detection method for pulse condition restoring method
CN107736880A (en) * 2017-10-24 2018-02-27 新绎健康科技有限公司 A kind of pulse analysis method and system
CN110113990A (en) * 2016-12-28 2019-08-09 欧姆龙株式会社 Pulse wave measuring apparatus, pulse wave measurement method and blood pressure measuring device
WO2020132899A1 (en) * 2018-12-25 2020-07-02 深圳市大富网络技术有限公司 Blood pressure and pulse detection system, and detection method and device therefor
CN113729648A (en) * 2021-09-08 2021-12-03 华东理工大学 Wearable pulse-taking bracelet system based on multiple pulse sensors
CN114305345A (en) * 2022-02-16 2022-04-12 平安科技(深圳)有限公司 Pulse condition identification method, system, device and storage medium
CN114431840A (en) * 2022-04-08 2022-05-06 北京大学深圳研究生院 Pulse acquisition device, pulse acquisition method and system

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102068241A (en) * 2005-04-24 2011-05-25 西安信唯信息科技有限公司 Pulse condition detection method for pulse condition restoring method
CN101049247A (en) * 2006-04-06 2007-10-10 香港理工大学 Three portions and nine pulse-takings pulse condition detector of pulse condition sensor of herbalist doctor, and pulse condition detection method
KR20080084007A (en) * 2007-03-14 2008-09-19 한국 한의학 연구원 Diagnosis system of deficient and forceful pulse
CN101668476A (en) * 2007-03-14 2010-03-10 韩国韩医学研究院 Diagnosis system of deficient and forceful pulse
CN110113990A (en) * 2016-12-28 2019-08-09 欧姆龙株式会社 Pulse wave measuring apparatus, pulse wave measurement method and blood pressure measuring device
CN107736880A (en) * 2017-10-24 2018-02-27 新绎健康科技有限公司 A kind of pulse analysis method and system
WO2020132899A1 (en) * 2018-12-25 2020-07-02 深圳市大富网络技术有限公司 Blood pressure and pulse detection system, and detection method and device therefor
CN113729648A (en) * 2021-09-08 2021-12-03 华东理工大学 Wearable pulse-taking bracelet system based on multiple pulse sensors
CN114305345A (en) * 2022-02-16 2022-04-12 平安科技(深圳)有限公司 Pulse condition identification method, system, device and storage medium
CN114431840A (en) * 2022-04-08 2022-05-06 北京大学深圳研究生院 Pulse acquisition device, pulse acquisition method and system

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
于婷婷: "用于脉搏波检测的背接触式触觉传感阵列单元", 《传感技术学报》, vol. 33, no. 4, pages 492 - 499 *

Also Published As

Publication number Publication date
CN115868944B (en) 2023-09-08

Similar Documents

Publication Publication Date Title
KR102579840B1 (en) Obtaining blood pressure signals using a pressure sensor array
US20200170577A1 (en) System and Methods for Adaptive Noise Quantification in Dynamic Biosignal Analysis
Jiang et al. Radial artery pulse waveform analysis based on curve fitting using discrete Fourier series
Guo et al. Wrist pulse signal acquisition and analysis for disease diagnosis: A review
JP2019106202A (en) Health state evaluation support system and capillary vessel data acquisition method
Magenes et al. Comparison of data mining techniques applied to fetal heart rate parameters for the early identification of IUGR fetuses
Hong et al. Detection of physical stress using multispectral imaging
CN111839488A (en) Non-invasive continuous blood pressure measuring device and method based on pulse wave
KR101912215B1 (en) Device for evaluating vascular elasticity
CN115868944A (en) Method for recognizing slow pulse characteristics by pulse envelope detection of cunguanchi three-part pressure sensing array
Shah et al. Extraction of dynamical information and classification of heart rate variability signals using scale based permutation entropy measures.
CN113288075A (en) Floating, middle and deep pulse position identification method and system
Jezewski et al. Extraction of fetal heart-rate signal as the time event series from evenly sampled data acquired using Doppler ultrasound technique
CN108742549B (en) Image information generation method and pulse wave measurement system
US20220211283A1 (en) Methods for blood pressure calibration selection and modeling methods thereof
CN115381412B (en) Method for identifying chordal pulse characteristics by combining pulse wave array signal envelope diagram with time domain signal
CN113576422B (en) Method, system and equipment for simplifying processing and quantizing physiological signal waveform
CN115251856A (en) Pulse strength judgment method and device, electronic equipment and storage medium
KR102147568B1 (en) Apparatus and method for determining smoothness of pulse wave signal
US20200375550A1 (en) Method, apparatus and computer program for measuring and analysis of blood pressure using photoplethysmography
CN116959742B (en) Blood glucose data processing method and system based on spherical coordinate kernel principal component analysis
Umasha et al. High-Dimensional Feature Space for Diabetes Diagnosis and Identification of Diabetic-Sensitive Features in Ayurvedic Nadi Signals
Massmann et al. Continuous signal quality estimation for robust heart rate extraction from photoplethysmographic signals
US20170188892A1 (en) Method and Apparatus for Using Adaptive Plethysmographic Signal Conditioning to Determine Patient Identity
JP7137550B2 (en) Arrhythmia determination device and arrhythmia determination program

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant