JP5916187B2 - Brain habit diagnosis apparatus and program - Google Patents

Description

  The present invention relates to a brain habit diagnosis apparatus and program for instructing to diagnose or improve a brain habit which is a state of brain activity associated with a habit of a living body.

  In general, it is known that a living body has various habits (physical behavior, behavioral patterns such as mental and psychological tendency and way of thinking).

  In the past, the habits of each living body were recognized only when they appeared on the table in the form of actions for a certain problem.

  On the other hand, brain activity has been detected by the reaction induced by the load of the task for the activity of each part of the brain. Known detection methods include brain waves, fMRI (Magnetic Resonance Imaging), PET (Positron CT), NIRS (Near-infrared Spectroscopy), and the like.

  For NIRS, non-invasive monitoring of brain function is performed using changes in the concentration of local oxidized hemoglobin (OxyHb), deoxidized hemoglobin (DeoxyHb), total hemoglobin (TotalHb), etc. NIRS is less restrictive than PET and fMRI and can be measured while moving. Therefore, there is an advantage that the cerebral circulation metabolism associated with the exercise task for moving the trunk and the walking task can be analyzed.

  The study of functional brain imaging using the localization of the index variation obtained from NIRS was first conducted in 1991 by a research group at the National Center for Neurology and Psychiatry in Japan. The result is Takashima S, Kato T, Hirano S, Mito T (1992) Observation of activation in local brain blood flow by means of near-infrared spectroscopy in Comprehensive Research Report Concerning Medical Care for Children (People) with Disabilities. Health and Welfare 179-181. And Kato T, Kamei A, Takashima S, Ozaki T (1993) Human visual cortical function during photic stimulation monitoring by means of near-infrared spectroscopy. By J Cereb Blood Flow Metab 13: 516-520. Reported.

  In this study, the NIRS device was used to demonstrate local activation for the first time using the localization of human visual cortical function during and after light stimulation. The first activation pattern of Hb (hemoglobin) demonstrated in the human brain in this study was an increase in OxyHb and TotalHb in the visual cortex and a slight increase in DeoxyHb during light stimulation.

  On the other hand, in the prefrontal cortex, the three Hb indices did not change.

  From this demonstration, the scalp on the cortex to be measured is sandwiched between a pair of probes around 3.0 cm to extract the local brain activity that causes a stronger change during the task than the effect of the scalp signal at rest. The basic technology for obtaining functional NIRS images was established.

  This study is based on the inventor's bold idea to identify local activity in the cortex as a functional pixel by distinguishing local brain activity during stimulation from resting (Kato T (2004) Principle and technique of NIRS imaging for human brain FORCE: fast-oxygen response in capillary event. See Proceedings of ISBET 1270: 88-99).

  This idea is based on optical CT (computed tomography) (Jobsis FF (1977) Method and apparatus for monitoring metabolism in body organs. United States Patent) 4281645.), DOT (diffusion optical tomography) (Barbour RL, Lubowsky J, Aronson R (1989) Method of imaging a random medium. US Patent no. 5137355. http://www.patents.com/us-5137355.html It was a functional NIRS-imaging method that was different from the measurement principle. By measuring temporal changes due to brain activity during the task, contamination of the scalp blood and movement artifacts included in the measurement at rest are significantly attenuated, so the biggest problem of optical CT technology is functional NIRS -imaging was solving.

  Detecting these brain activities and displaying their distribution has been widely studied and has become one of the techniques for mapping brain functions.

  The present inventor made an invention using a two-dimensional diagram as an index for evaluating brain function using the ratio of changes in OxyHb and DeoxyHb as a quantitative value, and obtained a patent (see Patent Document 1).

The inventor has invented an image method for enhancing white matter based on a contrast image of a brain photographed by an MRI apparatus, and filed a patent application (see Patent Document 2). By the invention using the white matter-enhanced image, it is possible to show the result of the growth process of the individual's brain so far.
Japanese Patent No. 4031438 International Publication No. 2009/005013

In the above conventional example, there are a problem related to the view of the brain and a problem related to brain measurement.
(The problem of how to view the brain)
(1) The growth process of the brain so far has a different shape for each individual as shown in the white matter enhanced image shown in Patent Document 2 filed by the present inventor. However, it was not the subject of brain measurements to see what kind of brain activity would differ as a result of different shapes.

  (2) In other words, brain function mapping, which has been widely performed, simply loads the task and extracts the response in the field, so what experience the subject has accumulated so far and what brain state Qualitatively or quantitatively was not judged at all to determine whether or not

  (3) Individuals, how the brain of the individual has accumulated so far, and what kind of brain state each part (for each brain address) is considered as "brain habits" Because there was no concept or word of “brain habit”, diagnosing brain habits could not be done qualitatively or quantitatively.

  (4) Since brain habit is the state of the subject's brain before the measurement, the measurement technology that captures the brain function reaction performed on the spot reflects the state of the subject's brain before the measurement. I never thought about it.

  (5) Since the brain habit was not judged, there was no objective indicator of whether the brain would actually change if the habit was changed.

  (6) It was impossible to give instructions on how to improve habits, and there was no instruction device for changing habits quantitatively and objectively based on brain habit data.

  (7) Since there was no device for determining the degree of brain habits by theme (breathing habits, listening habits, watching habits, speaking habits, thinking habits, etc.), whether or not the habits were acquired and what habits I could not objectively judge whether I was strong or not, bad brain habits, and good brain habits.

  (8) For example, regarding brain education, there was already no way to determine whether you have learned or not. I couldn't judge how many words I knew and how often I used them, and the habit of using them, one word at a time.

  (9) For example, regarding the exercise habits of the brain, there has been no way to determine whether there is a habit of having heavy luggage or a habit of not having heavier objects than a ballpoint pen.

(10) For example, in the field of dentistry, the possibility that the way of breathing affects the brain has been pointed out. Oral breathing habits that breathe through the mouth, which is usually a habit of breathing through the nose, are said to cause sleep problems, reduced concentration, and increased oral infections. However, it was not possible to distinguish from the brain the difference in the habits of mouth breathing and nasal breathing.
(Challenges related to brain measurement)
(1) Positron tomography (hereinafter referred to as PET), which has been used for conventional brain function mapping to diagnose brain function, is not only widely used due to radioactive exposure but also secondary to brain activity. Therefore, it cannot be said that the brain function is accurately determined to reflect a change in blood flow. Similarly, functional MRI (hereinafter fMRI) strongly reflects changes in blood flow outside the brain, so accurate measurement of brain function was difficult. A problem in the measurement of brain function NIRS is that the region characteristics of brain function images change depending on the selection of an index reflecting brain activity.

  (2) In many brain function NIRS studies, OxyHb increase and DeoxyHb decrease have been used as a conventional pattern of brain activity due to neural activity. However, how often the pattern of conventional brain activity occurs, what kind of index can be distinguished between the center and the periphery of the activity, and how often the center and the periphery of the activity can be separated It was not even subject to measurement.

  (3) Even if the conventional brain activity pattern occurs, the ratio of OxyHb increase and DeoxyHb decrease is not always constant, so it should be distinguished for each task.

  (4) In the invention of Patent Document 1 of the present inventor, since there was no concept of brain habit itself, “measurement of brain habit” was not at all targeted.

  (5) A brain measurement index reflecting brain habits was discovered, brain habits were diagnosed, and there was no brain habit improvement instruction device for the brain habits.

  The present invention has been made in order to solve the above-described problem, and provides a brain habit diagnosis apparatus and program for diagnosing or instructing to improve brain habits, which are states of brain activity associated with living body habits. The purpose is to provide.

The brain habit diagnosis apparatus according to the present invention includes: a light irradiation light receiving unit that irradiates a predetermined part of the brain related to brain activity associated with a living body habit, and receives light from within the living body; and the light irradiation light receiving unit. A brain habit for diagnosing a brain habit which is a state of brain activity associated with the habit of the living body using a near-infrared spectroscopic method. A diagnostic device,
The apparatus main body has a concentration change amount Δ [HbO ₂ ] of oxidized hemoglobin, a concentration change amount Δ [Hb] of deoxygenated hemoglobin, and a concentration change of oxidized hemoglobin based on light information from the light irradiation light receiving means. A calculation means for calculating various parameters derived on the basis of the concentration change amount of the amount Δ [HbO ₂ ] and the deoxidized hemoglobin Δ [Hb]; and a change amount Δ [HbO of the oxidized hemoglobin calculated by the calculation means. ₂ ], a determination means for determining the presence / absence or degree of brain habit based on the change Δ [Hb] of deoxygenated hemoglobin and various parameters, and the presence / absence or degree of brain habit determined by the determination means Display means,
It is characterized by this.

The calculation means calculates the concentration change ΔCBV of the blood volume in the predetermined part of the brain by the equation (1),
The concentration change amount ΔCOE of the oxygen exchange amount in a predetermined part of the brain is calculated by Equation (2),
In the two-dimensional diagram showing the relationship between the concentration change ΔCBV of the blood volume and the concentration change ΔCOE of the oxygen exchange amount, a scalar of the trajectory vector from the starting point of the measurement obtained by plotting the problem over time Calculate L value and phase K angle,
The determination means may determine the presence or absence of a brain habit and the degree based on the scalar L value and the phase K angle.
ΔCBV = Δ [Hb] + Δ [HbO ₂ ] ··· Equation (1)
ΔCOE = Δ [Hb] −Δ [HbO ₂ ]... Formula (2)

The calculating means uses a first axis as the blood volume concentration change amount ΔCBV, a second axis orthogonal to the first axis as the oxygen exchange amount concentration change amount ΔCOE, and the first axis as the first axis. On the other hand, the third axis set to −45 degrees is defined as the concentration change Δ [HbO ₂ ] of the oxidized hemoglobin, and the fourth axis set to +45 degrees with respect to the first axis Dividing into eight phase regions on the two-dimensional diagram with the concentration change Δ [Hb] of oxidized hemoglobin, and calculating the frequency for each phase region where the trajectory obtained by plotting the task over time is present And
The determination means may determine the presence or absence or degree of brain habits according to the frequency for each phase region.

  The calculation means may calculate a phase distribution rate, a highest frequency phase, and an average phase distribution rate for each phase region.

  The calculation unit may calculate the frequency and ratio of the phase region for each of a plurality of tasks.

The calculation means calculates the oxygen saturation Y on the OD plane from a two-dimensional diagram in which the horizontal axis represents the oxidized Hb amount in the region of interest ROI (O) and the vertical axis represents the deoxidized Hb amount in the ROI (D). May be calculated by the equation (3).
Oxygen saturation Y = 1 / (1 + tan (Y angle)) ··· Equation (3)

The calculation means uses the amount of oxidized Hb in the region of interest ROI (O) on the horizontal axis and the amount of deoxidized Hb in ROI (D) on the vertical axis, and the amount of change in oxidized Hb (ΔO) and deoxidation. From the two-dimensional diagram showing the relationship between the amount of change (ΔD) in the mold Hb, the change ΔY in the oxygen saturation Y is calculated by the equation (4) using the change ΔY angle in the tilt Y angle on the OD plane. But you can.
Change in oxygen saturation ΔY = -1 / (1 + sin2 (ΔY angle)) Equation (4)

The calculating means may calculate the oxygen adjustment area ORA by the equation (5) from the scalar L value of the orbital vector from the starting point of the measurement and the phase K angle.

  The determination means may determine whether there is a mouth breathing habit or a nasal breathing habit by comparing the case where the oral cavity of the living body is blocked with the case where the nasal cavity is blocked.

  The determination means may determine a habit of using words by letting the living body hear the words.

  The determination means may determine exercise habits by giving an object to the living body.

  You may further have the improvement support means which diagnoses the said biological body at a different time, and judges whether the brain habit has improved based on the change of the determination result of the brain habit by the said determination means.

  The improvement support means may determine an instruction for making a good brain habit based on a determination result of the brain habit by the determination means.

The light irradiation / light receiving means may be installed at a site corresponding to a brain classification address related to a brain habit to be diagnosed.
The system further includes one or more second light irradiation / light receiving means for irradiating light on a part other than the brain of the living body and receiving light from a part other than the brain of the living body, It is also possible to diagnose brain habits by simultaneously measuring
The part other than the brain is, for example, the chest of the living body.

  The program of the present invention is characterized by causing the apparatus main body of the brain habit diagnosis apparatus to execute processing.

  According to the present invention, the following excellent effects can be obtained.

  (1) The experience of the individual who is a living body and the brain of the individual, and what kind of brain state each part (per brain address) is considered as “brain habits” By setting the concept and word “brain habit” as a way of looking at the brain, brain habits can be qualitatively and quantitatively diagnosed.

  (2) Brain habits can be diagnosed from brain measurement indices reflecting brain habits, and instructions for improving brain habits can be given based on the brain habit data for quantitatively and objectively changing habits.

  (3) Since the degree of brain habits is determined for each theme (for example, breathing habits, listening habits, watching habits, speaking habits, thinking habits, etc.), whether or not habits have been acquired and what habits are strong It is possible to objectively determine whether or not there is a bad brain habit and good brain habit. As a result, bad habits of the brain can be changed into habits. Moreover, since it can measure quantitatively that the brain is growing for every theme, it can be trained again that it is not made a habit.

  (4) The person can be aware of the actual feeling learned objectively, and can grasp the current situation and take measures. For example, when it comes to brain education, there is no way to determine whether you have learned or not. It is possible to judge how many words you know and how often you use them, and how you use them.

  (5) For example, regarding the exercise habits of the brain, it can be determined whether there is a habit of having heavy luggage or a habit of not having heavier objects than a ballpoint pen on a daily basis. Since the distortion of the body changes the habit of stimulating the brain, it is possible to determine the habit of biting on the right and biting on the left.

  (6) By changing your breathing habits, your brain can work better. It is possible to determine whether or not the breathing habit has been improved by the dental treatment, and that the habit remains in the brain after the treatment. Differences in habitual breathing and nasal breathing can be distinguished from the brain, and during correction, the condition after correction can be observed.

  (7) When brain habits are determined and habits are changed, an objective indicator of whether the brain actually changes can be made.

  (8) By quantifying the lifestyle status in the brain and comprehensively considering it, a comprehensive brain habit score such as breathing habits, listening habits, watching habits, speaking habits, thinking habits can be calculated.

  (9) The lifestyle is related to aging, diabetes, etc., and the lifestyle is reflected in the brain habit. By measuring the brain habit, it is possible to determine and give an instruction to lead a life in which the brain is easy to grow.

  (10) It is possible to present a new brain habit index that reflects the cerebral oxygen regulation reaction CORE, which was not found in conventional brain function imaging, and its utilization method.

  (11) Phase K imaging, phase difference intensity imaging of scalar L values, and new brain function imaging reflecting both K phase difference and L intensity difference can be shown. In PET and fMRI, brain activity indicators measure intensity differences, but none of these indicators show phase differences.

  (12) It is possible to easily measure brain habits by monitoring the cerebral oxygen regulation reaction, which has been difficult with PET, fMRI, and brain function NIRS, which have been used for conventional brain function mapping.

  (13) It is possible to make a brain measurement target from a different viewpoint of the difference in brain activity that occurs in brains of different shapes among individual brains.

  (14) Since the brain habit is the state of the subject's brain before the measurement, the state of the subject's brain before the measurement can be detected from the measurement technique that captures the brain function reaction performed on the spot.

It is a block diagram which shows the structure of the brain habit diagnosis apparatus which concerns on the example of embodiment of this invention. It is explanatory drawing which attached | subjected the address for brain division to the schematic of the brain of a biological body. It is explanatory drawing which shows the position of the address according to 12 functional systems. It is a flowchart for demonstrating the procedure for diagnosing a brain habit. It is a two-dimensional diagram for explaining a CORE vector, a K angle, and an L value. It is explanatory drawing which shows eight phase divisions of the CORE vector which orbits on a two-dimensional diagram. (A)-(C) are schematic diagrams of ΔCBV and ΔCOE in capillaries. It is explanatory drawing for demonstrating the model of the locus | trajectory of a CORE vector depending on a phase K angle. It is explanatory drawing for demonstrating the model of the locus | trajectory of a CORE vector depending on a phase K angle. It is explanatory drawing for demonstrating the model of the locus | trajectory of a CORE vector depending on a phase K angle. It is explanatory drawing for demonstrating the model of the locus | trajectory of a CORE vector depending on a phase K angle. It is a schematic diagram for explaining an area ORA defined on the CORE vector plane. It is explanatory drawing for demonstrating arrangement | positioning of the biological probe in the case of diagnosing an exercise habit. (A) shows the phase distribution of CORE in M1 (primary motor area), and (B) shows the phase distribution of CORE in the peripheral part of M1. It is the table | surface which showed the value of the average phase K and the average L which the average CORE in a task shows in M1 (1st motor area) and a peripheral region. It is explanatory drawing which shows the phase difference image of K in a subject, the intensity | strength difference image of L, and the image of ORA. (A) is a schematic explanatory diagram when diagnosing nasal breathing habits, and (B) is a schematic explanatory diagram when diagnosing mouth breathing habits. (A) is a graph showing fluctuations in blood volume during mouth breathing and time-series fluctuations in oxygen exchange amount, (B) shows fluctuations in blood volume during nasal breathing, and time-series fluctuations in oxygen exchange amount Graph (C) is a graph showing the trajectory plotted on the CORE vector plane. It is a table | surface which shows K angle and L value at the time of nasal breathing of four test subjects and mouth breathing. It is a graph which shows the phase change at the time of nasal breathing of two test subjects, and mouth breathing. The trajectory on the two-dimensional diagram during mouth breathing and nasal breathing due to the difference in brain part is shown, (A) is the case of the inner BA10 of the frontal lobe, and (B) is the case of the outer BA10 of the frontal lobe. Shows the trajectory on the two-dimensional diagram of mouth breathing and nasal breathing by persons with different breathing habits, (A) for nasal breathers (average of 7 people), (B) for those with strong mouth breathing habits, (C) is a case of a mouth breathing habit person. (A) shows the locus on the two-dimensional diagram when the language understanding phase is shown, and (B) shows the locus on the two-dimensional diagram when the language non-understanding phase is shown. (A) is a graph showing the relationship with oxygen saturation Y, with the horizontal axis representing the amount of oxidized Hb in ROI (O), the vertical axis representing the amount of deoxidized Hb in ROI (D), and (B) A graph showing the relationship between the amount of change in oxidized Hb in ROI (ΔO), the amount of change in deoxidized Hb in ROI (ΔD), and the amount of change in Y angle (ΔY angle), (C) on the OD plane It is a graph for demonstrating that the change of oxygen saturation is determined by the coordinate change ((DELTA) O and (DELTA) D).

S: Brain habit diagnosis device 1: Probe for living body (light irradiation / light receiving means)
DESCRIPTION OF SYMBOLS 1a: Light emitting element 1b: Light receiving element 2: Apparatus main body 3: Light quantity adjustment part 4: Selection adjustment part 5: Signal amplification part 6: A / D conversion part 7: Control part 8: Storage part 9: Display part 10: Calculation part 11: Determination unit 12: Improvement support unit 13: Program

Hereinafter, embodiments of the present invention will be described with reference to the drawings.
(Outline of brain habit diagnosis device)
FIG. 1 is a block diagram showing a configuration of a brain habit diagnosis apparatus according to an embodiment of the present invention.

  As shown in FIG. 1, the brain habit diagnosis apparatus S according to the embodiment of the present invention irradiates light to a predetermined part of the brain related to brain activity associated with the habit of the living body, and receives light from within the living body. It has a living body probe 1 (light irradiation / light receiving means) and an apparatus main body 2 that inputs light information from the living body probe 1 and performs calculation, control, or storage. This is a diagnosis of brain habits, which is the state of brain activity associated with the habits.

  The living body probe 1 includes at least two or more light emitting elements (light emitting diodes) 1a for irradiating light to a predetermined measurement portion of the brain related to the brain activity associated with the habit of the living body, and transmitted light and reflected light from the measurement portion. Or at least 2 or more light receiving element (photodiode) 1b ... which receives the light after interacting with living bodies, such as scattered light, is comprised.

  The distance between the light emitting element 1a and the light receiving element 1b is about 1.5 cm to 3 cm.

  Further, in the biological probe 1, it is preferable to arrange a plurality of distances (in a matrix) between the light emitting element 1a and the light receiving element 1b at equal intervals, but the distance between the light emitting element 1a and the light receiving element 1b is not necessarily equal. Random arrangement is also possible without having to be an interval.

  The apparatus main body 2 includes a light amount adjustment unit 3 that adjusts the light emission amount of the light emitting elements 1a..., A selection adjustment unit 4 that selectively enables or disables arbitrary light receiving elements 1b. The signal amplifying unit 5 capable of gain control for amplifying the signals from the light receiving elements 1b..., The A / D converting unit 6 for converting the output of the signal amplifying unit 5 into numerical values, the control processing of each unit and the A / D converting unit 6 A control unit 7 that executes predetermined arithmetic processing based on the output, a storage unit 8 that is used to store the output of the A / D conversion unit 6, control data of each unit or calculation results, and the A / D conversion unit 6 and a display unit 9 for performing display based on the output results and the calculation results.

  The apparatus main body 2 may have a printing function for printing various data and a communication function for transmitting and receiving data via a communication network.

Based on the optical information from the biological probe 1, the control unit 7 changes the concentration change amount Δ [HbO ₂ ] of oxidized hemoglobin, the concentration change amount Δ [Hb] of deoxidized hemoglobin, and the concentration change amount of oxidized hemoglobin. A calculation unit 10 that calculates various parameters derived based on the concentration change amounts of Δ [HbO ₂ ] and deoxidized hemoglobin Δ [Hb], and an oxygenated hemoglobin concentration change amount Δ [HbO calculated by the calculation unit 10 ₂ ], a determination unit 11 that determines the presence or absence and degree of brain habit based on the change Δ [Hb] of deoxygenated hemoglobin and various parameters, and a brain habit by the determination unit 11 by diagnosing the living body at different times And an improvement support unit 12 that determines whether or not the brain habit has improved based on the change in the determination result.

  Moreover, the improvement support part 12 determines the instruction | indication for setting it as a favorable brain habit based on the determination result of the brain habit by the determination part 11, and the display part 9 displays the content of the instruction | indication.

  Two types of light emitting elements 1a ... of the biological probe 1 are prepared, one that emits light with a wavelength of 730 nm and one that emits light with a wavelength of 850 nm (note that the numerical value of the wavelength of this light is an example. It is not limited to the above, and three or more wavelengths may be combined). These are arranged alternately, for example, in the column direction. However, when considering other patterns, the arrangement should be such that the amount of received light can be measured in a balanced manner in consideration of the attenuation depending on the wavelength in the tissue. is important. All the light emitting elements 1a... Are connected to the light amount adjusting unit 3 of the apparatus main body 2, and the light emission amount can be adjusted as a whole or independently.

  On the other hand, all the light receiving elements 1b... Are connected to the signal amplifying unit 5 via the selection adjusting unit 4 of the apparatus body 2, and all or a part of the light receiving signals output from the respective light receiving elements 1b are selected. The signal is output to the signal amplifying unit 5 while being selected and adjusted by the adjusting unit 4 and is amplified there. The amplified received light signal is digitized by the A / D converter 6 and output to the controller 7.

  The control unit 7 applies a low-pass filter to the digital data input from the A / D conversion unit 6 to perform noise removal processing, and then stores the processing data (hereinafter referred to as “light reception amount”) in the storage unit 8 in a time table. Remember.

Moreover, the control part 7 performs the arithmetic processing demonstrated below based on the obtained received light quantity. First, the absorbance (OD ₇₃₀ ) at a wavelength of ₇₃₀ nm is calculated by the equation (6), and the absorbance (OD ₈₅₀ ) at a wavelength of 850 nm is calculated by the equation (7). Is stored in the storage unit 8.

O.D. ₇₃₀ = log 10 (I _{0 730} / I ₇₃₀ ) (6)
OD _.850 = log 10 (I _{0 850} / I ₈₅₀ ) (7)
I _{0 730} : Amount of light emitted at a wavelength of ₇₃₀ nm I ₇₃₀ : Amount of received light at a wavelength of ₇₃₀ nm I _{0 850} : Amount of emitted light at a wavelength of 850 nm I ₈₅₀ : Amount of received light at a wavelength of 850 nm Here, the concentration change amount of oxidized hemoglobin and deoxidized hemoglobin It is known from the known theory that there is a relationship of Equation (8) and Equation (9) between the concentration change amount and the absorbance change amount.

ΔO.D. ₇₃₀ = a1Δ [HbO ₂ ] + a1′Δ [Hb] (8)
ΔO.D. ₈₅₀ = a2Δ [HbO ₂ ] + a2′Δ [Hb] (9)
ΔO.D. ₇₃₀ : Change in absorbance at a wavelength of ₇₃₀ nm ΔOD. ₈₅₀ : Change in absorbance at a wavelength of 850 nm Δ [HbO ₂ ]: Change in concentration of oxidized hemoglobin Δ [Hb]: Change in concentration of deoxyhemoglobin a1, a1 ′, a2, a2 ′: Absorbance coefficient Accordingly, Equations (10) and (11) are obtained from this known simultaneous equation.

Δ [HbO ₂ ] = a {ΔO.D. ₇₃₀ − (a1 ′ / a2 ′) ΔO.D. ₈₅₀ } Expression (10)
Δ [Hb] = a (a 2 / a 2 ′) {(a 1 / a 2) ΔO.D. ₈₅₀ −ΔO.D. ₇₃₀ } (11)
a = a2 ′ / (a1a2′−a1′a2) ≈1 (1 or its neighboring value)
Therefore, in terms of calculated amount of change in absorbance wavelength 730nm (ΔO.D. ₇₃₀₎ and the amount of change in absorbance wavelength 850nm (ΔO.D _.850), concentration variation of oxidized hemoglobin (Δ [HbO _2]) The deoxyhemoglobin concentration change amount (Δ [Hb]) is calculated by the equation (10) and by the equation (11), and the calculation result is stored in the storage unit 8 in a time table. Note that the amount of change in total hemoglobin concentration (Δ [total-Hb]) is expressed by equation (12).
Δ [total-Hb] = Δ [HbO ₂ ] + Δ [Hb] (12)

  By the way, the change mode of each concentration change amount of the oxidized hemoglobin and the deoxidized hemoglobin in the capillary induced by stimulation to the tissue shows the following 9 patterns depending on the combination of the increase and decrease.

(1) Δ [HbO ₂ ] increase Δ [Hb] increase (2) Δ [HbO ₂ ] increase Δ [Hb] decrease (3) Δ [HbO ₂ ] increase Δ [Hb] zero (4) Δ [HbO ₂ ] Decrease Δ [Hb] increase (5) Δ [HbO ₂ ] decrease Δ [Hb] decrease (6) Δ [HbO ₂ ] decrease Δ [Hb] zero (7) Δ [HbO ₂ ] zero Δ [Hb] increase (8 ) Δ [HbO ₂ ] Zero Δ [Hb] decrease (9) Δ [HbO ₂ ] Zero Δ [Hb] Zero Actually, the metabolic activity of the tissue depends on the stimulus application conditions and the physiological state of the resting state. The above pattern changes over time. Δ [Hb] and Δ [HbO ₂ ] of the capillary fluctuate as bloodstream metabolic activity for taking oxygen into the tissue from oxidized hemoglobin in the capillary.

Therefore, in the present invention, in order to diagnose the brain habits of a living body, various parameters derived based on the concentration variation Δ [HbO ₂ ] of oxidized hemoglobin and the concentration variation Δ [Hb] of deoxidized hemoglobin are determined. Calculation is performed by the calculation unit 10 of the control unit 7.
(Body probe location and brain classification address)
The biological probe 1 is preferably installed at a site corresponding to a brain classification address related to the brain habit to be diagnosed.

  Here, the address for brain classification (the inventor and the applicant refer to the address as the brain address (registered trademark)) refers to a change in nerve activity according to a predetermined function and role in the brain. It is a number assigned separately for each part considering the increase in oxygen consumption and the function of regulating brain activity by activation and sedation.

  FIG. 2 is an explanatory diagram in which the brain classification address is added to the schematic diagram of the living body brain.

  The brain classification address is a word that represents each address of the brain map indicating the division of roles of the brain. The brain is finely named by cell type and groove shape. There are several other places where a general role is granted. However, there was no brain address classification including the concept of time that the brain grows while changing its shape.

  Therefore, the present inventor found that the brain is divided into 120 parts and changes in shape and function as it grows, and is named “brain address” so that even ordinary people can understand.

The brain classification addresses are roughly divided into 12 systems.
FIG. 3 is an explanatory diagram showing the positions of 12 functional system addresses.
・ Thinking brain addresses (brain addresses 10, 8, 9): Thinking and thinking ・ Planning / execution brain addresses (left and right brain addresses 6)
： Responsible for planning things and moving the body ・ Exercise brain address (Brain Address 4, 6): Use when moving the body ・ Communication brain address (Brain addresses 44, 45, 46)
: Communicates through speech, word manipulation and communication ・ Emotional brain address (brain address 11, amygdala): Relation to likes and dislikes and emotions ・ Social brain address (brain addresses 11, 10): For human interaction and position Consideration / understanding system brain address (brain number 22, 39, 40): Understanding things through words / sensory system brain address (brain number 3): Touching and touching with hands and skin / Memory system brain address (hippocampus) : Remember or remember ・ Knowledge system brain address (Brain address 20, 38): Stores the knowledge that is understood ・ Hearing system brain address (Brain addresses 41, 42, 22, 21)
): Listening to sounds using the ears ・ Visual system brain addresses (brain addresses 17, 18, 19, 7, 8, 9): Seeing using the eyes For the invention of this brain classification address, see this book Already filed and patented by the inventor (see Japanese Patent No. 4838845 for details).

(About main measurement parts of each brain habit)
Various brain habit measurement sites, that is, the installation site of the biological probe 1 are as follows based on the above-described brain classification addresses.
・ In the case of exercise habits, it includes 4 and 6 of the left and right brain addresses ・ In the case of thinking habits (thinking habits) and image thinking habits, it includes the left and right brain addresses 10, 8, and 9 ・ In the case of social habits In the case of hearing habits (listening habits), in the case of auditory habits (listening habits) In the case of sensory habits (touching habits), the left and right brain addresses 3 In the case of auditory habits (listening habits), the left and right brain addresses 41, 42, 22, 21 are included. In the case of breathing habits, the left and right brain addresses 10, 44, 45, 46, 6, 4, are included. Including in the case of language manipulation, speech recall and conversation habits (speaking habits), include left and right brain addresses 44, 45, 46. In the case of visual habits (watching habits), the left and right brain addresses 17, 18, Includes 19, 7, 8, 9 ・ Language understanding habits (word usage habits) In the case of image understanding habits, it includes left and right brain addresses 22, 39, 40. In the case of game habits, it includes left and right brain addresses, 17, 18, 19, 7, 4, 6. In the case of reading weeks, it is left and right. In the case of planning / execution habits, the left and right brain addresses 6 are included. In working memory habits, the left and right brain addresses 44, 45, 46, 39 are included. , 40 included

(About brain habits)
Here, as a new view of the brain, we define the concept and word “brain habits”.
“Brain habit” means the result of the subject (living body) before the measurement having experienced or worn it, or conversely, the state of the brain not experienced or unlearned. Since brain habits are not positive and negative, and there is an intermediate state between them, in practice, the level of brain habits can be regarded as on the spectrum.

  The brain habits are “the shape of the brain address” as to how far the shape of the brain address has grown for each brain address described above, and how much the corresponding brain address has been used in a certain task for each brain address. This can be understood from the two viewpoints of “experience in brain address”.

  The “brain address shape” can be obtained by creating a white matter-enhanced image and analyzing the branching of each brain address (see Patent Document 2). Depending on the stage of growth of the brain branching, if it is digitized and the degree of maturity is high, it can be determined that there is a lot of experience at the corresponding brain address.

  Even if scoring the form of “brain address” for brain habits, it is impossible to determine what information the form collects at that address or what information it was exposed to, so it can not be determined by the form of the brain address, so the brain response is measured To understand.

  “Experience of brain address” determines what kind of experience and experience the branching that grew for each brain address actually grew, and conversely, what kind of experience and experience did not experience, There was no device to train the brain address.

  Strictly speaking, brain habits can be measured for each of about 120 brain addresses, but can be divided into 12 functional systems for convenience, or brain addresses can be combined for the desired habituality.

  Specific examples of brain habits, “opposite unexperienced, unlearned brain state” are as follows.

  In the visual system, for example, a face photograph of a person who has never seen it, the first painting, the state of the brain when nothing is felt even when looking at the figurine, and the paintings that have met and curious many times The impression is completely different when you see it. This depends on how well the brain habits are made with human faces and paintings.

  In the auditory system and comprehension system, for example, Kanji that I have never seen and foreign languages I have never learned (Swahili, English, French, Spanish) Even if I grew up in Japan, I can't understand the meaning even if I listen to an immature language. On the other hand, even if you hear the same thing in your native language and Japanese, you can even understand.

  In the exercise system, for example, if a person who strikes a computer with two fingers on the left and right starts to strike with 10 fingers on the left and right, not only will the speed slow down, but it will also be annoying. This is not a brain habit using a PC with 10 fingers. On the other hand, the person who usually works with 10 fingers of a PC can hit it without looking at the keyboard. This is the difference between those who have exercise brain habits and those who do not.

  In the sensory system, for example, everyone closes their eyes and touches a spoon. However, if you are touching a giraffe made of sponge, you will not know until you are used to it. It is a difference in sensory habits that determines the shape by touching.

(Diagnosis procedure for brain habits)
Since brain habit sensitive indicators are also sensitive to brain growth, when oxygen is consumed in the brain, brain blood volume is also changed and adjusted as an indicator of brain activity. , Cerebral oxygen regulation (CORE) vector is set as the measurement target.
In other words, how the cerebral oxygen regulation (CORE) reaction changes, and since it was not the target of measurement by quantitative indicators, it is measured by the vector NIRS method using a two-dimensional vector plane.

  In brain function mapping that has been widely performed, tasks were simply loaded to extract the in-situ responses, but the conventional view was abandoned and all responses to the tasks were examined. Once again, it is understood from the perspective of how the brain has grown and the state of the brain has been created. In other words, brain activation is not a good activity or not, but a reaction pattern that occurs as a result of brain habits. For this purpose, the frequency and probability of two or more tasks are compared and analyzed.

  That is, “If the brain habits change, the probability distribution, phase distribution K angle, and intensity distribution (L, four Hb changes component vector) of the brain reaction change”.

Considering the trajectory of the CORE vector as a result of brain habits, the change in the trajectory of the vector CORE is evaluated as a change in brain habits.

  The vector NIRS method, which captures CORE, can perform real-time imaging of CORE, which differs for each part of the subject from the L intensity difference and K angle phase difference, from the phase difference and intensity difference.

  How often does it occur on the two-dimensional plane for the task, what kind of index can distinguish between the center and the periphery of the activity, and how often can the center and the periphery of the activity be separated? Check out.

  Diagnose brain habits and give instructions based on the comparison of two tasks and the dynamics of indicators such as phase distribution probability, phase frequency, average phase, and maximum frequency phase for two or more tasks that have the same task but with different loads.

  FIG. 4 is a flowchart for explaining a procedure for diagnosing brain habits.

First, the concentration change amount Δ [HbO ₂ ] of oxidized hemoglobin and the concentration change amount Δ [Hb] of deoxygenated hemoglobin are calculated by the calculation unit 10 (step S1).

Next, the concentration change ΔCBV of the blood volume at a predetermined part of the brain is calculated by equation (1),
The concentration change amount ΔCOE of the oxygen exchange amount in a predetermined part of the brain is calculated by the equation (2) (step S2).
ΔCBV = Δ [Hb] + Δ [HbO ₂ ]... Formula (1)
ΔCOE = Δ [Hb] −Δ [HbO ₂ ]... Formula (2)

Next, in the two-dimensional diagram showing the relationship between the blood concentration change ΔCBV and the oxygen exchange concentration change ΔCOE, the scalar vector of the trajectory vector from the starting point of measurement obtained by plotting the problem over time The L value and the phase K angle are calculated (step S3). Here, the phase K angle = Arctan (ΔCOE / ΔCBV).
Next, the determination unit 11 determines whether or not there is a brain habit based on the scalar L value and the index of the phase K angle (step S4).

  Thereafter, the improvement support unit 12 determines whether or not the brain habit is improved based on the change in the determination result of the brain habits by the determination unit 11, or the good brain habits based on the determination result of the brain habits. The display unit 9 displays the determination result of the determination unit 11 and the content of the instruction of the improvement support unit 12 (step S6).

(About definition of CORE vector, K angle, L value)
FIG. 5 is a two-dimensional diagram for explaining the CORE vector, K angle, and L value.
As shown in FIG. 5, the ΔO axis (change in the concentration of oxidized hemoglobin), the ΔD axis (the change in concentration of deoxygenated hemoglobin), the ΔCBV axis (the change in blood volume), the ΔCOE axis (oxygen) The polar coordinate plane consisting of the four axes of the change amount of the exchange amount is defined as a cerebral oxygen regulation vector plane (CORE vector plane).

  △ CBV axis is set to the horizontal axis, △ COE axis is set to the vertical axis perpendicular to the horizontal axis, △ O axis is set to -45 degrees with respect to △ CBV axis, △ D axis is set to △ CBV axis It is set to +45 degrees.

The CORE vector has the property of K angle indicating amplitude L and phase. The positive ΔCBV axis is assumed to be 0 degree, the intensity of oxygen metabolism is quantified as L value, and the phase is quantified as oxygen exchange degree K angle.
(Eight phase sections of the CORE vector trajectory on the two-dimensional diagram)
FIG. 6 is an explanatory diagram showing the eight phase segments of the CORE vector trajectory on the two-dimensional diagram.

  As shown in FIG. 6, the octant on the vector plane is divided into eight qualitative phase ranges (-4 to 4). Each phase region indicates the difference in the increase / decrease ratio of the CORE vector components (ΔO, ΔD, ΔCBV, ΔCOE).

  The K angle representing the degree of oxygen exchange and the eight phase sections are indicators for detecting changes in blood volume and oxygen regulation in the blood vessel within a certain time.

  Phases 1 to 4 (0 ≦ K ≦ 180) indicate the direction of hypoxia (ΔCOE increase). Phases −1 to 3 (−45 ≦ K ≦ 135) indicate deoxygenation (ΔD increase). Phases 2 to 2 (-90 ≦ K ≦ 90) indicate an increase in blood volume (ΔCBV increase). Phases -3 to 1 (-135 ≦ K ≦ 45) indicate an increase in arterial blood inflow (ΔO increase). Phases -4 to -1 (-180 ≦ K ≦ 0) indicate hyperoxygenation (ΔCOE decrease). Phases 3 to 4 and -3 to -4 (135 ≦ K ≦ 180, −180 ≦ K ≦ −135) indicate a decrease in blood volume (decrease in ΔCBV). For example, phase 1 (0 ≦ K ≦ 45) indicates a state in which the increase in deoxygenation is larger than the increase in inflow of arterial blood and the increase in blood volume is larger than the increase in hypoxia.

(Schematic diagram of ΔCBV and ΔCOE in capillaries)
7A to 7C are schematic diagrams of ΔCBV and ΔCOE in the capillary blood vessels.
As shown in FIG. 7A, when ΔCBV does not change, when nerve cells are activated and the demand for oxygen increases, oxygen moves from within the capillaries and ΔCOE increases, so the K angle also increases. . On the other hand, if there is no oxygen demand, as shown in FIG. 7B, ΔCOE decreases and ΔKV decreases as ΔCBV increases. It can be assumed that the center of the brain activity is closer to the state of FIG. 7C, and the state of FIG. Actually, activation (e) reflecting all ratios of CBV and COE by the mixture of FIG. 7 (A) and FIG. 7 (B) can be quantitatively evaluated using the K ratio.
(Model of CORE vector locus depending on phase K angle)
8 to 11 are explanatory diagrams for explaining a model of the locus of the CORE vector depending on the phase K angle. FIG. 8A shows the phase region 1, FIG. 8B shows the phase region 2, and FIG. 9A is the phase region 4, FIG. 9B is the phase region 4, FIG. 10A is the phase region 1, FIG. 10B is the phase region 2, and FIG. FIG. 11B shows the phase region -4, respectively.

  The time series model of four vector components is shown in contrast with the two-dimensional locus of CORE. The one-dimensional time-series data of all indices is always reflected in eight phases in two dimensions depending on the K angle. In actual measurement, ΔCBV (black) and ΔCOE (green) can be obtained from changes in ΔO (red) and ΔD (blue) and converted to CORE. Each vector component and each circle on CORE indicate the same time point. In each phase model, a is a trajectory with a constant K angle and L changes (red line), b is a trajectory with a constant K angle and L peak is half (thick line), and c is a model in which K changes and the CORE trajectory swells (dashed line) ) Respectively.

