US20140142453A1 - Method for extracting the feature of an abdominal breathing and a system using the same - Google Patents
Method for extracting the feature of an abdominal breathing and a system using the same Download PDFInfo
- Publication number
- US20140142453A1 US20140142453A1 US14/165,433 US201414165433A US2014142453A1 US 20140142453 A1 US20140142453 A1 US 20140142453A1 US 201414165433 A US201414165433 A US 201414165433A US 2014142453 A1 US2014142453 A1 US 2014142453A1
- Authority
- US
- United States
- Prior art keywords
- abdominal breathing
- function
- intrinsic mode
- feature
- functions
- 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.)
- Abandoned
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/08—Detecting, measuring or recording devices for evaluating the respiratory organs
- A61B5/0816—Measuring devices for examining respiratory frequency
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/08—Detecting, measuring or recording devices for evaluating the respiratory organs
- A61B5/087—Measuring breath flow
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
- A61B5/113—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb occurring during breathing
- A61B5/1135—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb occurring during breathing by monitoring thoracic expansion
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
- A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
- A61B5/6813—Specially adapted to be attached to a specific body part
- A61B5/6823—Trunk, e.g., chest, back, abdomen, hip
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/08—Detecting, measuring or recording devices for evaluating the respiratory organs
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
- A61B5/113—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb occurring during breathing
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Physiology (AREA)
- Pulmonology (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Dentistry (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
A method for extracting the feature of an abdominal breathing is disclosed, capable of extracting the feature of an abdominal breathing, without the requirement of a standard model of an abdominal breathing and the execution of a learning process being executed prior to the method for extracting the feature of an abdominal breathing. By means of computing a plurality of intrinsic mode functions corresponding to the abdominal breathing signal received, an Euler angle function and an instantaneous frequency function of each of the plurality of intrinsic mode functions, and comparing the plurality of instantaneous frequency function with a pre-determined zero-point threshold region, the method for extracting the feature of an abdominal breathing defines one of the plurality of instantaneous frequency function as an abdominal breathing feature function, which contains the feature of the abdominal breathing. In this way, the feature of an abdominal breathing is extracted.
Description
- 1. Field of the Invention
- The present invention relates to a method for extracting the feature of an abdominal breathing and an abdominal breathing feature extracting system, more particularly, to the method for extracting the feature of an abdominal breathing and the abdominal breathing feature extracting system capable of extracting the feature of an abdominal breathing, without the requirement of a standard model of an abdominal breathing and the execution of a learning process, which is executed prior to the method for extracting the feature of an abdominal breathing, letting the person to learn how to execute a standard abdominal breathing.
- 2. Description of Related Art
- In recent years, abdominal breathing plays an important role in rehabilitation and stress-released procedure. However, a person who needs to execute an abdominal breathing, such as a patient going to have a heart surgery, must be trained for a period of time by a teacher, such as a doctor, or a nurse. Besides, the results of the training, such as the correctness of the execution of the abdominal breathing, is not yet able to evaluated efficiently, as a standard model of an abdominal breathing is still required for the person to learn with. Moreover, the standard model is not necessarily suitable for full spectrum of people, such as from a children to an elderly.
- As mentioned above, a learning process must be executed, for the person to learn with the standard model, which will cost a lot of time, and a certain level of cost. And most important of all, the person must be the same place with the teacher teaching the abdominal breathing, raising a lot of difficulty for the person in transportation, especially for the elderly.
- As a result, a method for extracting the feature of an abdominal breathing and an abdominal breathing feature extracting system capable of extracting the feature of an abdominal breathing, without the requirement of a standard model of an abdominal breathing and the execution of a learning process, which is executed prior to the method for extracting the feature of an abdominal breathing, letting the person to learn how to execute a standard abdominal breathing.
- It is one object of the present invention is to provide a method for extracting the feature of an abdominal breathing, capable of extracting the feature of an abdominal breathing, without the requirement of a standard model of an abdominal breathing and the execution of a learning process being executed prior to the method for extracting the feature of an abdominal breathing.
- It is another object of the present invention is to provide an abdominal breathing feature extracting system, capable of extracting the feature of an abdominal breathing, without the requirement of a standard model of an abdominal breathing and the execution of a learning process being executed prior to the method for extracting the feature of an abdominal breathing.
- To achieve the object, the method for extracting the feature of an abdominal breathing of the present invention comprises the following steps of: receiving an abdominal breathing signal; executing an empirical mode decomposition process on the abdominal breathing signal, for computing a plurality of intrinsic mode functions; basing on the plurality of intrinsic mode functions and the results obtained after the execution of a Hilbert transform on the plurality of intrinsic mode functions, computing an Euler angle function of each of the plurality of intrinsic mode functions; taking the partial derivative of the Euler angle function of each of the plurality of intrinsic mode functions, with respect to time, for computing an instantaneous frequency function of each of the plurality of intrinsic mode functions, and extracting a plurality of maximum values of each of the plurality of instantaneous frequency functions; and according to a pre-determined sequence, comparing the plurality of maximum values of each of the plurality of instantaneous frequency functions with a zero-point threshold region, respectively, and when a result that all of the plurality of maximum values of one of the plurality of instantaneous frequency functions falling into the zero-point threshold region is obtained, the instantaneous frequency function having the plurality of maximum values is defined as an abdominal breathing feature function.
- To achieve the object, the abdominal breathing feature extracting system of the present invention comprises: a sensing module, for sensing an abdominal breathing signal; a computing module, coupled with the sensing module, for extracting an abdominal breathing feature function; and an abdominal breathing feature output module, coupled with the computing module, for outputting the abdominal breathing feature function; wherein, the computing module extracts the abdominal breathing feature function from the abdominal breathing signal by means of executing an abdominal breathing feature extracting method, and the abdominal breathing feature extracting method comprises the following steps of: receiving the abdominal breathing signal; executing an empirical mode decomposition process on the abdominal breathing signal, for computing a plurality of intrinsic mode functions; basing on the plurality of intrinsic mode functions and the results obtained after the execution of a Hilbert transform on the plurality of intrinsic mode functions, computing an Euler angle function of each of the plurality of intrinsic mode functions; taking the partial derivative of the Euler angle function of each of the plurality of intrinsic mode functions, with respect to time, for computing an instantaneous frequency function of each of the plurality of intrinsic mode functions, and extracting a plurality of maximum values of each of the plurality of instantaneous frequency functions; and according to a pre-determined sequence, comparing the plurality of maximum values of each of the plurality of instantaneous frequency functions with a zero-point threshold region, respectively, and when a result that all of the plurality of maximum values of one of the plurality of instantaneous frequency functions falling into the zero-point threshold region is obtained, the instantaneous frequency function having the plurality of maximum values is defined as the abdominal breathing feature function.
- Since the method for extracting the feature of an abdominal breathing can compute a plurality of intrinsic mode functions corresponding to an abdominal breathing signal received, and then compute an Euler angle function of each of the plurality of intrinsic mode functions. Later, an instantaneous frequency function of each of the plurality of intrinsic mode functions can be compute by taking the partial derivative of the Euler angle function of each of the plurality of intrinsic mode functions, with respect to time. Then, a plurality of maximum values of each of the plurality of instantaneous frequency functions can be extracted and compared with a zero-point threshold region. At final, with respect to the result of the comparison, an abdominal breathing feature function corresponding to the abdominal breathing signal received can be defined. Therefore, the method for extracting the feature of an abdominal breathing of the present invention is capable of extracting the feature of an abdominal breathing directly from the abdominal breathing signal received, without the requirement of a standard model of an abdominal breathing and the execution of a learning process being executed prior to the method for extracting the feature of an abdominal breathing.
- In addition, by including an computing module executing the method for extracting the feature of an abdominal breathing, a sensing module used for receiving the abdominal breathing signal, and an abdominal breathing feature output module used for outputting the abdominal breathing feature function, the abdominal breathing feature extracting system of the present invention is capable of extracting the feature of an abdominal breathing, without the requirement of a standard model of an abdominal breathing and the execution of a learning process being executed prior to the method for extracting the feature of an abdominal breathing.
- Other objects, advantages, and novel features of the invention will become more apparent from the following detailed description when taken in conjunction with the accompanying drawings.
-
FIG. 1 is a flowchart of the method for extracting the feature of an abdominal breathing according to first embodiment of the present invention. -
FIG. 2 is perspective view of an abdominal breathing feature extracting system. -
FIG. 3A displays the thoracic breathing signals acquired by a belt located on the chest of the person executing the thoracic breathing. -
FIG. 3B displays one of the plurality of intrinsic mode functions corresponding to the thoracic breathing signals displayed inFIG. 3A . -
FIG. 4A displays the abdominal breathing signals acquired by a belt located on the abdominal of the person executing the thoracic breathing. -
FIG. 4B displays one of the plurality of intrinsic mode functions corresponding to the abdominal breathing signals displayed inFIG. 4A . -
FIG. 5A displays the thoracic breathing signals acquired by a belt located on the chest of the person executing the abdominal breathingFIG. 5B displays one of the plurality of intrinsic mode functions corresponding to the thoracic breathing signals displayed inFIG. 5A . -
FIG. 6A displays the abdominal breathing signals acquired by a belt located on the abdominal of the person executing the abdominal breathing. -
FIG. 6B displays one of the plurality of intrinsic mode functions corresponding to the abdominal breathing signals displayed inFIG. 6A . -
FIG. 7A toFIG. 7E display 5 instantaneous frequency functions, of 5 of the plurality of intrinsic mode functions corresponding to the abdominal breathing signals displayed inFIG. 6A . -
FIG. 8 is a perspective view of the abdominal breathing feature extracting system according to second embodiment of the present invention. -
FIG. 9 is a flowchart of the method for extracting the feature of an abdominal breathing executed by the computing module of the abdominal breathing feature extracting system according to second embodiment of the present invention. - As shown in
FIG. 1 , which is a flowchart of the method for extracting the feature of an abdominal breathing according to first embodiment of the present invention, the method for extracting the feature of an abdominal breathing comprises the following steps of: - (A) receiving an abdominal breathing signal;
- (B) executing an empirical mode decomposition process on the abdominal breathing signal, for computing a plurality of intrinsic mode functions;
- (C) basing on the plurality of intrinsic mode functions and the results obtained after the execution of a Hilbert transform on the plurality of intrinsic mode functions, computing an Euler angle function of each of the plurality of intrinsic mode functions;
- (D) taking the partial derivative of the Euler angle function of each of the plurality of intrinsic mode functions, with respect to time, for computing an instantaneous frequency function of each of the plurality of intrinsic mode functions, and extracting a plurality of maximum values of each of the plurality of instantaneous frequency functions; and
- (E) according to a pre-determined sequence, comparing the plurality of maximum values of each of the plurality of instantaneous frequency functions with a zero-point threshold region, respectively, and when a result that all of the plurality of maximum values of one of the plurality of instantaneous frequency functions falling into the zero-point threshold region is obtained, the instantaneous frequency function having the plurality of maximum values is defined as an abdominal breathing feature function.
- In the following, the operation of each of the above steps of the method for extracting the feature of an abdominal breathing will be described in detail, in accompanied with figures.
- Please refer to
FIG. 2 , which is a perspective view of an abdominal breathing feature extracting system. A person, whose abdominal breathing features are going to be extracted, wears twobelts - In the period of time when the method for extracting the feature of an abdominal breathing is executed, both of the displacement of the chest and the displacement of the abdominal of the person are acquired by the two
belts belts portable devices 23, through any kind of possible signal transmission means, such as wire connection, for example, cables, or wireless connection. The computer and/orportable devices 23 store a program corresponding to the method for extracting the feature of an abdominal breathing according to first embodiment of the present invention, in the storage unit thereof. - Once the computer and/or
portable devices 23 of the abdominal breathing feature extracting system receives the abdominal breathing signals, which is the step (A) of the method for extracting the feature of an abdominal breathing according to first embodiment of the present invention, the step (B) will be executed immediately. As shown inFIG. 1 , in the step (B), an empirical mode decomposition process is executed on the acquired abdominal breathing signals, for computing a plurality of intrinsic mode functions of each of the abdominal breathing signals. - In the present embodiment, there are 2 breathing signals being acquired by the two
belts - For example, as shown in
FIG. 2 , once the breathing signals acquired by thebelt 21 having been received by the computer and/orportable devices 23, the computer and/orportable devices 23 execute the EMD process on the breathing signals, such as the one displayed inFIG. 3A , in order to obtain a plurality of intrinsic mode functions, so-called IMFs, corresponding to the thoracic breathing signals. In the case shown inFIG. 3A , there are 16 IMFs obtained, or “decomposed” from the thoracic breathing signals displayed inFIG. 3A . - After all of the IMFs corresponding to the thoracic breathing signals displayed in
FIG. 3A is obtained, only the 11th of these IMFs is displayed inFIG. 3B , as an example. It should be noticed that, the selection of which of these IMFs to be displayed inFIG. 3B is arbitrary, and the selection of these IMFs should not limit the scope of the present invention. - In the following 8 figures, i.e. through
FIG. 3A toFIG. 6B , the 4 breathing signals and the 4 IMFs obtained after the execution of the EMD process and the selection are respectively displayed. As described above, the thoracic breathing signals acquired by thebelt 21 located on the chest of the person executing the thoracic breathing is displayed inFIG. 3A . The IMF obtained after the execution of the EMD process on the breathing signals displayed inFIG. 3A and the selection is displayed inFIG. 3B . - In same manner, the abdominal breathing signals acquired by the
belt 22 located on the abdominal of the person executing the thoracic breathing is displayed inFIG. 4A , while the IMF obtained after the execution of the EMD process on the breathing signals displayed inFIG. 4A and the selection is displayed inFIG. 4B . Besides, the thoracic breathing signals acquired by thebelt 21 located on the chest of the person executing the abdominal breathing is displayed inFIG. 5A , while the IMF obtained after the execution of the EMD process on the breathing signals displayed inFIG. 5A and the selection is displayed inFIG. 5B . - In addition, the abdominal breathing signals extracted by the
belt 22 located on the abdominal of the person executing the abdominal breathing is displayed inFIG. 6A , while the IMF obtained after the execution of the EMD process on the breathing signals displayed inFIG. 6A and the selection is displayed inFIG. 6B . - Since the execution of the EMD process on a signal, such as the breathing signal of the present invention, is well-known in the field, for example, the non-linear and non-stationary data processing field, detailed description regarding the execution of the EMD process will be omitted hereinafter.
- After the plurality of IMFs of an abdominal breathing signal, such as the one displayed in
FIG. 6A , being obtained, the step (C) of the method for extracting the feature of an abdominal breathing according to first embodiment of the present invention is executed. In the present embodiment, there are 16 IMFs corresponding to the abdominal breathing signals displayed inFIG. 6A . However, the number of the plurality of IMFs is not limited, there could be 1 to 16 IMFs corresponding to the abdominal breathing signals, depending on the characteristic of the IMFs. In addition, each of the plurality of IMFs has a characteristic frequency, and the values of the characteristic frequency of different IMFs are different from each other. - As shown in
FIG. 1 , in the step (C), an Euler angle function of each of the plurality of IMFs is computed basing on the plurality of IMFs and the results obtained after the execution of a Hilbert transform on the plurality of IMFs. In the present embodiment, since there are 16 IMFs, 16 results will be obtained after the execution of a Hilbert transform on these 16 IMFs, respectively. Therefore, 16 Euler angle functions will be obtained after the execution of the step (C). - Since the execution of the Hilbert transform on a function is well-known in the field, such as the signal processing field, detailed description regarding the execution of the Hilbert transform will be omitted hereinafter.
- After executing the Hilbert transform on all of the 16 IMFs, 16 Euler angle functions will be obtained. In the present embodiment, the Euler angle function is a function describing the functional relation between the Euler angle and time in an Euler formula. The Euler formula can be represented as:
-
re jθ =r(cos θ+j sin θ) Formula (1) - where, r represents the magnitude, and θ represents the Euler angle.
- Therefore, since the Euler angle function is computed basing on an intrinsic mode function and the result obtained after the execution of a Hilbert transform on the IMFs, the Euler angle functions can be represented as:
-
r a e jθa =IMFsa +jH(IMFsa) Formula (2) - where, ra represents the magnitude of ath intrinsic mode function, θa represents the Euler angle of ath intrinsic mode function, and H represents the Hilbert transform.
- As described above, 16 Euler angle functions are obtained after the execution of the step (C). Then, in the step (D) of the method for extracting the feature of an abdominal breathing according to first embodiment of the present invention, an instantaneous frequency function of each of the plurality of intrinsic mode function is computed, by taking the partial derivative of the Euler angle function of each of the plurality of intrinsic mode functions, with respect to time. The instantaneous frequency functions can be represented as:
-
- where, ωa represents the instantaneous frequency functions of ath intrinsic mode function, and θa represents the Euler angle of ath intrinsic mode function.
- As described above, since there are 16 IMFs corresponding to the breathing signals acquired by the
belt 22 located on the abdominal of the person executing the abdominal breathing is displayed inFIG. 6A , 16 instantaneous frequency functions will be obtained after taking the partial derivative of the 16 Euler angle functions, with respect to time. Please refer toFIG. 7A toFIG. 7E , 5 of the 16 instantaneous frequency functions are displayed. It should be noticed that, the selection of these 5 instantaneous frequency functions to be displayed in these 5 figures is arbitrary, and the selection of these instantaneous frequency functions should not limit the scope of the present invention. - After all of these instantaneous frequency functions have been obtained, a plurality of maximum values of each of these instantaneous frequency functions, in the present embodiment, 16 of them, is extracted. Besides, as can be seen in these 5 figures, the plurality of maximum values of each of these 5 instantaneous frequency functions is local maximum values.
- At last, in the step (E) of the method for extracting the feature of an abdominal breathing according to first embodiment of the present invention, the plurality of maximum values of each of the plurality of instantaneous frequency functions is compared with a zero-point threshold region, respectively, according to a pre-determined sequence. Besides, in the present embodiment, the pre-determined sequence is a sequence staring from the IMFs having a higher characteristic frequency value thereof, and terminating with the IMFs having a lower characteristic frequency value. In other words, the plurality of instantaneous frequency functions is sequenced, according to the characteristic frequency value of the IMFs corresponding thereto.
- As a result, in the present embodiment, the plurality of maximum values of the instantaneous frequency function, corresponding to the IMFs having a highest characteristic frequency value, is compared with a zero-point threshold region first, i.e. the instantaneous frequency function displayed in
FIG. 7A . - As shown in
FIG. 7A , the zero-point threshold region 71 is a region where the value of the instantaneous frequency function is between −0.5 and 0.5. In addition, as the zero-point threshold region 71 is between −0.5 and 0.5, the zero-point threshold region 71 is also called as the zero-point threshold ±0.5 region. As clearly shown inFIG. 7A , most of the plurality of maximum values of the instantaneous frequency function is higher than the top-limit of the zero-point threshold region 71. That is, the result that all of the plurality of maximum values of the instantaneous frequency function falling into the zero-point threshold region cannot be obtained. - Then, turning to the next instantaneous frequency function in the pre-determined sequence, such as the one displayed in
FIG. 7B . As clearly shown in the figure, some of the plurality of maximum values of the instantaneous frequency function is higher than the top-limit of the zero-point threshold region 71. That is, the result that all of the plurality of maximum values of the instantaneous frequency function falling into the zero-point threshold region cannot be obtained. - Next, in the same manner, turn to the next instantaneous frequency function in the pre-determined sequence, such as the one displayed in
FIG. 7C . As clearly shown in the figure, the highest one of the plurality of maximum values of the instantaneous frequency function, i.e. the global maximum value, is higher than the top-limit of the zero-point threshold region 71. That is, the result that all of the plurality of maximum values of the instantaneous frequency function falling into the zero-point threshold region cannot be obtained. - Next, in the same manner, turn to the next instantaneous frequency function in the pre-determined sequence, such as the one displayed in
FIG. 7D . As clearly shown in the figure, the highest one of the plurality of maximum values of the instantaneous frequency function, i.e. the global maximum value raised around 93 seconds, is higher than the top-limit of the zero-point threshold region 71. That is, the result that all of the plurality of maximum values of the instantaneous frequency function falling into the zero-point threshold region cannot be obtained. - However, in the next instantaneous frequency function in the pre-determined sequence, such as the one displayed in
FIG. 7E , the highest one of the plurality of maximum values of the instantaneous frequency function is lower than the top-limit of the zero-point threshold region 71, and the lowest of the plurality of maximum values of the instantaneous frequency function is higher than the bottom-limit of the zero-point threshold region 71. That is, the result that all of the plurality of maximum values of the instantaneous frequency function falling into the zero-point threshold region is obtained. - At this time, i.e. when a result that all of the plurality of maximum values of one of the plurality of instantaneous frequency functions falling into the zero-point threshold region is obtained, the instantaneous frequency function having the plurality of maximum values is defined as an abdominal breathing feature function. In this embodiment, the instantaneous frequency function displayed in
FIG. 7E is defined as an abdominal breathing feature function, corresponding to the abdominal breathing signals displayed inFIG. 6A . Moreover, in the present embodiment, the distribution of the plurality of minimum values of the abdominal breathing feature function in the time dimension is the frequency of the abdominal breathing. - It should be noticed that, by executing the method for extracting the feature of an abdominal breathing according to first embodiment of the present invention, all of the abdominal breathing feature functions corresponding to other abdominal breathing signals can be respectively obtained, such as the abdominal breathing signal displayed in
FIG. 5A , the abdominal breathing signal displayed inFIG. 6A . - Please refer to
FIG. 8 , which is a perspective view of the abdominal breathing feature extracting system according to second embodiment of the present invention. As shown in the figure, the abdominal breathing feature extracting system according to second embodiment of the present invention comprises: asensing module 81, acomputing module 82 coupled with thesensing module 81, and an abdominal breathingfeature output module 83 coupled with thecomputing module 82. Wherein, thesensing module 81 is used for sensing an abdominal breathing signal, while thecomputing module 82 is used for extracting an abdominal breathing feature function. In addition, the abdominal breathingfeature output module 83 is used for outputting the abdominal breathing feature function. - In addition, the
computing module 82 extracts the abdominal breathing feature function from the abdominal breathing signal, by means of executing an abdominal breathing feature extracting method. As shown inFIG. 9 , the abdominal breathing feature extracting method comprises the following steps of: - (A) receiving the abdominal breathing signal;
- (B) executing an empirical mode decomposition process on the abdominal breathing signal, for computing a plurality of intrinsic mode functions;
- (C) basing on the plurality of intrinsic mode functions and the results obtained after the execution of a Hilbert transform on the plurality of intrinsic mode functions, computing an Euler angle function of each of the plurality of intrinsic mode functions;
- (D) taking the partial derivative of the Euler angle function of each of the plurality of intrinsic mode functions, with respect to time, for computing an instantaneous frequency function of each of the plurality of intrinsic mode functions, and extracting a plurality of maximum values of each of the plurality of instantaneous frequency functions; and
- (E) according to a pre-determined sequence, comparing the plurality of maximum values of each of the plurality of instantaneous frequency functions with a zero-point threshold region, respectively, and when a result that all of the plurality of maximum values of one of the plurality of instantaneous frequency functions falling into the zero-point threshold region is obtained, the instantaneous frequency function having the plurality of maximum values is defined as the abdominal breathing feature function.
- In the present embodiment, the
sensing module 81 of the abdominal breathing feature extracting system according to second embodiment of the present invention is integrated to a belt. For example, a belt for a person to wear on his/her body, preferably around the abdominal. In addition, as shown inFIG. 8 , thesensing modules 81 includes apiezoelectric unit 811 for extracting the displacement of the person's body, and an analog-to-digital convertingunit 812 electrically connected with the piezoelectric unit, in order to convert the analog signal output from the piezoelectric unit into a corresponding digital signal. Besides, for minimizing the size and weight of the sensing module, the sensing module is a system on chip (SoC). - In addition, as shown in
FIG. 8 , the abdominal breathingfeature output module 83 includes awireless transmission unit 831, for outputting the abdominal breathing feature function to a remote server (not shown in the figure) through a wireless connection, such as the 3G or the WiFi technology. - Since the abdominal breathing feature extracting method executed by the
computing module 82 of the abdominal breathing feature extracting system according to second embodiment of the present invention is the same abdominal breathing feature extracting method described in the first embodiment of the present invention, detailed description regarding the operation of the abdominal breathing feature extracting method is omitted hereinafter, for simplify the description of the present embodiment. - Although the present invention has been explained in relation to its preferred embodiment, it is to be understood that many other possible modifications and variations can be made without departing from the spirit and scope of the invention as hereinafter claimed.
Claims (16)
1. A method for extracting the feature of an abdominal breathing, comprising the following steps of:
receiving an abdominal breathing signal;
executing an empirical mode decomposition process on the abdominal breathing signal, for computing a plurality of intrinsic mode functions;
basing on the plurality of intrinsic mode functions and the results obtained after the execution of a Hilbert transform on the plurality of intrinsic mode functions, computing an Euler angle function of each of the plurality of intrinsic mode functions;
taking the partial derivative of the Euler angle function of each of the plurality of intrinsic mode functions, with respect to time, for computing an instantaneous frequency function of each of the plurality of intrinsic mode functions, and extracting a plurality of maximum values of each of the plurality of instantaneous frequency functions; and
according to a pre-determined sequence, comparing the plurality of maximum values of each of the plurality of instantaneous frequency functions with a zero-point threshold region, respectively, and when a result that all of the plurality of maximum values of one of the plurality of instantaneous frequency functions falling into the zero-point threshold region is obtained, the instantaneous frequency function having the plurality of maximum values is defined as an abdominal breathing feature function.
2. The method as claimed in claim 1 , wherein the Euler angle function is a function describing the functional relation between the Euler angle and time in an Euler formula, and the Euler formula is represented as:
re jθ =r(cos θ+j sin θ)
re jθ =r(cos θ+j sin θ)
where, r represents the magnitude, and θ represents the Euler angle.
3. The method as claimed in claim 1 , wherein the plurality of maximum values of each of the plurality of instantaneous frequency functions is local maximum values.
4. The method as claimed in claim 1 , wherein the number of the plurality of intrinsic mode functions is between 1 and 16, and each of the plurality of intrinsic mode functions has a characteristic frequency, and the values of the characteristic frequency of different intrinsic mode functions are different from each other.
5. The method as claimed in claim 4 , wherein the pre-determined sequence is a sequence staring from the intrinsic mode function having a higher characteristic frequency value thereof, and terminating with the intrinsic mode function having a lower characteristic frequency value.
6. The method as claimed in claim 1 , wherein the zero-point threshold region is a region where the value of the instantaneous frequency function is between −0.5 and 0.5.
7. The method as claimed in claim 1 , wherein the distribution of the plurality of minimum values of the abdominal breathing feature function in the time dimension is the frequency of the abdominal breathing.
8. An abdominal breathing feature extracting system, comprising:
a sensing module, for sensing an abdominal breathing signal;
a computing module, coupled with the sensing module, for extracting an abdominal breathing feature function; and
an abdominal breathing feature output module, coupled with the computing module, for outputting the abdominal breathing feature function;
wherein, the computing module extracts the abdominal breathing feature function from the abdominal breathing signal by means of executing an abdominal breathing feature extracting method, and the abdominal breathing feature extracting method comprises the following steps of:
receiving the abdominal breathing signal;
executing an empirical mode decomposition process on the abdominal breathing signal, for computing a plurality of intrinsic mode functions;
basing on the plurality of intrinsic mode functions and the results obtained after the execution of a Hilbert transform on the plurality of intrinsic mode functions, computing an Euler angle function of each of the plurality of intrinsic mode functions;
taking the partial derivative of the Euler angle function of each of the plurality of intrinsic mode functions, with respect to time, for computing an instantaneous frequency function of each of the plurality of intrinsic mode functions, and extracting a plurality of maximum values of each of the plurality of instantaneous frequency functions; and
according to a pre-determined sequence, comparing the plurality of maximum values of each of the plurality of instantaneous frequency functions with a zero-point threshold region, respectively, and when a result that all of the plurality of maximum values of one of the plurality of instantaneous frequency functions falling into the zero-point threshold region is obtained, the instantaneous frequency function having the plurality of maximum values is defined as the abdominal breathing feature function.
9. The system as claimed in claim 8 , wherein the sensing module is integrated to a belt, and the sensing module includes a piezoelectric unit and an analog-to-digital converting unit.
10. The system as claimed in claim 8 , wherein the sensing module is a system on chip.
11. The system as claimed in claim 8 , wherein the abdominal breathing feature output module includes a wireless transmission unit.
12. The system as claimed in claim 8 , wherein the Euler angle function is a function describing the functional relation between the Euler angle and time in an Euler formula, and the Euler formula is represented as:
re jθ =r(cos θ+j sin θ)
re jθ =r(cos θ+j sin θ)
where, r represents the magnitude, and θ represents the Euler angle.
13. The system as claimed in claim 8 , wherein the number of the plurality of intrinsic mode functions is between 1 and 16, and each of the plurality of intrinsic mode functions has a characteristic frequency, and the values of the characteristic frequency of different intrinsic mode functions are different from each other.
14. The system as claimed in claim 8 , wherein the pre-determined sequence is a sequence staring from the intrinsic mode function having a higher characteristic frequency value thereof, and terminating with the intrinsic mode function having a lower characteristic frequency value.
15. The system as claimed in claim 8 , wherein the zero-point threshold region is a region where the value of the instantaneous frequency function is between −0.5 and 0.5.
16. The system as claimed in claim 8 , wherein the distribution of the plurality of minimum values of the abdominal breathing feature function in the time dimension is the frequency of the abdominal breathing.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US14/165,433 US20140142453A1 (en) | 2011-08-08 | 2014-01-27 | Method for extracting the feature of an abdominal breathing and a system using the same |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US13/204,956 US8834385B2 (en) | 2011-08-08 | 2011-08-08 | Method for extracting the feature of an abdominal breathing and a system using the same |
US14/165,433 US20140142453A1 (en) | 2011-08-08 | 2014-01-27 | Method for extracting the feature of an abdominal breathing and a system using the same |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US13/204,956 Division US8834385B2 (en) | 2011-08-08 | 2011-08-08 | Method for extracting the feature of an abdominal breathing and a system using the same |
Publications (1)
Publication Number | Publication Date |
---|---|
US20140142453A1 true US20140142453A1 (en) | 2014-05-22 |
Family
ID=47677963
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US13/204,956 Expired - Fee Related US8834385B2 (en) | 2011-08-08 | 2011-08-08 | Method for extracting the feature of an abdominal breathing and a system using the same |
US14/165,433 Abandoned US20140142453A1 (en) | 2011-08-08 | 2014-01-27 | Method for extracting the feature of an abdominal breathing and a system using the same |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US13/204,956 Expired - Fee Related US8834385B2 (en) | 2011-08-08 | 2011-08-08 | Method for extracting the feature of an abdominal breathing and a system using the same |
Country Status (1)
Country | Link |
---|---|
US (2) | US8834385B2 (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10357163B1 (en) * | 2012-06-01 | 2019-07-23 | Vital Connect, Inc. | Respiratory rate detection using decomposition of ECG |
JP6140562B2 (en) * | 2013-07-24 | 2017-05-31 | 株式会社槌屋 | Respiration measurement method and respiration measurement device |
WO2015049150A1 (en) * | 2013-10-01 | 2015-04-09 | Koninklijke Philips N.V. | Improved signal selection for obtaining a remote photoplethysmographic waveform |
CN106108904B (en) * | 2016-06-23 | 2019-06-07 | 华中科技大学 | A kind of contactless human body respiration parameter method for real-time measurement and system |
CN111248922B (en) * | 2020-02-11 | 2022-05-17 | 中国科学院半导体研究所 | Human body respiration condition acquisition paste based on accelerometer and gyroscope and preparation method thereof |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4827943A (en) * | 1986-09-23 | 1989-05-09 | Advanced Medical Technologies, Inc. | Portable, multi-channel, physiological data monitoring system |
US6241683B1 (en) * | 1998-02-20 | 2001-06-05 | INSTITUT DE RECHERCHES CLINIQUES DE MONTRéAL (IRCM) | Phonospirometry for non-invasive monitoring of respiration |
US6893404B2 (en) * | 2000-04-07 | 2005-05-17 | Remo Ehf | Breathing movement measurements and apparatus |
US20060041201A1 (en) * | 2004-08-23 | 2006-02-23 | Khosrow Behbehani | System, software, and method for detection of sleep-disordered breathing using an electrocardiogram |
US7041062B2 (en) * | 2000-09-15 | 2006-05-09 | Friendly Sensors Ag | Device and method for producing respiration-related data |
US20070032733A1 (en) * | 2004-01-16 | 2007-02-08 | David Burton | Method and apparatus for ECG-derived sleep disordered breathing monitoring, detection and classification |
US7674230B2 (en) * | 2004-11-22 | 2010-03-09 | Widemed Ltd. | Sleep monitoring using a photoplethysmograph |
US20100130873A1 (en) * | 2008-04-03 | 2010-05-27 | Kai Sensors, Inc. | Non-contact physiologic motion sensors and methods for use |
US20100152600A1 (en) * | 2008-04-03 | 2010-06-17 | Kai Sensors, Inc. | Non-contact physiologic motion sensors and methods for use |
US20100204601A1 (en) * | 2009-02-10 | 2010-08-12 | Tanita Corporation | Respiration type evaluation apparatus |
-
2011
- 2011-08-08 US US13/204,956 patent/US8834385B2/en not_active Expired - Fee Related
-
2014
- 2014-01-27 US US14/165,433 patent/US20140142453A1/en not_active Abandoned
Patent Citations (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4827943A (en) * | 1986-09-23 | 1989-05-09 | Advanced Medical Technologies, Inc. | Portable, multi-channel, physiological data monitoring system |
US6241683B1 (en) * | 1998-02-20 | 2001-06-05 | INSTITUT DE RECHERCHES CLINIQUES DE MONTRéAL (IRCM) | Phonospirometry for non-invasive monitoring of respiration |
US6893404B2 (en) * | 2000-04-07 | 2005-05-17 | Remo Ehf | Breathing movement measurements and apparatus |
US7041062B2 (en) * | 2000-09-15 | 2006-05-09 | Friendly Sensors Ag | Device and method for producing respiration-related data |
US20070032733A1 (en) * | 2004-01-16 | 2007-02-08 | David Burton | Method and apparatus for ECG-derived sleep disordered breathing monitoring, detection and classification |
US7801593B2 (en) * | 2004-08-23 | 2010-09-21 | Board Of Regents, The University Of Texas System | System, software, and method for detection of sleep-disordered breathing using an electrocardiogram |
US7343198B2 (en) * | 2004-08-23 | 2008-03-11 | The University Of Texas At Arlington | System, software, and method for detection of sleep-disordered breathing using an electrocardiogram |
US20080319326A1 (en) * | 2004-08-23 | 2008-12-25 | The University Of Texas At Arlington | System, software, and method for detection of sleep-disordered breathing using an eltrocardiogram |
US20060041201A1 (en) * | 2004-08-23 | 2006-02-23 | Khosrow Behbehani | System, software, and method for detection of sleep-disordered breathing using an electrocardiogram |
US7674230B2 (en) * | 2004-11-22 | 2010-03-09 | Widemed Ltd. | Sleep monitoring using a photoplethysmograph |
US20100130873A1 (en) * | 2008-04-03 | 2010-05-27 | Kai Sensors, Inc. | Non-contact physiologic motion sensors and methods for use |
US20100152600A1 (en) * | 2008-04-03 | 2010-06-17 | Kai Sensors, Inc. | Non-contact physiologic motion sensors and methods for use |
US20100240999A1 (en) * | 2008-04-03 | 2010-09-23 | Kai Medical, Inc. | Systems and methods for point in time measurement of physiologic motion |
US20100249630A1 (en) * | 2008-04-03 | 2010-09-30 | Kai Medical, Inc. | Systems and methods for respiratory rate measurement |
US20100249633A1 (en) * | 2008-04-03 | 2010-09-30 | Kai Medical, Inc. | Systems and methods for determining regularity of respiration |
US20100292568A1 (en) * | 2008-04-03 | 2010-11-18 | Kai Medical, Inc. | Systems and methods for measurement of depth of breath and paradoxical breathing |
US8454528B2 (en) * | 2008-04-03 | 2013-06-04 | Kai Medical, Inc. | Non-contact physiologic motion sensors and methods for use |
US20100204601A1 (en) * | 2009-02-10 | 2010-08-12 | Tanita Corporation | Respiration type evaluation apparatus |
Also Published As
Publication number | Publication date |
---|---|
US8834385B2 (en) | 2014-09-16 |
US20130041277A1 (en) | 2013-02-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20140142453A1 (en) | Method for extracting the feature of an abdominal breathing and a system using the same | |
US11389083B2 (en) | Method, device and system for sensing neuromuscular, physiological, biomechanical, and musculoskeletal activity | |
US10993651B2 (en) | Exercise guidance method and exercise guidance device | |
Furness et al. | Acute injuries in recreational and competitive surfers: incidence, severity, location, type, and mechanism | |
O’Donnell et al. | The validity and reliability of the GymAware linear position transducer for measuring counter-movement jump performance in female athletes | |
US20150201867A1 (en) | Electronic free-space motion monitoring and assessments | |
US20170105682A1 (en) | Heart rate monitor with time varying linear filtering | |
US8113843B2 (en) | Apparatus for analyzing operations and method for analyzing operations | |
US20200121987A1 (en) | Smart gym | |
González-Villanueva et al. | Design of a wearable sensing system for human motion monitoring in physical rehabilitation | |
US11763696B2 (en) | Systems and methods for facilitating mind-body-emotion state self-adjustment and functional skills development by way of biofeedback and environmental monitoring | |
US10271776B2 (en) | Computer aided analysis and monitoring of mobility abnormalities in human patients | |
Lee et al. | A smartphone-centric system for the range of motion assessment in stroke patients | |
Louro et al. | Stability of patterns of behavior in the butterfly technique of the elite swimmers | |
Montes et al. | Using Hexoskin wearable technology to obtain body metrics during trail hiking | |
US20240122499A1 (en) | Wearable inertial sensor system and methods | |
US11705748B2 (en) | Wearable gesture recognition device for medical screening and associated operation method and system | |
Trajković et al. | Interrater and Intrarater Reliability of the EasyForce Dynamometer for Assessment of Maximal Shoulder, Knee and Hip Strength | |
TWI419675B (en) | Method for extracting the feature of an abdominal breathing and a system using the same | |
US20160213286A1 (en) | Electronic device, monitoring and feedback system on thoracoabdominal motion and method thereof | |
Romero-Morales et al. | Ankle Taping Effectiveness for the Decreasing Dorsiflexion Range of Motion in Elite Soccer and Basketball Players U18 in a Single Training Session: A Cross-Sectional Pilot Study | |
CN113241148A (en) | Fitness scheme generation method and device, terminal equipment and storage medium | |
JP5373010B2 (en) | Method for extracting abdominal breathing feature function and system using this method | |
CN111353345B (en) | Method, apparatus, system, electronic device, and storage medium for providing training feedback | |
Martínez-Martí et al. | Validation of instrumented insoles for measuring height in vertical jump |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |