US20080146901A1 - Optical measurement instrument for living body - Google Patents
Optical measurement instrument for living body Download PDFInfo
- Publication number
- US20080146901A1 US20080146901A1 US11/954,627 US95462707A US2008146901A1 US 20080146901 A1 US20080146901 A1 US 20080146901A1 US 95462707 A US95462707 A US 95462707A US 2008146901 A1 US2008146901 A1 US 2008146901A1
- Authority
- US
- United States
- Prior art keywords
- signal
- separation
- calculation unit
- signals
- living body
- 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
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/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/02028—Determining haemodynamic parameters not otherwise provided for, e.g. cardiac contractility or left ventricular ejection fraction
-
- 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/6814—Head
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
- A61B5/7207—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2134—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on separation criteria, e.g. independent component analysis
- G06F18/21343—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on separation criteria, e.g. independent component analysis using decorrelation or non-stationarity, e.g. minimising lagged cross-correlations
Definitions
- the present invention relates to an optical measurement technology, and particularly to an optical measurement method and instrument for receiving optical signals from a living body to obtain information on the living body. More particularly, the present invention relates to removing noise signals superimposed on the information on the living body to highly precisely obtain desired living body signals.
- noise signals which are caused by several factors, included in measurement signals that are obtained by measuring information inside the living body using light (e.g., Patent Document 1).
- An averaging method, a band-pass filter, and the like are used to reduce the effects from the noise signals.
- a technique using the averaging method reduces the noise signals based on a precondition that a brain activity responds to the same assignment in the same manner.
- a technique using the band-pass filter reduces the noise signals based on a precondition that a brain activity signal and a noise signal exist in different frequency bands. These techniques process a plurality of noise signals. However, a technique is also proposed that is specialized in processing noise signals caused by special factors.
- noise signals caused by pulsations can be reduced by means of a band-cut filter, because their frequency is readily identified.
- a technique is proposed that reduces pulsation noise signals included in optical measurement signals based on the pulsation signals measured at an area, such as an ear (Patent Document 2).
- a technique is also proposed that reduces effects from the pulsation noise signals by extracting the pulsation noise signals from optical measurement signals itself (Non-Patent Document 1).
- Non-Patent Document 1 Maria Angela Franceschini et al., NeuroImage 21 (2004) 372-386
- Non-Patent Document 2 Sato et al., NeuroImage 21 (2004) 1554-1562
- Non-Patent Document 3 Hirosaka et al., Physica D 194 (2004), Pages 320-332
- a subject being tested When a subject being tested is given a certain task in order to obtain a desired biological signal, for example, when the subject being tested is asked to play a simple voluntary movement at a given time, a displacement in an area being tested caused by the voluntary movement could lead to the noise.
- the noise signals superimposed on a desired biological signal very likely show changes similar to those of the desired biological signal such as a brain activity signal.
- the obtained signals can be erroneously interpreted.
- a noise associated with a body motion, or a body motion noise is among them.
- the body motion noise synchronizes with the brain activity, the targeted signal. Therefore, it is difficult to separate them. However, it is absolutely necessary to discriminate them.
- the present invention is implemented by an optical measurement instrument for a living body that comprises: a plurality of light irradiation means for irradiating an object being tested with light; a plurality of light reception means for detecting the light that is irradiated to the object being tested and propagates inside the object being tested; and a calculation unit for calculating hemodynamics of a plurality of factors inside the object being tested from the detection signal detected by the light reception means, wherein the calculation unit comprises: a separation calculation unit for separating a plurality of signals from the detection signal; a selection unit for calculating a correlation coefficient between a plurality of the signals separated by the separation calculation unit to select a signal in which the correlation coefficient is a predetermined value or greater as a noise signal; and a reconstruction calculation unit for reconfiguring the detection signal excluding the noise signal selected by the selection unit.
- a plurality of factors can be detected from optical signals measured through the use of a plurality of light sources. For example, lights of 690 nm and 830 nm can be used to detect oxygenated hemoglobin and deoxygenated hemoglobin (Non-Patent Document 2).
- the separation calculation unit calculates a plurality independent components (separation signals), and mixing coefficient matrices that indicate how each of the independent components contributes to a plurality of measurement signals by means of, for example, an independent component analysis method (Non-Patent Document 3).
- the independent component analysis method can recover original measurement signals from the independent components and mixture coefficient matrices.
- the selection unit selects separation signals based on a certain selection criterion. For example, when trying to select noise components from the separation signals separated by the independent component analysis method, the selection unit obtains a correlation coefficient between a plurality of measurement factors for each of the mixture coefficient matrices corresponding to each separation signal that is binarized by the positive and negative of the value, and selects a separation signal indicating a correlation coefficient greater than a certain value as a noise component.
- the reconstruction calculation unit uses only the separation signals selected by the selection unit or other separation signals to reconfigures signals. For example, when the separation signals obtained by the independent component analysis are selected as the noise components, the reconstruction calculation unit performs configuration by setting the mixture coefficient matrices of the selected separation components to zero. This allows the noises to be removed.
- the present invention allows the noise components, such as the body motion noises that are difficult to be separated and removed, to be separated and removed, and provides signals with few noise components.
- FIG. 1 is a schematic diagram of an embodiment of the present invention
- FIG. 2 is an exemplary interface unit in the embodiment
- FIG. 3 is exemplary signals measured by the interface unit of FIG. 2 ;
- FIG. 4 is exemplary separations signals calculated by a separation calculation unit
- FIG. 5 is exemplary mixing coefficient matrices calculated by the separation calculation unit
- FIG. 6 is an exemplary display of the separation signals as a result of selection by a selection unit using the mixing coefficient matrices
- FIG. 7 is an example in which the mixing coefficient matrices calculated by the separation calculation unit are binarized.
- FIG. 8 is an exemplary display of the separation signals as a result of selection by the selection unit using the binarized mixing coefficient matrices.
- FIG. 1 A specific block diagram is shown in FIG. 1 as an embodiment.
- An interface unit ( 110 ) for optical measurement of a living body is attached to part of or the whole of the head of a subject being tested.
- a light irradiation unit ( 101 ) irradiates the living body with mixed lights of 690 nm and 830 nm wavelengths through an optical fiber ( 111 ) coupled to the interface unit ( 110 ).
- the lights of each wavelength, which passed through the living body via the optical fiber ( 111 ) coupled to the interface unit ( 110 ) are detected by a light detection unit ( 102 ), and the detection result is stored in a storage unit ( 103 ).
- light of other wavelengths may also be irradiated, and the number of combination of wavelengths employed may be three or more.
- the obtained optical signals may be employed to the processing described in the following sections, after they are subjected to some sort of calculation processing.
- the optical signals are separated by a separation calculation unit ( 121 ) based on the method and setting selected by a user.
- the signals separated by the separation calculation unit (referred to as separation signals) are manually or automatically selected by a selection unit ( 122 ) based on the method and setting selected by the user.
- the signals selected by the selection unit ( 122 ) are reconfigured by a reconstruction calculation unit ( 123 ) based on the method and setting selected by the user.
- a parameter input and a result display at the separation calculation unit, selection unit, and reconstruction calculation unit are performed by a display unit.
- the parameter refers to, for example, a setting value for the signal separation used by the separation calculation unit, and a threshold of the correlation coefficient used by the selection unit.
- Measure points at this time include a measure point 1 ( 200 a ), a measure point 2 ( 200 b ), a measure point 3 ( 200 c ), and a measure point 4 ( 200 d ), totaling four points.
- the measure point refers to a middle point between an irradiation point and a detection point.
- the present embodiment provides an example in which concentration changes in oxygenated hemoglobin and deoxygenated hemoglobin are calculated from optical signals (Non-Patent Document 2). Examples of signals calculated from the obtained optical signals are shown in FIG. 3 .
- FIG. 1 a measure point 1
- 200 b a measure point 2
- 200 c a measure point 3
- 200 d measure point 4
- FIG. 3 shows oxygenated hemoglobin concentration changes ( 201 to 204 ) measured at four measure points, and deoxygenated hemoglobin concentration changes ( 301 to 304 ) measured at the same measure points.
- a lateral axis indicates time (second), while a vertical axis indicates hemoglobin concentration change amount (mM-mm).
- oxygenated hemoglobin concentration changes at a measure point 1 202 oxygenated hemoglobin concentration changes at a measure point 2 203 oxygenated hemoglobin concentration changes at a measure point 3 204 oxygenated hemoglobin concentration changes at a measure point 4 301 deoxygenated hemoglobin concentration changes at a measure point 1 302 deoxygenated hemoglobin concentration changes at a measure point 2 303 deoxygenated hemoglobin concentration changes at a measure point 3 304 deoxygenated hemoglobin concentration changes at a measure point 4
- a TDDICA method (Non-Patent Document 3) may be used as one of the separation techniques. This is one of the techniques referred to as the independent component analysis.
- a measurement signal X comprised of time series data simultaneously obtained at respective four measure points during the measurement is expressed as follows:
- T indicates the number of data.
- a set of signals comprised of four statistically independent time series signal components (original signal S) is expressed as follows:
- the mixing coefficient matrix A and original signal S are estimated from only the independent condition of the original signal component and the measurement signal x.
- the number of samplings refers to the number of data in the time direction.
- the estimation is performed by whitening and rotation.
- the whitening means conversion into a linear independent signal by mixture of the measurement signals. It is given using the covariance matrix C of the measurement signal X as follows:
- This rotation matrix is determined by the above independence condition.
- it is difficult to achieve this. Therefore, minimizing the size of correlation is set as a condition instead.
- the noise applied to each sensor may b handled as the positive.
- following expression is used in which an observation noise is added to the right side thereof.
- pseudo whitening is performed instead of whitening, and the pseudo whitening is performed in the following manner using covariance matrix G.
- the matrix G must be estimated separately. It is possible to estimate the matrix G from a residual when one section in a rest state is fitted by a polynomial expression. A description has been made as to a case in which the number of measurement points is four. However, the separation calculation can be conducted by the same procedure even when the number of measurement points differs from that.
- FIG. 4 An example of separation signals is shown in FIG. 4. 211 to 214 are separation signals calculated from the oxygenated hemoglobin concentration changes obtained at four measure points. 311 to 314 are separation signals calculated from deoxygenated hemoglobin concentration changes obtained at four measure points. Description of reference numerals used in FIG. 4 is as follows.
- separation signal 1 of oxygenated hemoglobin concentration changes 212 separation signal 2 of oxygenated hemoglobin concentration changes 213 separation signal 3 of oxygenated hemoglobin concentration changes 214 separation signal 4 of oxygenated hemoglobin concentration changes 311 separation signal 1 of deoxygenated hemoglobin concentration changes 312 separation signal 2 of deoxygenated hemoglobin concentration changes 313 separation signal 3 of deoxygenated hemoglobin concentration changes 314 separation signal 4 of deoxygenated hemoglobin concentration changes
- FIG. 5 An example of mixing coefficient matrices is shown in FIG. 5.
- 231 to 234 indicate the mixing coefficient matrices 211 to 214 , respectively. Color shading indicates the size of the coefficient. Description of reference numerals used in FIG. 5 is as follows.
- mixing coefficient matrix of the separation signal 1 of oxygenated hemoglobin concentration changes 232 mixing coefficient matrix of the separation signal 2 of oxygenated hemoglobin concentration changes 233 mixing coefficient matrix of the separation signal 3 of oxygenated hemoglobin concentration changes 234 mixing coefficient matrix of the separation signal 4 of oxygenated hemoglobin concentration changes 234 a mixing coefficient matrix at a measure point 1 of the separation signal 1 of oxygenated hemoglobin concentration changes 234 b mixing coefficient matrix at a measure point 2 of the separation signal 1 of oxygenated hemoglobin concentration changes 234 c mixing coefficient matrix at a measure point 3 of the separation signal 1 of oxygenated hemoglobin concentration changes 234 d mixing coefficient matrix at a measure point 4 of the separation signal 1 of oxygenated hemoglobin concentration changes 331 mixing coefficient matrix of the separation signal 1 of deoxygenated hemoglobin concentration changes 332 mixing coefficient matrix of the separation signal 2 of deoxygenated hemoglobin concentration changes 333 mixing coefficient matrix of the separation signal 3 of deoxygenated hemoglobin concentration changes 334 mixing coefficient matrix of the separation signal 4 of deoxygenated hemoglobin concentration
- the mixing coefficient matrix 234 of the separation signal 214 will be described in more detail.
- 234 a to 234 d corresponds to the measure points 1 to 4 , respectively. More specifically, it is indicated that the separation signal 214 is strongly mixed at the measure point 2 ( 234 b ), and hardly mixed at the measure point 3 ( 234 c ).
- filtering processing may be performed before the separation calculation.
- a filtering processing is effective that reduces frequency band differing from the frequency band.
- the noise is referred to as a body motion noise hereinafter.
- a body motion noise it is a noise that occurs when a body moves. While the occurrence mechanism of the body motion noise is not known sufficiently yet, it is thought that the body motion noise occurs due to a change in the contact state between a measurement interface and a subject to be measured caused by the body motion, a dynamic change in the internal state of the subject to be measured, and the like.
- the signal of brain activity evoked by body motion is a target signal
- it is extremely difficult to separate the target signal from a body motion signal because the body motion signal occurs in synchronization with the brain activity as the target signal.
- similarity is not guaranteed between different measurement factors (here, oxygenated hemoglobin and deoxygenated hemoglobin). Therefore, it is difficult to select the body motion noise from the correlation of the calculated separation signals.
- the correlation between the oxygenated hemoglobin and deoxygenated hemoglobin is calculated from the mixing coefficient matrices shown in FIG. 5 , and an example is shown in FIG. 6 in which a correlation coefficient with correlation higher than a certain value (here 0.5) is represented in black, while that with correlation lower than the certain value is represented in white.
- a correlation coefficient with correlation higher than a certain value here 0.5
- that with correlation lower than the certain value is represented in white.
- column 421 to 424 correspond to the mixing coefficient matrices 231 to 234 , respectively, and rows 431 to 434 correspond to mixing coefficient matrices 331 to 334 , respectively.
- 441 indicates that the correlation coefficient between the mixing coefficient matrix 231 and mixing coefficient matrix 331 is greater than 0.5
- 442 indicates that the correlation between the mixing coefficient matrix 232 and mixing coefficient matrix 332 is greater than 0.5.
- the result indicates that the separation components 211 , 311 and separation components 212 , 312 are body motion noises.
- the mixing coefficient matrix has a high value at a particular measure point (measure point 3 for 233 and 333 in FIG. 5 ) and has a low value at other measure points. Similarity of each value between oxygenated hemoglobin and deoxygenated hemoglobin is thought to be low at the measure point where the value is low because of the complexity of the living body.
- binarization using the positive and negative of the value is performed. More specifically, binarization is performed using zero as a threshold value.
- An example of binarized mixing coefficient matrices is shown in FIG. 7. 241 to 244 and 341 to 344 show 231 to 234 and 331 to 334 of FIG. 5 that are binarized, respectively.
- Correlation coefficients are calculated from these binarized mixing coefficient matrices, and an example of display is shown in FIG. 8 in which correlation coefficients that are greater than a certain threshold (0.5 here) are selected.
- columns 521 to 524 correspond to binarized mixing coefficient matrices 241 to 244 , respectively, and rows 531 to 534 correspond to mixing coefficient matrices 341 to 344 , respectively.
- 541 indicates that the correlation coefficient between the binarized mixing coefficient matrix 241 and mixing coefficient matrix 341 is greater than a certain value
- 542 indicates that the correlation coefficient between the mixing coefficient matrix 242 and mixing coefficient matrix 342 is smaller than a certain value.
- 442 is selected as a body motion noise.
- 542 is not selected as the body motion noise. Accordingly, separation signals 212 and 312 are thought to have corresponded to the brain activity.
- the user may select a body motion noise based on this result or may automatically select it.
- 0.5 is used here as a threshold of the correlation coefficient, other values may be used depending on the measurement signals.
- selection is made from among all of the separation signals.
- the selection unit may extract separation components having high correlation with the time for the assignment in advance to make selection from among the previously extracted separation components.
- synchronism of signal change with the assignment that is performed repeatedly may be used. Specifically, when the assignment is repeated, for example, for five times, measurement signals are extracted at times in synchronization therewith, and a correlation coefficient of the extracted five measurement signal fragments is used as an indicator of the synchronism.
- the similarity in binarized mixing coefficient matrix between oxygenated hemoglobin and deoxygenated hemoglobin is used as a selection criterion in the present embodiment
- the difference in phase between oxygenated hemoglobin changes and deoxygenated hemoglobin changes during brain activity may be used as the selection criterion.
- an instantaneous phase is calculated through Hilbert transform or the like, and when the phase difference between oxygenated hemoglobin and deoxygenated hemoglobin is smaller than a certain value, it may be determined that the body motion noise exists, while when the phase difference is in the range of a certain value different from the foregoing value, it may be determined that brain activity exists.
- the separation calculation and signal selection are performed by use of the whole measurement signals in the present embodiment.
- the separation calculation and signal selection may also be performed by dividing the measurement signals into parts.
- the signal selection may be performed using the result of the separation calculation that is performed for part of a plurality of the measurement signals.
- the obtained result of a plurality of separation calculations may be grouped by a method, such as clustering, and subsequently each group may be subjected to the selection processing.
- the reconstruction calculation unit ( 123 ) performs signal reconstruction based on the signal selection performed by the selection unit. Specifically, when the selection is for the purpose of removing noises, zero is substituted into the mixing coefficient matrix of the separation signal, and reconstruction is performed by a reconstruction method using the independent component analysis method. When the selection is for the purpose of extracting a target signal, zero is substituted into the mixing coefficient matrix of an unselected separation signal, and the reconstruction is performed by the reconstruction method using the independent reconstruction analysis method. It should be noted that although 0 is substituted into the mixing coefficient matrix of the separation signal for the purpose of removing noises in the above, a value that is greater than 0 and smaller than 1 may be multiplied.
- the present invention can be employed for the optical measurement instrument for a living body for measuring the living body.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- General Health & Medical Sciences (AREA)
- Surgery (AREA)
- Molecular Biology (AREA)
- Medical Informatics (AREA)
- Heart & Thoracic Surgery (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Cardiology (AREA)
- Physiology (AREA)
- Artificial Intelligence (AREA)
- Signal Processing (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Psychiatry (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
Description
- The present application claims priority from Japanese application JP2006-337737 filed on Dec. 15, 2006, the content of which is hereby incorporated by reference into this application.
- The present invention relates to an optical measurement technology, and particularly to an optical measurement method and instrument for receiving optical signals from a living body to obtain information on the living body. More particularly, the present invention relates to removing noise signals superimposed on the information on the living body to highly precisely obtain desired living body signals.
- Besides targeted signals, there are also noise signals, which are caused by several factors, included in measurement signals that are obtained by measuring information inside the living body using light (e.g., Patent Document 1). An averaging method, a band-pass filter, and the like are used to reduce the effects from the noise signals.
- A technique using the averaging method reduces the noise signals based on a precondition that a brain activity responds to the same assignment in the same manner.
- A technique using the band-pass filter reduces the noise signals based on a precondition that a brain activity signal and a noise signal exist in different frequency bands. These techniques process a plurality of noise signals. However, a technique is also proposed that is specialized in processing noise signals caused by special factors.
- For example, noise signals caused by pulsations (pulsation noise signals) can be reduced by means of a band-cut filter, because their frequency is readily identified. A technique is proposed that reduces pulsation noise signals included in optical measurement signals based on the pulsation signals measured at an area, such as an ear (Patent Document 2). A technique is also proposed that reduces effects from the pulsation noise signals by extracting the pulsation noise signals from optical measurement signals itself (Non-Patent Document 1).
- Moreover, a technique is also proposed that divides optical measurement signals measured at many points into a plurality of signals by a signal analysis method, and extracts noise signals and target signals (Patent document 3).
- However, when the noise signals show a similar variation as the target biological signals, conventional techniques have not been capable of removing the noise signals.
- When a subject being tested is given a certain task in order to obtain a desired biological signal, for example, when the subject being tested is asked to play a simple voluntary movement at a given time, a displacement in an area being tested caused by the voluntary movement could lead to the noise. At this time, the noise signals superimposed on a desired biological signal very likely show changes similar to those of the desired biological signal such as a brain activity signal.
- In such a case, the obtained signals can be erroneously interpreted. For example, a noise associated with a body motion, or a body motion noise is among them.
- While the occurrence mechanism of the body motion noise is not sufficiently known, changes in the contact state between a measurement interface and the subject being tested caused by the body motion, dynamic changes in the state inside the measurement target caused by the body motion, and the like are conceivable.
- When the brain activity induced by the voluntary movement is a targeted signal, the body motion noise synchronizes with the brain activity, the targeted signal. Therefore, it is difficult to separate them. However, it is absolutely necessary to discriminate them.
- It should be noted that although it has become likely that a signal separation method, such as an independent component analysis method, can separate these signals, a clear discrimination reference based on which to discriminate between the body motion noise and brain activity is not provided yet.
- The present invention is implemented by an optical measurement instrument for a living body that comprises: a plurality of light irradiation means for irradiating an object being tested with light; a plurality of light reception means for detecting the light that is irradiated to the object being tested and propagates inside the object being tested; and a calculation unit for calculating hemodynamics of a plurality of factors inside the object being tested from the detection signal detected by the light reception means, wherein the calculation unit comprises: a separation calculation unit for separating a plurality of signals from the detection signal; a selection unit for calculating a correlation coefficient between a plurality of the signals separated by the separation calculation unit to select a signal in which the correlation coefficient is a predetermined value or greater as a noise signal; and a reconstruction calculation unit for reconfiguring the detection signal excluding the noise signal selected by the selection unit.
- A plurality of factors can be detected from optical signals measured through the use of a plurality of light sources. For example, lights of 690 nm and 830 nm can be used to detect oxygenated hemoglobin and deoxygenated hemoglobin (Non-Patent Document 2). The separation calculation unit calculates a plurality independent components (separation signals), and mixing coefficient matrices that indicate how each of the independent components contributes to a plurality of measurement signals by means of, for example, an independent component analysis method (Non-Patent Document 3). The independent component analysis method can recover original measurement signals from the independent components and mixture coefficient matrices.
- The selection unit selects separation signals based on a certain selection criterion. For example, when trying to select noise components from the separation signals separated by the independent component analysis method, the selection unit obtains a correlation coefficient between a plurality of measurement factors for each of the mixture coefficient matrices corresponding to each separation signal that is binarized by the positive and negative of the value, and selects a separation signal indicating a correlation coefficient greater than a certain value as a noise component.
- The reconstruction calculation unit uses only the separation signals selected by the selection unit or other separation signals to reconfigures signals. For example, when the separation signals obtained by the independent component analysis are selected as the noise components, the reconstruction calculation unit performs configuration by setting the mixture coefficient matrices of the selected separation components to zero. This allows the noises to be removed.
- The present invention allows the noise components, such as the body motion noises that are difficult to be separated and removed, to be separated and removed, and provides signals with few noise components.
- Other objects, features and advantages of the invention will become apparent from the following description of the embodiments of the invention taken in conjunction with the accompanying drawings.
-
FIG. 1 is a schematic diagram of an embodiment of the present invention; -
FIG. 2 is an exemplary interface unit in the embodiment; -
FIG. 3 is exemplary signals measured by the interface unit ofFIG. 2 ; -
FIG. 4 is exemplary separations signals calculated by a separation calculation unit; -
FIG. 5 is exemplary mixing coefficient matrices calculated by the separation calculation unit; -
FIG. 6 is an exemplary display of the separation signals as a result of selection by a selection unit using the mixing coefficient matrices; -
FIG. 7 is an example in which the mixing coefficient matrices calculated by the separation calculation unit are binarized; and -
FIG. 8 is an exemplary display of the separation signals as a result of selection by the selection unit using the binarized mixing coefficient matrices. - A specific block diagram is shown in
FIG. 1 as an embodiment. An interface unit (110) for optical measurement of a living body is attached to part of or the whole of the head of a subject being tested. A light irradiation unit (101) irradiates the living body with mixed lights of 690 nm and 830 nm wavelengths through an optical fiber (111) coupled to the interface unit (110). The lights of each wavelength, which passed through the living body via the optical fiber (111) coupled to the interface unit (110), are detected by a light detection unit (102), and the detection result is stored in a storage unit (103). Here, light of other wavelengths may also be irradiated, and the number of combination of wavelengths employed may be three or more. Moreover, the obtained optical signals may be employed to the processing described in the following sections, after they are subjected to some sort of calculation processing. - Next, the optical signals are separated by a separation calculation unit (121) based on the method and setting selected by a user. The signals separated by the separation calculation unit (referred to as separation signals) are manually or automatically selected by a selection unit (122) based on the method and setting selected by the user. The signals selected by the selection unit (122) are reconfigured by a reconstruction calculation unit (123) based on the method and setting selected by the user.
- A parameter input and a result display at the separation calculation unit, selection unit, and reconstruction calculation unit are performed by a display unit. Here, the parameter refers to, for example, a setting value for the signal separation used by the separation calculation unit, and a threshold of the correlation coefficient used by the selection unit.
- Now, detailed description will be provided. First, an example of the interface unit is shown in
FIG. 2 as 200. Measure points at this time include a measure point 1 (200 a), a measure point 2 (200 b), a measure point 3 (200 c), and a measure point 4 (200 d), totaling four points. Here, the measure point refers to a middle point between an irradiation point and a detection point. The present embodiment provides an example in which concentration changes in oxygenated hemoglobin and deoxygenated hemoglobin are calculated from optical signals (Non-Patent Document 2). Examples of signals calculated from the obtained optical signals are shown inFIG. 3 .FIG. 3 shows oxygenated hemoglobin concentration changes (201 to 204) measured at four measure points, and deoxygenated hemoglobin concentration changes (301 to 304) measured at the same measure points. A lateral axis indicates time (second), while a vertical axis indicates hemoglobin concentration change amount (mM-mm). These are measurement signals. Description of reference numerals used inFIG. 3 is as follows. - 201 oxygenated hemoglobin concentration changes at a measure point 1
202 oxygenated hemoglobin concentration changes at a measure point 2
203 oxygenated hemoglobin concentration changes at a measure point 3
204 oxygenated hemoglobin concentration changes at a measure point 4
301 deoxygenated hemoglobin concentration changes at a measure point 1
302 deoxygenated hemoglobin concentration changes at a measure point 2
303 deoxygenated hemoglobin concentration changes at a measure point 3
304 deoxygenated hemoglobin concentration changes at a measure point 4 - Next, one of the specific calculation techniques performed by the separation calculation unit (121) will be described. A TDDICA method (Non-Patent Document 3) may be used as one of the separation techniques. This is one of the techniques referred to as the independent component analysis. A measurement signal X comprised of time series data simultaneously obtained at respective four measure points during the measurement is expressed as follows:
-
- wherein T indicates the number of data.
- A set of signals comprised of four statistically independent time series signal components (original signal S) is expressed as follows:
-
- Assuming that the measurement signal X is expressed in a linear mixing signal that is expressed as follows:
-
X=AS [Expression 3] - in which the set of signals S is multiplied by a mixing coefficient matrix A (indicating to what extent respective original signals S are included in the respective measurement signal X),
- the mixing coefficient matrix A and original signal S are estimated from only the independent condition of the original signal component and the measurement signal x. Here, the number of samplings refers to the number of data in the time direction. The estimation is performed by whitening and rotation. The whitening means conversion into a linear independent signal by mixture of the measurement signals. It is given using the covariance matrix C of the measurement signal X as follows:
-
X 0 =C −1/2 X [Expression 4] - The estimate U of the original signal s comprised of independent components is given as follows by using a proper rotation matrix R:
-
U=RX0 [Expression 5] - This rotation matrix is determined by the above independence condition. Here, that a delay correlation becomes zero at K delay times·k (k=1, 2, . . . , k) is determined as independence. Actually, it is difficult to achieve this. Therefore, minimizing the size of correlation is set as a condition instead. This is achieved by the rotation matrix R that minimizes a loss function expressed by the following expression:
-
- where x′ indicates a dislocation, while <•>t indicates an average time. The estimate of the mixing coefficient matrix is given by the following expression:
-
W=C 1/2 R t [Expression 7] - The noise applied to each sensor may b handled as the positive. In such as case, in stead of the [expression 3], following expression is used in which an observation noise is added to the right side thereof.
-
X=AS+N [Expression 8] - At this time, pseudo whitening is performed instead of whitening, and the pseudo whitening is performed in the following manner using covariance matrix G.
-
A 0=(C−G)−1/2 X [Expression 9] - Subsequent processing is the same as the case where sensor noise is not considered. The matrix G must be estimated separately. It is possible to estimate the matrix G from a residual when one section in a rest state is fitted by a polynomial expression. A description has been made as to a case in which the number of measurement points is four. However, the separation calculation can be conducted by the same procedure even when the number of measurement points differs from that.
- U: The separation signal and mixing coefficient matrix R are obtained by calculations above mentioned.
- An example of separation signals is shown in
FIG. 4. 211 to 214 are separation signals calculated from the oxygenated hemoglobin concentration changes obtained at four measure points. 311 to 314 are separation signals calculated from deoxygenated hemoglobin concentration changes obtained at four measure points. Description of reference numerals used inFIG. 4 is as follows. - 211 separation signal 1 of oxygenated hemoglobin concentration changes
212 separation signal 2 of oxygenated hemoglobin concentration changes
213 separation signal 3 of oxygenated hemoglobin concentration changes
214 separation signal 4 of oxygenated hemoglobin concentration changes
311 separation signal 1 of deoxygenated hemoglobin concentration changes
312 separation signal 2 of deoxygenated hemoglobin concentration changes
313 separation signal 3 of deoxygenated hemoglobin concentration changes
314 separation signal 4 of deoxygenated hemoglobin concentration changes - An example of mixing coefficient matrices is shown in
FIG. 5. 231 to 234 indicate the mixingcoefficient matrices 211 to 214, respectively. Color shading indicates the size of the coefficient. Description of reference numerals used inFIG. 5 is as follows. - 231 mixing coefficient matrix of the separation signal 1 of oxygenated hemoglobin concentration changes
232 mixing coefficient matrix of the separation signal 2 of oxygenated hemoglobin concentration changes
233 mixing coefficient matrix of the separation signal 3 of oxygenated hemoglobin concentration changes
234 mixing coefficient matrix of the separation signal 4 of oxygenated hemoglobin concentration changes
234 a mixing coefficient matrix at a measure point 1 of the separation signal 1 of oxygenated hemoglobin concentration changes
234 b mixing coefficient matrix at a measure point 2 of the separation signal 1 of oxygenated hemoglobin concentration changes
234 c mixing coefficient matrix at a measure point 3 of the separation signal 1 of oxygenated hemoglobin concentration changes
234 d mixing coefficient matrix at a measure point 4 of the separation signal 1 of oxygenated hemoglobin concentration changes
331 mixing coefficient matrix of the separation signal 1 of deoxygenated hemoglobin concentration changes
332 mixing coefficient matrix of the separation signal 2 of deoxygenated hemoglobin concentration changes
333 mixing coefficient matrix of the separation signal 3 of deoxygenated hemoglobin concentration changes
334 mixing coefficient matrix of the separation signal 4 of deoxygenated hemoglobin concentration changes - The mixing
coefficient matrix 234 of theseparation signal 214 will be described in more detail. 234 a to 234 d corresponds to the measure points 1 to 4, respectively. More specifically, it is indicated that theseparation signal 214 is strongly mixed at the measure point 2 (234 b), and hardly mixed at the measure point 3 (234 c). - Here, for highly precise signal separation, filtering processing may be performed before the separation calculation. When the frequency band for a target signal is known, a filtering processing is effective that reduces frequency band differing from the frequency band.
- Next, the selection of the separations signals by the selection unit (122) will be described. Here, a case is considered in which a noise related to the motion of a body (not brain activity) is selected. The noise is referred to as a body motion noise hereinafter. In other words, it is a noise that occurs when a body moves. While the occurrence mechanism of the body motion noise is not known sufficiently yet, it is thought that the body motion noise occurs due to a change in the contact state between a measurement interface and a subject to be measured caused by the body motion, a dynamic change in the internal state of the subject to be measured, and the like. When the signal of brain activity evoked by body motion is a target signal, it is extremely difficult to separate the target signal from a body motion signal, because the body motion signal occurs in synchronization with the brain activity as the target signal. However, it is absolutely necessary to discriminate them. For the change with time of the body motion noise, similarity is not guaranteed between different measurement factors (here, oxygenated hemoglobin and deoxygenated hemoglobin). Therefore, it is difficult to select the body motion noise from the correlation of the calculated separation signals.
- However, it is thought that the extent of the influence the body motion noise has on each measure point is equal between different measurement factors. Therefore, it is possible to select the body motion noise from the correlation of the calculated mixing coefficient matrices.
- Accordingly, the correlation between the oxygenated hemoglobin and deoxygenated hemoglobin is calculated from the mixing coefficient matrices shown in
FIG. 5 , and an example is shown inFIG. 6 in which a correlation coefficient with correlation higher than a certain value (here 0.5) is represented in black, while that with correlation lower than the certain value is represented in white. Description of reference numerals used inFIG. 6 is as follows. - 421 column corresponding to the mixing coefficient matrix of the separation signal 1 of oxygenated hemoglobin concentration changes
422 column corresponding to the mixing coefficient matrix of the separation signal 2 of oxygenated hemoglobin concentration changes
423 column corresponding to the mixing coefficient matrix of the separation signal 3 of oxygenated hemoglobin concentration changes
424 column corresponding to the mixing coefficient matrix of the separation signal 4 of oxygenated hemoglobin concentration changes
431 row corresponding to the mixing coefficient matrix of the separation signal 1 of deoxygenated hemoglobin changes
432 row corresponding to the mixing coefficient matrix of the separation signal 2 of deoxygenated hemoglobin changes
433 row corresponding to the mixing coefficient matrix of the separation signal 3 of deoxygenated hemoglobin changes
434 row corresponding to the mixing coefficient matrix of the separation signal 4 of deoxygenated hemoglobin changes
441 correlation coefficient between the mixing coefficient matrix of the separation signal 1 of oxygenated hemoglobin concentration changes and the mixing coefficient matrix of the separation signal 1 of deoxygenated hemoglobin concentration changes that is higher than 0.5 and is represented in black.
442 correlation between the mixing coefficient matrix of the separation signal 2 of oxygenated hemoglobin concentration changes and the mixing coefficient matrix of the separation signal 1 of deoxygenated hemoglobin concentration changes that is higher than 0.5 and is represented in black. - In
FIG. 6 ,column 421 to 424 correspond to the mixingcoefficient matrices 231 to 234, respectively, androws 431 to 434 correspond to mixingcoefficient matrices 331 to 334, respectively. In other words, 441 indicates that the correlation coefficient between the mixingcoefficient matrix 231 and mixingcoefficient matrix 331 is greater than 0.5, while 442 indicates that the correlation between the mixingcoefficient matrix 232 and mixingcoefficient matrix 332 is greater than 0.5. - Therefore, the result indicates that the
separation components separation components - Incidentally, it is known that the influence that changes due to the brain activity have on each measure point tends to be similar between the oxygenated hemoglobin and deoxygenated hemoglobin. Therefore, it is likely that the brain activity is also included in the result shown in
FIG. 6 in addition to the body motion noises. - Therefore, attention is directed to the difference between the brain activity and body motion noise. In other words, when a certain separation signal is a brain activity, and it exists in a limited area, the mixing coefficient matrix has a high value at a particular measure point (measure point 3 for 233 and 333 in
FIG. 5 ) and has a low value at other measure points. Similarity of each value between oxygenated hemoglobin and deoxygenated hemoglobin is thought to be low at the measure point where the value is low because of the complexity of the living body. - However, since the body motion noise is physically caused, the similarity of the mixing coefficient matrix between the oxygenated hemoglobin and deoxygenated hemoglobin is thought to be high even at the measure point where the value is low.
- Therefore, in order to extract the characteristic of the point where the value of the mixing coefficient matrix is low, binarization using the positive and negative of the value is performed. More specifically, binarization is performed using zero as a threshold value. An example of binarized mixing coefficient matrices is shown in
FIG. 7. 241 to 244 and 341 to 344show 231 to 234 and 331 to 334 ofFIG. 5 that are binarized, respectively. - Description of reference numerals used in
FIG. 7 is as follows.
241 binarized mixing coefficient matrix of separation signal 1 of oxygenated hemoglobin concentration changes
242 binarized mixing coefficient matrix of separation signal 2 of oxygenated hemoglobin concentration changes
243 binarized mixing coefficient matrix of separation signal 3 of oxygenated hemoglobin concentration changes
244 binarized mixing coefficient matrix of separation signal 4 of oxygenated hemoglobin concentration changes
341 binarized mixing coefficient matrix of separation signal 1 of deoxygenated hemoglobin concentration changes
342 binarized mixing coefficient matrix of separation signal 2 of deoxygenated hemoglobin concentration changes
343 binarized mixing coefficient matrix of separation signal 3 of deoxygenated hemoglobin concentration changes
344 binarized mixing coefficient matrix of separation signal 4 of deoxygenated hemoglobin concentration changes - Correlation coefficients are calculated from these binarized mixing coefficient matrices, and an example of display is shown in
FIG. 8 in which correlation coefficients that are greater than a certain threshold (0.5 here) are selected. - Description of reference numerals used in
FIG. 8 is as follows. - 521 column corresponding to the binarized mixing coefficient matrix of the separation signal 1 of oxygenated hemoglobin concentration changes
522 column corresponding to the binarized mixing coefficient matrix of the separation signal 2 of oxygenated hemoglobin concentration changes
523 column corresponding to the binarized mixing coefficient matrix of the separation signal 3 of oxygenated hemoglobin concentration changes
524 column corresponding to the binarized mixing coefficient matrix of the separation signal 4 of oxygenated hemoglobin concentration changes
531 row corresponding to the binarized mixing coefficient matrix of the separation signal 1 of deoxygenated hemoglobin concentration changes
532 row corresponding to the binarized mixing coefficient matrix of the separation signal 2 of deoxygenated hemoglobin concentration changes
533 row corresponding to the binarized mixing coefficient matrix of the separation signal 3 of deoxygenated hemoglobin concentration changes
534 row corresponding to the binarized mixing coefficient matrix of the separation signal 4 of deoxygenated hemoglobin concentration changes
541 correlation coefficient between the binarized mixing coefficient matrix of the separation signal 1 of oxygenated hemoglobin concentration changes and binarized mixing coefficient matrix of the separation signal 1 of deoxygenated hemoglobin concentration changes that is greater than 0.5 and represented in black
542 correlation coefficient between the binarized mixing coefficient matrix of the separation signal 2 of oxygenated hemoglobin concentration changes and binarized mixing coefficient matrix of the separation signal 2 of deoxygenated hemoglobin concentration changes that is smaller than 0.5 and is represented in white - In
FIG. 8 , in the same way as inFIG. 6 ,columns 521 to 524 correspond to binarized mixingcoefficient matrices 241 to 244, respectively, androws 531 to 534 correspond to mixingcoefficient matrices 341 to 344, respectively. In other words, 541 indicates that the correlation coefficient between the binarizedmixing coefficient matrix 241 and mixingcoefficient matrix 341 is greater than a certain value, and 542 indicates that the correlation coefficient between the mixingcoefficient matrix 242 and mixingcoefficient matrix 342 is smaller than a certain value. - In
FIG. 6 , 442 is selected as a body motion noise. However, inFIG. 8 , 542 is not selected as the body motion noise. Accordingly, separation signals 212 and 312 are thought to have corresponded to the brain activity. - This means that the present technique allows the selection of the body motion noise as well as the brain activity noise at the same time.
- In the selection of the separation signal, the user may select a body motion noise based on this result or may automatically select it. Although 0.5 is used here as a threshold of the correlation coefficient, other values may be used depending on the measurement signals.
- In the present embodiment, selection is made from among all of the separation signals. However, in order to pay attention to a change in an assignment given to the subject being examined during the measurement, the selection unit may extract separation components having high correlation with the time for the assignment in advance to make selection from among the previously extracted separation components.
- Furthermore, as an indicator of the similarity to the time for the assignment at this time, synchronism of signal change with the assignment that is performed repeatedly may be used. Specifically, when the assignment is repeated, for example, for five times, measurement signals are extracted at times in synchronization therewith, and a correlation coefficient of the extracted five measurement signal fragments is used as an indicator of the synchronism.
- Moreover, while the similarity in binarized mixing coefficient matrix between oxygenated hemoglobin and deoxygenated hemoglobin is used as a selection criterion in the present embodiment, the difference in phase between oxygenated hemoglobin changes and deoxygenated hemoglobin changes during brain activity may be used as the selection criterion. Specifically, an instantaneous phase is calculated through Hilbert transform or the like, and when the phase difference between oxygenated hemoglobin and deoxygenated hemoglobin is smaller than a certain value, it may be determined that the body motion noise exists, while when the phase difference is in the range of a certain value different from the foregoing value, it may be determined that brain activity exists.
- Furthermore, in order to derive the selection criterion, the separation calculation and signal selection are performed by use of the whole measurement signals in the present embodiment. However, the separation calculation and signal selection may also be performed by dividing the measurement signals into parts.
- At this time, the signal selection may be performed using the result of the separation calculation that is performed for part of a plurality of the measurement signals. At this time, the obtained result of a plurality of separation calculations may be grouped by a method, such as clustering, and subsequently each group may be subjected to the selection processing.
- Next, the reconstruction calculation unit (123) performs signal reconstruction based on the signal selection performed by the selection unit. Specifically, when the selection is for the purpose of removing noises, zero is substituted into the mixing coefficient matrix of the separation signal, and reconstruction is performed by a reconstruction method using the independent component analysis method. When the selection is for the purpose of extracting a target signal, zero is substituted into the mixing coefficient matrix of an unselected separation signal, and the reconstruction is performed by the reconstruction method using the independent reconstruction analysis method. It should be noted that although 0 is substituted into the mixing coefficient matrix of the separation signal for the purpose of removing noises in the above, a value that is greater than 0 and smaller than 1 may be multiplied.
- The present invention can be employed for the optical measurement instrument for a living body for measuring the living body.
- It should be further understood by those skilled in the art that although the foregoing description has been made on embodiments of the invention, the invention is not limited thereto and various changes and modifications may be made without departing from the spirit of the invention and the scope of the appended claims.
Claims (9)
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2006-337737 | 2006-12-15 | ||
JP2006337737A JP5034477B2 (en) | 2006-12-15 | 2006-12-15 | Biological light measurement device |
Publications (2)
Publication Number | Publication Date |
---|---|
US20080146901A1 true US20080146901A1 (en) | 2008-06-19 |
US8311600B2 US8311600B2 (en) | 2012-11-13 |
Family
ID=39528331
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US11/954,627 Active 2031-07-11 US8311600B2 (en) | 2006-12-15 | 2007-12-12 | Optical measurement instrument for living body |
Country Status (2)
Country | Link |
---|---|
US (1) | US8311600B2 (en) |
JP (1) | JP5034477B2 (en) |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090324034A1 (en) * | 2008-06-30 | 2009-12-31 | Nellcor Puritan Bennett Ireland | Systems and methods for ridge selection in scalograms of signals |
US20090326871A1 (en) * | 2008-06-30 | 2009-12-31 | Nellcor Puritan Bennett Ireland | Systems and methods for artifact detection in signals |
US20100014761A1 (en) * | 2008-06-30 | 2010-01-21 | Nellcor Puritan Bennett Llc | Methods And Systems For Discriminating Bands In Scalograms |
USD626562S1 (en) | 2008-06-30 | 2010-11-02 | Nellcor Puritan Bennett Llc | Triangular saturation pattern detection indicator for a patient monitor display panel |
USD626561S1 (en) | 2008-06-30 | 2010-11-02 | Nellcor Puritan Bennett Llc | Circular satseconds indicator and triangular saturation pattern detection indicator for a patient monitor display panel |
US8216136B2 (en) | 2009-03-05 | 2012-07-10 | Nellcor Puritan Bennett Llc | Systems and methods for monitoring heart rate and blood pressure correlation |
US8290730B2 (en) | 2009-06-30 | 2012-10-16 | Nellcor Puritan Bennett Ireland | Systems and methods for assessing measurements in physiological monitoring devices |
JP2013192905A (en) * | 2012-03-23 | 2013-09-30 | Hitachi Ltd | Biophotonic measuring apparatus |
US20150045640A1 (en) * | 2012-02-20 | 2015-02-12 | Hamamatsu Photonics K.K. | Concentration measurement device and concentration measurement method |
USD778928S1 (en) * | 2013-12-02 | 2017-02-14 | Symantec Corporation | Display screen with graphical user interface |
US10194870B2 (en) | 2015-05-27 | 2019-02-05 | Covidien Lp | Systems and methods for optimizing autoregulation measurements |
US10219705B2 (en) | 2015-05-08 | 2019-03-05 | Covidien Lp | System and method for identifying autoregulation zones |
US10226188B2 (en) | 2013-08-23 | 2019-03-12 | Covidien Lp | Systems and methods for monitoring blood pressure |
JP2019037713A (en) * | 2017-08-29 | 2019-03-14 | 日本光電工業株式会社 | Pulse oximeter |
US10271779B2 (en) | 2015-06-30 | 2019-04-30 | Covidien Lp | System and method of monitoring autoregulation |
US10292663B2 (en) | 2015-06-30 | 2019-05-21 | Covidien Lp | System and method of monitoring autoregulation |
US10383579B2 (en) | 2014-10-16 | 2019-08-20 | Covidien Lp | System and method for monitoring autoregulation |
US10463292B2 (en) | 2015-10-16 | 2019-11-05 | Covidien Lp | System and method for identifying autoregulation zones |
US10499818B2 (en) | 2015-10-19 | 2019-12-10 | Covidien Lp | System and method for providing blood pressure safe zone indication during autoregulation monitoring |
US10736578B2 (en) | 2016-07-14 | 2020-08-11 | Covidien Lp | Systems and methods of monitoring autoregulation |
US10932724B2 (en) | 2015-06-17 | 2021-03-02 | Covidien Lp | Systems and methods for monitoring autoregulation using a confidence level |
US11096588B2 (en) | 2015-10-06 | 2021-08-24 | Covidien Lp | System and method for monitoring autoregulation utilizing normalized regional oxygen saturation values |
US11419506B2 (en) | 2016-08-22 | 2022-08-23 | Covidien Lp | System and method for identifying blood pressure zones during autoregulation monitoring |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5466389B2 (en) * | 2008-10-30 | 2014-04-09 | 株式会社日立メディコ | Biological light measurement apparatus and measurement noise estimation method |
CN102333489A (en) * | 2009-03-19 | 2012-01-25 | 株式会社岛津制作所 | Photobiological measuring device and analyzing method |
JP5909388B2 (en) * | 2012-03-09 | 2016-04-26 | 株式会社日立製作所 | Biological light measurement device |
JP6028100B2 (en) * | 2013-07-11 | 2016-11-16 | 株式会社日立製作所 | Biological light measurement data analysis apparatus, analysis method, and program therefor |
JP2017138254A (en) | 2016-02-05 | 2017-08-10 | 国立研究開発法人海洋研究開発機構 | System and method for estimating resources |
JP6412074B2 (en) * | 2016-09-07 | 2018-10-24 | 浜松ホトニクス株式会社 | Concentration measuring device and concentration measuring method |
EP3731755A4 (en) | 2017-12-29 | 2021-08-04 | Raydiant Oximetry, Inc. | Systems, devices, and methods for performing trans-abdominal fetal oximetry and/or trans-abdominal fetal pulse oximetry using independent component analysis |
JP6675693B2 (en) * | 2019-03-22 | 2020-04-01 | 国立研究開発法人海洋研究開発機構 | Resource estimation system and resource estimation method |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4753246A (en) * | 1986-03-28 | 1988-06-28 | The Regents Of The University Of California | EEG spatial filter and method |
US20040039267A1 (en) * | 2000-10-16 | 2004-02-26 | Shingo Kawasaki | Organism optical measurement instrument |
US6701170B2 (en) * | 2001-11-02 | 2004-03-02 | Nellcor Puritan Bennett Incorporated | Blind source separation of pulse oximetry signals |
US20040106856A1 (en) * | 2002-11-18 | 2004-06-03 | Masahiro Kimura | Optical measuring apparatus and method |
US20040171919A1 (en) * | 2003-02-28 | 2004-09-02 | Hitachi, Ltd. | Living body light measuring device |
WO2005046483A1 (en) * | 2003-11-12 | 2005-05-26 | Hitachi Medical Corporation | Biological light measuring apparatus |
US20060006343A1 (en) * | 2004-07-08 | 2006-01-12 | Naoki Tanaka | Optical measurement system for living bodies and method thereof |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3588880B2 (en) | 1995-11-17 | 2004-11-17 | 株式会社日立製作所 | Biological light measurement device |
JP3863095B2 (en) | 2002-11-25 | 2006-12-27 | 本田技研工業株式会社 | Optical measuring device |
JP4631510B2 (en) * | 2005-03-31 | 2011-02-16 | 株式会社島津製作所 | Brain function information monitoring device |
-
2006
- 2006-12-15 JP JP2006337737A patent/JP5034477B2/en not_active Expired - Fee Related
-
2007
- 2007-12-12 US US11/954,627 patent/US8311600B2/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4753246A (en) * | 1986-03-28 | 1988-06-28 | The Regents Of The University Of California | EEG spatial filter and method |
US20040039267A1 (en) * | 2000-10-16 | 2004-02-26 | Shingo Kawasaki | Organism optical measurement instrument |
US6701170B2 (en) * | 2001-11-02 | 2004-03-02 | Nellcor Puritan Bennett Incorporated | Blind source separation of pulse oximetry signals |
US20040106856A1 (en) * | 2002-11-18 | 2004-06-03 | Masahiro Kimura | Optical measuring apparatus and method |
US20040171919A1 (en) * | 2003-02-28 | 2004-09-02 | Hitachi, Ltd. | Living body light measuring device |
WO2005046483A1 (en) * | 2003-11-12 | 2005-05-26 | Hitachi Medical Corporation | Biological light measuring apparatus |
US20070142719A1 (en) * | 2003-11-12 | 2007-06-21 | Shingo Kawasaki | Optical measurement apparatus for living body |
US20060006343A1 (en) * | 2004-07-08 | 2006-01-12 | Naoki Tanaka | Optical measurement system for living bodies and method thereof |
Cited By (38)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8295567B2 (en) | 2008-06-30 | 2012-10-23 | Nellcor Puritan Bennett Ireland | Systems and methods for ridge selection in scalograms of signals |
US8289501B2 (en) | 2008-06-30 | 2012-10-16 | Nellcor Puritan Bennett Ireland | Methods and systems for discriminating bands in scalograms |
US20100014761A1 (en) * | 2008-06-30 | 2010-01-21 | Nellcor Puritan Bennett Llc | Methods And Systems For Discriminating Bands In Scalograms |
USD626562S1 (en) | 2008-06-30 | 2010-11-02 | Nellcor Puritan Bennett Llc | Triangular saturation pattern detection indicator for a patient monitor display panel |
USD626561S1 (en) | 2008-06-30 | 2010-11-02 | Nellcor Puritan Bennett Llc | Circular satseconds indicator and triangular saturation pattern detection indicator for a patient monitor display panel |
US8077297B2 (en) | 2008-06-30 | 2011-12-13 | Nellcor Puritan Bennett Ireland | Methods and systems for discriminating bands in scalograms |
US20090326871A1 (en) * | 2008-06-30 | 2009-12-31 | Nellcor Puritan Bennett Ireland | Systems and methods for artifact detection in signals |
US9113815B2 (en) | 2008-06-30 | 2015-08-25 | Nellcor Puritan Bennett Ireland | Systems and methods for ridge selection in scalograms of signals |
US8827917B2 (en) | 2008-06-30 | 2014-09-09 | Nelleor Puritan Bennett Ireland | Systems and methods for artifact detection in signals |
USD736250S1 (en) | 2008-06-30 | 2015-08-11 | Covidien Lp | Portion of a display panel with an indicator icon |
US8483459B2 (en) | 2008-06-30 | 2013-07-09 | Nèllcor Puritan Bennett Ireland | Systems and methods for ridge selection in scalograms of signals |
US20090324034A1 (en) * | 2008-06-30 | 2009-12-31 | Nellcor Puritan Bennett Ireland | Systems and methods for ridge selection in scalograms of signals |
US8216136B2 (en) | 2009-03-05 | 2012-07-10 | Nellcor Puritan Bennett Llc | Systems and methods for monitoring heart rate and blood pressure correlation |
US8932219B2 (en) | 2009-03-05 | 2015-01-13 | Nellcor Puritan Bennett Ireland | Systems and methods for monitoring heart rate and blood pressure correlation |
US8290730B2 (en) | 2009-06-30 | 2012-10-16 | Nellcor Puritan Bennett Ireland | Systems and methods for assessing measurements in physiological monitoring devices |
US10405783B2 (en) | 2012-02-20 | 2019-09-10 | Hamamatsu Photonics K.K. | Concentration measurement device and concentration measurement method |
EP3092950A1 (en) * | 2012-02-20 | 2016-11-16 | Hamamatsu Photonics K.K. | Concentration measurement device and concentration measurement method |
US9700248B2 (en) * | 2012-02-20 | 2017-07-11 | Hamamatsu Photonics K.K. | Concentration measurement device and concentration measurement method |
US20150045640A1 (en) * | 2012-02-20 | 2015-02-12 | Hamamatsu Photonics K.K. | Concentration measurement device and concentration measurement method |
JP2013192905A (en) * | 2012-03-23 | 2013-09-30 | Hitachi Ltd | Biophotonic measuring apparatus |
US11160464B2 (en) | 2013-08-23 | 2021-11-02 | Covidien Lp | Systems and methods for monitoring blood pressure |
US10226188B2 (en) | 2013-08-23 | 2019-03-12 | Covidien Lp | Systems and methods for monitoring blood pressure |
USD778928S1 (en) * | 2013-12-02 | 2017-02-14 | Symantec Corporation | Display screen with graphical user interface |
US10383579B2 (en) | 2014-10-16 | 2019-08-20 | Covidien Lp | System and method for monitoring autoregulation |
US10219705B2 (en) | 2015-05-08 | 2019-03-05 | Covidien Lp | System and method for identifying autoregulation zones |
US10194870B2 (en) | 2015-05-27 | 2019-02-05 | Covidien Lp | Systems and methods for optimizing autoregulation measurements |
US10932724B2 (en) | 2015-06-17 | 2021-03-02 | Covidien Lp | Systems and methods for monitoring autoregulation using a confidence level |
US10271779B2 (en) | 2015-06-30 | 2019-04-30 | Covidien Lp | System and method of monitoring autoregulation |
US10292663B2 (en) | 2015-06-30 | 2019-05-21 | Covidien Lp | System and method of monitoring autoregulation |
US11096588B2 (en) | 2015-10-06 | 2021-08-24 | Covidien Lp | System and method for monitoring autoregulation utilizing normalized regional oxygen saturation values |
US10463292B2 (en) | 2015-10-16 | 2019-11-05 | Covidien Lp | System and method for identifying autoregulation zones |
US10499818B2 (en) | 2015-10-19 | 2019-12-10 | Covidien Lp | System and method for providing blood pressure safe zone indication during autoregulation monitoring |
US11653840B2 (en) | 2015-10-19 | 2023-05-23 | Covidien Lp | System and method for providing blood pressure safe zone indication during autoregulation monitoring |
US10736578B2 (en) | 2016-07-14 | 2020-08-11 | Covidien Lp | Systems and methods of monitoring autoregulation |
US11419506B2 (en) | 2016-08-22 | 2022-08-23 | Covidien Lp | System and method for identifying blood pressure zones during autoregulation monitoring |
JP2019037713A (en) * | 2017-08-29 | 2019-03-14 | 日本光電工業株式会社 | Pulse oximeter |
JP7020820B2 (en) | 2017-08-29 | 2022-02-16 | 日本光電工業株式会社 | Pulse oximeter |
US11331046B2 (en) | 2017-08-29 | 2022-05-17 | Nihon Kohden Corporation | Pulse oximeter |
Also Published As
Publication number | Publication date |
---|---|
US8311600B2 (en) | 2012-11-13 |
JP5034477B2 (en) | 2012-09-26 |
JP2008148797A (en) | 2008-07-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8311600B2 (en) | Optical measurement instrument for living body | |
Liu et al. | Exploring the time–frequency content of high frequency oscillations for automated identification of seizure onset zone in epilepsy | |
CN111340142B (en) | Epilepsia magnetoencephalogram spike automatic detection method and tracing positioning system | |
EP1685801B1 (en) | Biological light measuring apparatus | |
US9131889B2 (en) | Cognitive impairment determination apparatus, cognitive impairment determination system and program | |
EP3281581B1 (en) | Annotation of a wavefront | |
JP5948831B2 (en) | Measuring device, measuring method, program, and recording medium | |
US10380737B2 (en) | Detection of microaneurysms | |
US20150025340A1 (en) | Calibration curve creating method and apparatus for the same, and blood component calibration apparatus | |
US7693689B2 (en) | Noise-component removing method | |
KR102482796B1 (en) | Device and method for non-contact hrv analysis using remote ppg | |
Pu et al. | Novel tailoring algorithm for abrupt motion artifact removal in photoplethysmogram signals | |
KR101601041B1 (en) | Pattern Classification Apparatus and Method for fMRI | |
Karthika et al. | Retinal image analysis using contourlet transform and multistructure elements morphology by reconstruction | |
CN107510462B (en) | Blood oxygen measuring method and device | |
Vintila et al. | Pupil response as an indicator of hazard perception during simulator driving | |
CN110602978A (en) | System and method for extracting physiological information from video sequence | |
CN109620198B (en) | Cardiovascular index detection and model training method and device | |
Besheli et al. | A sparse representation strategy to eliminate pseudo-HFO events from intracranial EEG for seizure onset zone localization | |
CN110547768B (en) | Near-infrared brain function imaging quality control method and control system | |
KR102509809B1 (en) | Method for noise reduction of functional near-infrared spectroscopy | |
O'Regan et al. | Automatic detection of EEG artefacts arising from head movements using gyroscopes | |
KR20160100140A (en) | The method and apparatus for measuring bio signal | |
KR102574085B1 (en) | Apparatus and method for predicting blood glucose concentration | |
US20110279775A1 (en) | Apparatus and Method for Diagnosing Disease Involving Optic Nerve |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: HITACHI, LTD., JAPAN Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:KATURA, TAKUSIGE;TANAKA, NAOKI;MAKI, ATSUSHI;REEL/FRAME:020233/0472;SIGNING DATES FROM 20071115 TO 20071119 Owner name: HITACHI, LTD., JAPAN Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:KATURA, TAKUSIGE;TANAKA, NAOKI;MAKI, ATSUSHI;SIGNING DATES FROM 20071115 TO 20071119;REEL/FRAME:020233/0472 |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
FEPP | Fee payment procedure |
Free format text: PAYOR NUMBER ASSIGNED (ORIGINAL EVENT CODE: ASPN); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |
|
FPAY | Fee payment |
Year of fee payment: 4 |
|
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 8TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: M1552); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY Year of fee payment: 8 |