When CORE swells and orbits like c as compared with a, even if the peak of the L value is constant, the K angle is constant for a, but the K angle changes for c. In a, the point in time at which the maximum value of L is shown is the same as the point in time at which the maximum change value of the one-dimensional model is shown. However, in c, the time point indicating the maximum value of the L value is different from the time point indicating the maximum change of each index. The CORE model shows a one-dimensional waveform when a rotational motion occurs clockwise within each phase. In the counterclockwise rotational movement, the time change of the one-dimensional waveform occurs in reverse. A difference in the K angle indicates a phase difference. A difference in L value indicates a difference in intensity.
(Schematic diagram of area ORA defined on CORE vector plane)
FIG. 12 is a schematic diagram for explaining the area ORA defined on the CORE vector plane. As shown in FIG. 12, L obtained from the coordinates of an arbitrary moving point obtains the area swept by K angles from the positive ΔCBV axis. The area is also positive or negative depending on the sign of the K angle. The positive change shown in ORA1 is defined as the oxygen exchange increasing direction, and the negative change shown in ORA2 is defined as the oxygen exchange decreasing direction (that is, oxygen passing direction). ORA reflects how much the CORE trajectory has moved positively or negatively.

In imaging using the phase K and the amplitude L, an image showing the difference in phase was created using the K angle. From the viewpoint of observing the rotational motion around the origin at the start of the task, the oxygen regulation area (ORA) on the CORE vector plane was defined as the following equation (5) using Lt and K angle radians Kradt. .

  ORA indicates the cumulative area of the CORE trajectory, where Lt after t seconds from the measurement start point (origin) is swept by Kt degrees from the ΔCBV axis at K = 0 °. When ORA is positive, it indicates that the area has expanded in the oxygen exchange increasing direction (0 <K), and when ORA is negative, it indicates that the area has expanded in the oxygen passage direction (K <0). , We imaged the ORA during the task (from t = 0s to t = 36s). Since the K angle component is also included in ORA, a phase difference image is shown. This is an ORA imaging that reflects both the phase difference and the intensity difference of activity in order to evaluate the subject's ORA variation.

(CORE phase distribution rate, highest frequency phase, average phase distribution rate)
The increasing response rate of the cumulative sum of each vector component in the task and the phase distribution rate of the eight phase sections indicated by CORE were calculated using the cumulative sum of the vector components in the task. The increase response rate was calculated from the ratio of the number of trials in which the cumulative sum of each vector component in the number of analysis trials of each task showed a positive value. The average increase response rate at the peripheral site was calculated by calculating the average increase response rate of each channel.

  The CORE phase distribution ratio was calculated from the ratio of the number of trials of each phase determined by the classification of the K angle of the single-trial CORE in the number of analysis trials of each task. Adding the eight phase distribution ratios gives 100%. The average phase distribution ratio in the peripheral region was calculated by calculating the average of the phase distribution ratio of each channel. Among the phase distribution ratios of M1 and surrounding parts, the phase showing the highest phase distribution ratio was defined as the “highest frequency phase” in that part.

  The difference between the increasing response rate of each vector component in the motor area M1 and the average increasing response rate of each vector component in the surrounding region was tested by two-dimensional layout analysis of variance and multiple comparison (Bonferroni). The difference between the phase distribution rate at M1 and the average phase distribution rate at the peripheral site was tested using an independent two-group T test for each phase.

  Here, the “CORE phase distribution ratio” defines eight quadrants on the CORE vector plane as qualitative eight phase sections (phases -4 to 4) indicating combinations of increase / decrease of four vector components.

  Calculation is made using the cumulative sum of each vector component in the task.

  The cumulative sum (ΔOcs, ΔDcs, ΔCBVcs, ΔCOEcs) of vector components at an arbitrary time is obtained as follows.

  △ O cumulative sum (△ Ocs)

△CBVの累積和（△CBVcs） △ D cumulative sum (△ Dcs)

△ Cumulative sum of CBV (△ CBVcs)

次に、単試行の任意時間帯でのCOREcsが示す位相ＫをArctan （△COEcs/△CBVcs）から求める。 △ Cumulative sum of COE (△ COEcs)

Next, the phase K indicated by COREcs in an arbitrary time zone of a single trial is obtained from Arctan (ΔCOEcs / ΔCBVcs).

[CORE phase distribution ratio] = 100x [number of trials of each phase indicated by single trial COREcs] / [total number of trials]
Therefore, the CORE phase distribution ratio was calculated from the ratio of the number of trials for each phase indicated by the single-trial COREcs in the number of analysis trials for each task. When the distribution ratio of the eight phase segments is added, it becomes 100%.

  The average of the phase distribution ratio of each channel can be obtained and calculated.

  The “average phase” is obtained from Arctan (ΔCOEcs / ΔCBVcs) as a phase K indicated by COREcs in an arbitrary time zone of a single trial.

The CORE vector plotted on the CORE vector plane using △ CBVcs and △ COEcs is COREcs. The total average of the phases indicated by COREcs for each trial was defined as “average phase”.
= [Number of trials of each phase indicated by single trial COREcs] x [Number of phases] / [Number of trials of all tasks]
“Average phase” = (1 × [Frequency of phase 1] + 2 × [Frequency of phase 2] + 3 × [Frequency of phase 3] + 4 × [Frequency of phase 4] + (− 1) × [Frequency of phase−1] + ( -2) x [frequency of phase-2] + (-3) x [frequency of phase-3] + (-4) x [frequency of phase-4]) / [total number of trials]
In addition, “the highest frequency phase” is the number of trials of each phase indicated by the COREcs of each phase single trial with respect to the total number of trials in each phase distribution ratio of the measured part.
[CORE phase distribution ratio] = 100 × [number of trials in each phase indicated by single-trial COREcs] / [number of trials of all tasks] indicates the phase showing the maximum frequency.

The phase distribution rate of COREcs has the highest frequency phase, which reveals the nature of brain habits that show a specific phase distribution depending on the region and task.
(About judgment of brain habits)
The presence or absence or degree of brain habits determined by the determination unit 11 can be determined by whether oxygen is efficiently consumed and supplied efficiently when brain cells work.

  A) When there is a brain habit, oxygen consumption occurs without increasing the blood volume in a brain regulation reaction with good oxygen consumption efficiency.

  B) When there is no brain habit, oxygen consumption does not occur even though the blood volume is increased by the regulation reaction of the brain with poor oxygen consumption efficiency.

  C) When the brain habits occur, the brain address changes from B) to A).

  That is, FIG. 7A is a state where brain habits are strong (high state), FIG. 7B is a state where there is no brain habits (low state), and FIG. It can be determined that the brain has no habits and is in the middle of a certain state.

  In regions where brain habits are strong (high), ΔCOE increases even if ΔCBV is constant, so K increases (see FIG. 7A). The closer the K angle approaches 90 degrees, the higher the amount of oxygen exchange with respect to the change in blood cell volume, which is thought to reflect the increase in oxygen consumption efficiency depending on the task. In this case, a reaction occurs in which the oxygen concentration in the capillaries decreases.

  On the other hand, as shown in FIG. 7 (B), ΔCOE does not increase even when ΔCBV increases due to the inflow of arterial blood with high oxygen saturation at a site where there is no brain habit with strong OxyHb passage response (low). The K angle decreases as the amount of oxygen exchange with increasing blood volume becomes negative. The closer the K angle is to -90 degrees, the higher the amount of oxygen that is not used for oxygen exchange (passing rate) relative to the amount of change in blood cell volume, the oxygen supply efficiency into the blood vessel depending on the task increases, and oxygen consumption This is thought to reflect the decline in efficiency. In this case, a reaction occurs in which the oxygen concentration in the blood vessel increases.

  As shown in FIG. 7C, the physiological significance of the K angle as the degree of oxygen exchange can quantitatively show all the ratios of ΔCOE fluctuations to ΔCBV fluctuations that vary depending on the strength of brain habits. is there. As the K angle is closer to 0 degrees, there is no change in ΔCOE, and ΔCBV is increasing. Conversely, the closer the K angle is to 180 degrees, there is no change in ΔCOE, indicating a reaction in which ΔCBV decreases.

  Therefore, △ COE increase is quantitatively evaluated as 0 <K angle <180, and △ CBV increase is -90 <K angle <90. It will be observed that brain habits will be established as they get closer.

  The L value or the four vector indicators that make up the L value, the greater the L value, the greater the change in Gron, and the inefficient use of energy by changing the blood volume. It can be determined that the higher the state, the lower the brain habit.

Actually, eight quadrants on the CORE vector plane are defined as eight qualitative phase sections (phases -4 to 4) indicating combinations of increase / decrease of four vector components. In the CORE vector plane, a positive ΔCBV indicates a ΔCBV increase (0 <ΔD + ΔO), and a negative value indicates a ΔCBV decrease (ΔD + ΔO <0). On the △ COE axis, a positive value (0 <△ D- △ O) indicating progress toward a hypoxic change is a negative value (△ D- △ O < 0) means △ COE decrease.
△ O is larger, lower brain habits △ D is larger, brain habits are higher △ CBV is larger, brain habits are lower △ COE is larger, brain habits are higher, so when combining L and K When there is no brain habit, L increases and K indicates hyperoxia, and when there is brain habit, L does not increase and K indicates hypoxia.

(If there are multiple issues)
When there are two tasks such as the frequency A of the A task and the frequency b of the B task for the plurality of tasks, the calculation unit 10
A brain habitual degree = a / (a + b) × 100 (%), B brain habitual degree = b / (a + b) × 100 (%).

  The determination unit 11 determines that the addictive displacement is strong when the degree of brain habit is close to 100%, and that the addictive difference is poor when it is close to 50%.

As the brain habit lateralization index, the calculation unit 10
Brain habit lateralization index = (a−b) / (a + b) is calculated, and when the habitual lateralization index is close to 1, −1, the determination unit 11 has stronger addictiveness and is close to 0 It is judged that the habitual difference is poor. In the case of 50% and 0, the determination unit 11 determines that both have low addictiveness if the low oxygen phase is strong, and both have strong addictiveness if not low oxygen.
In addition, the determination part 11 may compare 3 or more habits similarly.

(Instructions of the improvement support unit 12)
The improvement support unit 12 gives an instruction to perform the task B when, for example, the brain habit level of the task A is 75% or more, and issues an instruction to perform the task A when the degree is 25% or less. In the case of, instructions for half of the tasks A and B are issued.

  When the measurement is performed for the second time and thereafter, the calculation unit 10 calculates the degree of improvement of brain habit = (current degree of brain habit) − (previous degree of brain habit) or the rate of improvement of brain habit = (current degree of brain habit. -The previous brain habit level) / (previous brain habit level) is calculated, and the improvement support unit 12 issues an instruction based on those values.

Next, a specific diagnosis of brain habits will be described.
(In the case of exercise habits)
FIG. 13 is an explanatory diagram for explaining the arrangement of biological probes when diagnosing exercise habits. In FIG. 13, numerals (1 to 7) written between the irradiation probe (◯) that is the light emitting element 1a and the detection probe (●) that is the light receiving element 1b indicate channel numbers. The probe was installed so that 1ch corresponds to the left hemisphere M1 (precentral knob: *) (Murakoshi et al., 2006). The surrounding 2-7ch was arranged to surround 1ch. In the figure, precentral sulcus (preCS), central sulcus (CS), and postcentral sulcus (postCS) are shown.

For example, when deciding whether or not you have a habit of having heavy luggage or routinely having a habit of not having an object heavier than a ballpoint pen, prepare two or more types of dumbbells or luggage and lift them.
In the habit of not having a heavier object than the ballpoint pen, L becomes stronger and K indicates oxygen consumption. In such a case, the determination unit 11 is determined to have no exercise habits, and the improvement support unit 12 presents exercise training menus such as push-ups and abdominal muscles.
Even if a heavy object is lifted, the increase in L is poor, and if the phase K angle does not indicate hypoxia, it can be determined that exercise habits are strong.

(Phase frequency in exercise habit measurement)
The brain habit can be determined by one measurement, but a plurality of measurements may be tried to improve accuracy.

Here, the phase frequency refers to the cumulative frequency indicating the corresponding phase section in N trials. That is,
Phase frequency f = [Cumulative frequency indicating the relevant phase category] / [N trials]
For example, dumbbells are calculated for each phase change of 0 kg, 4.5 kg, and 9.5 kg (0 kg: 35 trials, 4.5 kg: 35 trials, 9.5 kg: 30 trials).

  FIG. 14A shows the phase distribution of CORE in M1 (primary motor area), and FIG. 14B shows the phase distribution of CORE in the peripheral part of M1.

  As shown in FIG. 14 (A), in M1, phase region 2 showing an increase in ΔCOE and an increase in ΔCBV is the highest at the highest frequency phase, 4.5kg at 54.3%, 9.5 kg at 46.7%, significantly higher than the surrounding area. (Significance probability p <0.01).

  On the other hand, as shown in FIG. 14 (B), in the peripheral region of M1, phase regions -1 and -2 showing an increase in ΔCBV without an increase in ΔCOE are widely distributed at the maximum high frequency phase, totaling 57.1% (4.5 kg) ), 67.7% (9.5 kg) in total, which was significantly higher than M1 in each phase range (significance probability p <0.05).

  That is, since the ΔCOE decreased with respect to the increase in ΔCBV in the peripheral region, it was shown that there is a high probability that ΔCBV and ΔCOE deviate.

  Here, since the subjects who are 4.5 kg and 9.5 kg and the phase is not positive and the minus is lifting oxygen without using oxygen, the determination unit 11 determines that the subject has exercise habits. Conversely, subjects with a high positive frequency are determined to have poor exercise habits.

  The improvement support unit 12 issues an instruction to acquire an exercise habit according to the above phase distribution.

  FIG. 15 is a table showing the values of the average phase K and average L indicated by the average CORE in the task in M1 (primary motor area) and the peripheral region. Here, ORA represents the value of the change during the task obtained from the above-described equation (5).

  As can be seen from FIG. 15, the K of M1 increased from 39.1 degrees to 56.6 degrees as the dumbbell load increased. On the other hand, the maximum high-frequency phase K decreased significantly from -32.7 degrees to -56.9 degrees in the surrounding area. As for the phase difference between the M1 and the K angle of the peripheral part, a significant difference of 105.7 degrees at 4.5 kg and 109.4 degrees at 9.5 kg was observed (p <0.01).

  A significant intensity difference between M1 and L in the surrounding region was observed only at 0 kg, not 4.5 kg and 9.5 kg. L increased significantly due to the increase in dumbbell weight in both M1 and the surrounding region. Compared to 0kg L of M1, the strength of L of the 9.5kg task increased significantly to 33.3 times and 3.5 times of 4.5kg. The intensity of L of the 9.5kg task in the surrounding area was significantly increased to 8.4 times 0kg and 2.4 times 4.5kg.

  ORA increased significantly with M1 at 4 kg and 9.5 kg from the surrounding area. In addition, an increase in dumbbell weight caused significant ORA changes in both M1 and surrounding sites. Since ORA contains both K and L, it was possible to emphasize the difference between the increase and decrease in the surroundings and M1.

FIG. 16 is an explanatory diagram showing a K phase difference image, an L intensity difference image, and an ORA image in question.
As shown in FIG. 16, it was imaged that oxygen metabolism is highest in the primary motor area using K and ORA. a is a phase difference image of K and can be said to be an oxygen exchange degree quantitative image. The darker the red, the larger the K-angle. K increased only in M1. On the other hand, the K angle decreased in the peripheral region. b is an intensity difference image of L, and the darker the red, the larger the L increase amount. Although L increased with the dumbbell weight increase, M1 was not the maximum increase site and was not limited to M1. The ORA of c increased only in M1, and increased oxygen exchange compared to the surrounding site. On the other hand, the surrounding region showed a negative change indicating passing.

(If breathing habits)
In order to distinguish between mouth breathing habit and nasal breathing habit using this brain habit diagnosis device S, using nasal tape or oral tape for 5-10 minutes, two different breathing methods are performed, The two groups of trajectories are statistically evaluated in two dimensions, and a determination is made based on the presence / absence of the difference and the phase frequency. The calculation unit 10 calculates the brain habit level described above, and the determination unit 12 determines based on the calculated brain habit level. For example, when low oxygen is shown, it shows that it becomes hard and it is determined that the habit is poor.

(Determination and instruction of breathing habits)
In breathing habits, two breathing methods, mouth breathing and nasal breathing are taken as tasks, and the measurement results are compared.
The criteria and instructions for mouth breathing and nasal breathing are as follows, for example.
1) If mouth breathing habits are strong, there is little burden on the brain and it is difficult to show hypoxia regulation Instruction: Instruct frequent nasal breathing exercises. At the same time, it gives instructions for checking nose biting.

2) When the habit of mouth breathing is low, the burden on the brain is strong and hypoxic regulation is shown.
Instruction: Continue this state

3) If both tasks are highly addictive, both do not show hypoxic regulation.
Instruction: Re-measure breathing status on another day

4) If both subjects are less addictive, both show hypoxic regulation.
Instructions: Instruct nasal breathing training from time to time. Re-measure breathing status on another day

  5) Check if there is a phase change in two measurements. Observe as a change in brain habits whether the phase changes from hypoxia to high oxygen or not. As strength decreases, addictiveness increases. If the phase changes clockwise, the addictiveness becomes stronger. If it changes counterclockwise, it becomes less addictive.

  In general, adverse effects of mouth breathing habits have been reported such as adverse effects on the oral cavity, malocclusion, gingivitis, maxillofacial growth and development, immune abnormalities, sleep disorders, restless ADHD, and emotional changes. .

FIG. 17A is a schematic explanatory diagram when diagnosing nasal breathing habits, and FIG. 17B is a schematic explanatory diagram when diagnosing mouth respiratory habits.
As shown in FIGS. 17A and 17B, the living body probe 1 is installed on the chest and head of the living body (subject), and simultaneously, respiratory measurement and brain measurement are performed.

  As shown in FIG. 17A, when diagnosing nasal breathing habits, the living body probe 1 is placed at the position of the BA 10 outside the brain, and the living body probe for measuring the number of breaths is placed on the chest. Then, the oral cavity is closed with a tape or the like.

  As shown in FIG. 17B, when diagnosing mouth breathing habits, the living body probe 1 is placed at the position of the BA 10 inside the brain, and the living body probe for measuring the number of breaths is placed on the chest. Then, the nasal cavity is closed with tape or the like.

  FIG. 18A is a graph showing fluctuations in blood volume during mouth breathing and time-series fluctuations in oxygen exchange volume, and FIG. 18B is a time series fluctuation in blood volume during nasal breathing and fluctuations in oxygen exchange volume. (C) is a graph showing the trajectory plotted on the CORE vector plane.

  As shown in FIGS. 18 (A) and 18 (B), there was no significant difference in the one-dimensional graph, but as shown in FIG. 18 (C), the two-dimensional graph clearly shows nasal breathing and mouth breathing. A significant difference was observed. That is, in mouth breathing, phase range 3 was shown at 30 seconds from the start of breathing, but not in nasal breathing. This means that oxygen is not used in nasal breathing. It shows that you can breathe with energy saving.

  FIG. 19 is a table showing the K angle and L value during nasal breathing and mouth breathing of four subjects.

  As shown in FIG. 19, for example, the average phase of nasal respiration of subject 6 is 159.3 and the average phase of mouth respiration is -160.6, nasal respiration indicates a low oxygen phase, and oral respiration indicates high oxygen. The phase difference between the two tasks is 159.3 − (− 160.6) = 329.9, which indicates that the subject is a mouth breathing habit.

  On the other hand, the average phase of nasal breathing of subject 1 is −30.5 and the average phase of mouth breathing is 102.9, which indicates a low oxygen phase in nasal breathing and high oxygen in mouth breathing. The phase difference between the two tasks is −30.5− (102.9) = − 133.4, which indicates that he is a person with nasal breathing habits.

  Here, if the subject 6 becomes the subject 1 after training, it can be seen that the breathing habits have improved by the phase of -30.5-159.3 = -189.8.

  FIG. 20 is a graph showing the phase change during nasal breathing and mouth breathing of two subjects.

  As shown in FIG. 20, subject A is shifted to phase 2 by mouth breathing, and subject B is shifted to phase -2 by mouth breathing. Accordingly, the determination unit 11 determines that the subject A has a strong nasal breathing habit and the subject B has rather a mouth breathing habit.

  FIG. 21 shows the trajectory on the two-dimensional diagram during mouth breathing and nasal breathing due to the difference in brain part, where (A) shows the case of the inner BA10 of the frontal lobe and (B) shows the case of the outer BA10 of the frontal lobe.

  As shown in FIG. 21, in mouth breathing, phase range 3 was shown at 30 seconds from the start of breathing but not in nasal breathing. Although the reaction was different in two places, mouth breathing in the two sites showed a phase range 2 and a phase range 3 and exhibited hypoxia.

  FIG. 22 shows the trajectory on the two-dimensional diagram of mouth breathing and nasal breathing by persons with different breathing habits. (A) is a nasal breather (average of 7 people), (B) is a strong mouth breathing habit In the case of (C), the case of mouth breathing habits is shown.

  As shown in FIG. 22, a significant difference was clearly recognized in two dimensions, and the phase range 3 was shown in the mouth breathing habit 30 seconds at the start of breathing, but it was not shown in the nose breathing habit. It shows that nasal breathing habits can breathe with energy saving without using oxygen.

(In the case of language usage habits)
If this brain habit diagnosis apparatus S is used, it is possible to determine the habit of language use from the brain habit. For example, prepare a pool of 100 words and ask the subject to hear one word at a time.
The calculation unit 10 calculates the K angle and the L value based on the index from the Vector-NIRS being listened to, so as to distinguish the hypoxic phase from the rest.

  Thus, the addictiveness of 100 pools is determined from the ratio of the number of hypoxia X / 100 = understanding words.

  The improvement support unit 12 gives an instruction to re-learn words that did not show hypoxia and improve habituation.

  FIG. 23A shows a trajectory on the two-dimensional diagram when the language understanding phase is shown, and FIG. 23B shows a trajectory on the two-dimensional diagram when the language non-understanding phase is shown.

  The determination unit 11 determines the degree of language usage habits from the frequencies indicating the language understanding phase and the non-understanding phase.

When phase range-2 is shown as shown in FIG. 23B, the determination unit 11 determines that the language usage habit is low, and the improvement support unit 12 issues an instruction to require strong training.
As shown in FIG. 23A, when the phase range 2 is shown, the determination unit 11 determines that the language usage frequency is high, and the improvement support unit 12 issues an instruction to learn a new word.

(For other brain habits)
The determination unit 11 can determine the degree of reading habits based on brain habits. For example, reading the same book for a predetermined time and comparing the data of people who have read the book with those who have not read the book, the reading habit level is determined for the book.

(About other parameters)
FIG. 24A is a graph showing the relationship between oxygen saturation Y and the horizontal axis representing the amount of oxidized Hb in ROI (O) and the vertical axis representing the amount of deoxidized Hb in ROI (D). Is a graph showing the relationship between the amount of change in oxidized Hb in ROI (ΔO), the amount of change in deoxidized Hb in ROI (ΔD) and the amount of change in Y angle (ΔY angle), (C) is OD It is a graph for demonstrating that the change of oxygen saturation is decided by the coordinate change ((DELTA) O and (DELTA) D) on a plane.
In FIG. 24 (A), assuming the upper tilt Y angle,
Oxygen saturation Y = 1 / (1 + tan (Y angle)) ··· Equation (3)
It becomes.

In FIG. 24B, the change in oxygen saturation Y depends on the ΔY angle.
Change in oxygen saturation ΔY = -1 / (1 + sin2 (ΔY angle)) Equation (4)
As obtained.
Thus, the change in oxygen saturation can occur if either the amount of change in oxidized Hb (ΔO) or the amount of change in deoxidized Hb (ΔD) changes.
Phase change amount = brain habit displacement degree is considered, and phase change can be calculated as ΔY.

(Difference between left and right brain in brain habit diagnosis)
In exercise habits, breathing habits, and other brain habits, the right and left brains are distinguished in consideration of the division of the role of the brain address.

  As for exercise habits, the right and left body use different brain addresses. The dominant hand is exactly an exercise habit, so in the exercise habit measurement of right hand tooth brushing and left hand tooth brushing, including the left and right primary motor areas (brain number 4), [activation of the right brain], [ The [right brain exercise habit value] and [left brain exercise habit value] related to exercise can be calculated from the frequency, intensity, and phase of activation of the left brain.

Furthermore, [Exercise habit right brain left brain difference value] = ([Right brain exercise habit value]-[Left brain exercise habit value]) / ([Right brain exercise habit value] + [Left brain exercise habit value]), ―1 ≦ [Exercise habit right brain Left brain difference] ≦ 1
From this, it can be determined whether the exercise habit depends on the right brain or the left brain.

  Many people are dominant in the left brain, but some people use the left and right brains or the right brain.

  Also, using hiragana is dominant in the left brain, and kanji is often used in both the left and right brains. Therefore, language habit measurement measures whether the right brain or the left brain has language habituation.

  In language manipulation, speech recall and conversation habits (speaking habits), including left and right brain addresses 44, 45, and 46, language manipulation is performed from the frequency, intensity, and phase of [right brain activation] and [left brain activation]. [Right brain language operation habit value] and [Left brain language operation habit value] can be calculated.

  Furthermore, [Language Manipulation Right Brain Left Value] = ([Right Language Maneuver Value]-[Left Brain Maneuver Value]) / ([Right Language Maneuver Value] + [Left Brain Maneuver Value]) ―1 ≦ [Language manipulation habit right brain left brain difference value] ≦ 1, it can be determined whether the language manipulation habit depends on the right brain or the left brain.

  Similarly, in language understanding habits, language understanding habits (word usage habits) and image understanding habits include right and left brain addresses 22, 39, and 40 [right brain activation], [left brain activation] [Right brain language comprehension habit value] and [Left brain language comprehension habit value] related to language understanding can be calculated from the frequency, intensity, and phase

  Furthermore, [Language understanding habit right brain left brain difference value] = ([Right brain language understanding habit value]-[Left brain language understanding habit value]) / ([Right brain language understanding habit value] + [Left brain language understanding habit value]) ―1 ≦ [Language comprehension habit right brain left brain difference value] ≦ 1, it can be determined whether the language comprehension habit and the image comprehension habit depend on the right brain or the left brain.

  If you have an understanding habit, or if you know what you have heard, make a mechanism that emits light and sound signals to make it easier to change instructions.

  In the visual habit (viewing habit), the left and right brain addresses 17, 18, 19, 7, 8, and 9 are measured simultaneously, and in the hearing habit (hearing habit), the left and right brain addresses 41, 42, 22, 21 are measured. In the sensory habits, the left and right brain addresses 3 are measured simultaneously. In the thinking habit (thinking habits), the left and right brain addresses 10, 8, and 9 are measured simultaneously. In social habits, left and right brain addresses 11, 10 can be measured simultaneously, and the left and right habits of each address can be compared.

  In breathing habits, right and left habituations can be compared for each address by simultaneously measuring left and right brain addresses 10, 44, 45, 46, 6, 4 and so on. However, when comparing breathing methods, the left and right brains may be compared, or only one side can be evaluated.

(program)
The program 13 according to the embodiment of the present invention shown in FIG. 1 is characterized by causing the control unit 7 of the brain habit diagnosis apparatus S to execute the above processing.

  The program 13 may be recorded on a recording medium such as a magnetic disk, a CD-ROM, or a semiconductor memory, or may be downloaded via a communication network.

  The present invention is not limited to the above-described embodiments, and various modifications can be made within the scope of the technical matters described in the claims.

For example, the types, combinations, and the like of various physiological indexes (change amounts, parameters) described in the specification and the drawings are examples, and are not limited thereto.
Further, depending on the brain habit to be diagnosed, one or more other parts that irradiate light other than the living body brain (the above-mentioned chest or other parts) and receive light from other parts of the living body brain The living body probe 1 may be further provided, and brain habits may be diagnosed by simultaneously measuring a living body brain and a part other than the brain.

  The present invention is used for diagnosing or instructing to improve a brain habit, which is a state of brain activity associated with a habit of a living body.

Claims (15)

Irradiating light to a predetermined part of the brain related to brain activity associated with living body habits, receiving light from the living body and receiving light information from the light receiving and receiving means, calculating, A brain habit diagnosis apparatus that diagnoses a brain habit that is a state of brain activity associated with the living body's habits, using a near-infrared spectroscopy, having a device body that performs control or storage,
The apparatus main body includes a concentration change amount Δ [HbO ₂ ] of oxidized hemoglobin, a concentration change amount Δ [Hb] of deoxidized hemoglobin, and a concentration change of oxidized hemoglobin based on light information from the light irradiation / receiving means. Calculation means for calculating various parameters derived based on the amount Δ [HbO ₂ ] and the deoxyhemoglobin Δ [Hb] concentration change amount, and the oxidized hemoglobin concentration change Δ [HbO calculated by the calculation means ₂ ], based on the amount of change Δ [Hb] of deoxygenated hemoglobin and various parameters, a determination means for determining the presence / absence and degree of brain habits, and the presence / absence and degree of brain habits determined by the determination means are displayed. possess and display means,
The calculation means calculates the concentration change ΔCBV of the blood volume in the predetermined part of the brain by the equation (1),
The concentration change amount ΔCOE of the oxygen exchange amount in a predetermined part of the brain is calculated by Equation (2),
In the two-dimensional diagram showing the relationship between the concentration change ΔCBV of the blood volume and the concentration change ΔCOE of the oxygen exchange amount, a scalar of the trajectory vector from the starting point of the measurement obtained by plotting the problem over time Calculate L value and phase K angle,
The determination means determines the presence or absence or degree of brain habit based on the scalar L value and the phase K angle,
ΔCBV = Δ [Hb] + Δ [HbO ₂ ]... Formula (1)
ΔCOE = Δ [Hb] −Δ [HbO ₂ ]... Formula (2)
The calculating means uses a first axis as the blood volume concentration change amount ΔCBV, a second axis orthogonal to the first axis as the oxygen exchange amount concentration change amount ΔCOE, and the first axis as the first axis. On the other hand, the third axis set to −45 degrees is defined as the concentration change Δ [HbO ₂ ] of the oxidized hemoglobin, and the fourth axis set to +45 degrees with respect to the first axis Dividing into eight phase regions on the two-dimensional diagram with the concentration change Δ [Hb] of oxidized hemoglobin, and calculating the frequency for each phase region where the trajectory obtained by plotting the task over time is present And
The determination means determines the presence or absence and degree of brain habits according to the frequency for each phase region,
A brain habit diagnosis apparatus characterized by that. The brain habit diagnosis apparatus according to claim 1 , wherein the calculating unit calculates a phase distribution rate, a highest frequency phase, and an average phase distribution rate for each phase region. The calculating means, brain habit diagnostic apparatus according to claim 1 or 2, characterized in that to calculate the frequency and the ratio of the phase zone for each of a plurality of issues. The calculation means calculates the oxygen saturation Y on the OD plane from a two-dimensional diagram in which the horizontal axis represents the oxidized Hb amount in the region of interest ROI (O) and the vertical axis represents the deoxidized Hb amount in the ROI (D). Is calculated by the equation (3) as the inclination Y angle of
The brain habit diagnosis apparatus according to any one of claims 1 to 3 , wherein
Oxygen saturation Y = 1 / (1 + tan (Y angle)) ··· Equation (3) The calculation means uses the amount of oxidized Hb in the region of interest ROI (O) on the horizontal axis and the amount of deoxidized Hb in ROI (D) on the vertical axis, and the amount of change in oxidized Hb (ΔO) and deoxidation. From the two-dimensional diagram showing the relationship of the amount of change (ΔD) of the type Hb, the change ΔY of the oxygen saturation Y is calculated by the equation (4) using the change ΔY angle of the inclination Y angle on the OD plane.
The brain habit diagnosis apparatus according to any one of claims 1 to 4 , wherein the brain habit diagnosis apparatus is characterized.
Change in oxygen saturation ΔY = -1 / (1 + sin2 (ΔY angle)) Equation (4) The calculation means calculates an oxygen adjustment area ORA from the scalar L value of the orbital vector from the origin of measurement start and the phase K angle by the equation (5).
The brain habit diagnosis apparatus according to any one of claims 1 to 5 , wherein the brain habit diagnosis apparatus is described above.

The determination means determines whether there is a mouth breathing habit or a nasal breathing habit by comparing the case where the oral cavity of the living body is blocked with the case where the nasal cavity is blocked,
The brain habit diagnosis apparatus according to any one of claims 1 to 6 , The determination means determines the habit of using words by letting the living body speak words,
The brain habit diagnosis apparatus according to any one of claims 1 to 6 , The determination means determines exercise habits by giving an object to the living body.
The brain habit diagnosis apparatus according to any one of claims 1 to 6 , wherein the brain habit diagnosis apparatus is characterized. And further comprising improvement support means for diagnosing the living body at different times and determining whether or not the brain habit is improved based on a change in a result of determination of the brain habit by the determination means. The brain habit diagnosis apparatus according to any one of 1 to 9 . The improvement support means determines instructions for making good brain habits based on the determination result of the brain habits by the determining means,
The display means displays the content of the instruction;
The brain habit diagnosis apparatus according to claim 10 . The light irradiation light receiving means is installed at a site corresponding to the address for brain division related to the brain habit to be diagnosed,
The brain habit diagnosis apparatus according to any one of claims 1 to 11 , wherein the brain habit diagnosis apparatus is characterized. The system further includes one or more second light irradiation / light receiving means for irradiating light on a part other than the brain of the living body and receiving light from a part other than the brain of the living body, Diagnosing brain habits by simultaneously measuring
The brain habit diagnosis apparatus according to any one of claims 1 to 12 , wherein the brain habit diagnosis apparatus is described above. The brain habit diagnosis apparatus according to claim 13 , wherein the part other than the brain is a breast of the living body. The program which makes the process of the apparatus main body of the brain habit diagnosis apparatus as described in any one of the said Claims 1 thru | or 14 be performed.

Images