WO2019136745A1 - Magnetic resonance image based brain age test method and apparatus - Google Patents

Magnetic resonance image based brain age test method and apparatus Download PDF

Info

Publication number
WO2019136745A1
WO2019136745A1 PCT/CN2018/072603 CN2018072603W WO2019136745A1 WO 2019136745 A1 WO2019136745 A1 WO 2019136745A1 CN 2018072603 W CN2018072603 W CN 2018072603W WO 2019136745 A1 WO2019136745 A1 WO 2019136745A1
Authority
WO
WIPO (PCT)
Prior art keywords
brain
magnetic resonance
age
resonance image
weighted
Prior art date
Application number
PCT/CN2018/072603
Other languages
French (fr)
Chinese (zh)
Inventor
罗怡珊
Original Assignee
深圳博脑医疗科技有限公司
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 深圳博脑医疗科技有限公司 filed Critical 深圳博脑医疗科技有限公司
Priority to PCT/CN2018/072603 priority Critical patent/WO2019136745A1/en
Priority to CN201880000008.4A priority patent/CN110337670B/en
Publication of WO2019136745A1 publication Critical patent/WO2019136745A1/en

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/143Segmentation; Edge detection involving probabilistic approaches, e.g. Markov random field [MRF] modelling

Definitions

  • the present application belongs to the field of image processing technologies, and in particular, to a brain age testing method based on magnetic resonance images, a brain age testing device based on magnetic resonance images, an electronic device, and a computer readable storage medium.
  • the brain age (referred to as brain age) is used to measure the development of the human brain.
  • people also have a physiological age, which is used to describe the time from the date of birth. Normally, the physiological age is not equal to the brain age. For example, when a person's brain develops a lesion, it causes a sharp aging of the brain, a memory loss and a slow response, which makes the brain older than the physiological age. When people exercise regularly and maintain physical and mental pleasure for a long time, it will delay the aging of the brain. That makes the brain age less than the physiological age.
  • the determination of brain age can raise people's awareness of brain health, intervene in brain health in advance, and delay brain aging.
  • the embodiments of the present application provide a brain age testing method based on magnetic resonance images, a brain age testing device based on magnetic resonance images, an electronic device, and a computer readable storage medium, which can implement testing of human brain age.
  • a first aspect of the present application provides a brain age testing method based on a magnetic resonance image, comprising:
  • the brain age estimation model is based on a normalized volume of the training sample individual and a brain atrophy value And brain age training is available.
  • a second aspect of the present application provides a brain age testing device based on a magnetic resonance image, comprising:
  • An image acquisition unit configured to acquire a T1-weighted brain magnetic resonance image of the individual to be tested
  • An image analyzing unit configured to determine a brain structure, a brain lobe, and a brain tissue segmentation map of the test subject based on the T1-weighted brain magnetic resonance image
  • a parameter calculation unit configured to calculate a normalized volume of the brain structure, and calculate a brain atrophy value of the brain lobe based on the brain tissue segmentation map;
  • a brain age acquisition unit configured to input the normalized volume and the brain atrophy value into a brain age estimation model, and obtain a brain age of the test subject, wherein the brain age estimation model is based on the individuality of the training sample individual A volume, brain atrophy value and brain age training are obtained.
  • a third aspect of the present application provides an electronic device including a memory, a processor, and a computer program stored in the memory and operable on the processor, wherein the processor executes the computer program The steps of the method as described above are achieved.
  • a fourth aspect of the present application provides a computer readable storage medium storing a computer program, wherein the computer program is executed by a processor to implement the steps of the method as described above.
  • the present application provides a brain age testing method based on magnetic resonance images, first obtaining a T1-weighted brain magnetic resonance image of an individual to be tested, and secondly, analyzing the T1-weighted brain magnetic resonance image to obtain a brain structure, Brain lobe and brain tissue segmentation map, and calculating a normalized volume of the test subject and a brain atrophy value according to the obtained brain structure, brain lobe and brain tissue segmentation map, the normalized volume and the brain
  • the atrophy value is input into the brain age estimation model, and the brain age of the individual to be tested is obtained, thereby realizing the test of the brain age of the individual to be tested.
  • the brain size and brain atrophy value are used as brain age test parameters.
  • the user can estimate his or her brain age, so as to understand the state of his brain health, and facilitate the intervention of his brain health in advance to delay.
  • Brain aging through the technical solutions provided by this application can estimate brain degeneration and improve people's awareness of brain health.
  • FIG. 1 is a schematic flowchart showing an implementation process of a brain age testing method based on magnetic resonance images according to Embodiment 1 of the present application;
  • FIG. 2 is a schematic diagram showing an implementation flow of a brain age testing method based on magnetic resonance images provided in Embodiment 2 of the present application;
  • FIG. 3 is a schematic structural diagram of a brain age testing device based on magnetic resonance images provided in Embodiment 3 of the present application;
  • FIG. 4 is a schematic structural diagram of a brain age testing device based on magnetic resonance images provided in Embodiment 4 of the present application;
  • FIG. 5 is a schematic diagram of an electronic device according to Embodiment 5 of the present application.
  • the brain age testing method based on the magnetic resonance image provided by the embodiment of the present application is applicable to an electronic device.
  • the above electronic device includes, but is not limited to, a desktop computer, a tablet computer, a cloud server, a mobile phone terminal, and the like.
  • the brain age testing method in the embodiment of the present application includes:
  • Step S101 acquiring a T1-weighted brain magnetic resonance image of the individual to be tested
  • the T1 weighted brain magnetic resonance image of the individual to be tested is first acquired, so that the related parameters about brain development of the test subject can be obtained from the T1 weighted brain magnetic resonance image, thereby utilizing the The parameter estimates the brain age of the individual to be tested, wherein the T1-weighted brain magnetic resonance image of the individual to be tested is a three-dimensional image.
  • Step S102 determining a brain structure, a brain lobe, and a brain tissue segmentation map of the test subject according to the T1 weighted brain magnetic resonance image;
  • the brain tissue segmentation map includes a gray matter segmentation map, a white matter segmentation map, and a cerebrospinal fluid segmentation map.
  • the above T1-weighted brain magnetic resonance image can be registered with a preset brain template library.
  • the brain template library is pre-stored, and the brain template library includes T1 weighted brain magnetic resonance images of two or more different brains, that is, including two or more different templates, preferably, to ensure the brain structure of the individual to be tested,
  • the accuracy of the brain leaf and brain tissue segmentation map which may include T1 weighted brain magnetic resonance images of people with different brain health, gender, age and age between 10 and 90 years old, and the number of images is more than 5 One.
  • the brain structure and the brain leaf can be manually segmented in each T1 weighted brain magnetic resonance image (each template) in the brain template library, and the brain structure and brain corresponding to each template in the brain template library are obtained.
  • Leaf, the above brain structure may include brain parenchyma, cerebellum, hippocampus, amygdala, ventricle, lateral ventricle, thalamus, caudate nucleus, putamen, globus pallidus, nucleus accumbens, midbrain, pons, cerebral ventricle, etc.
  • the brain tissue is automatically segmented in each T1 weighted brain magnetic resonance image in the brain template library, and manual manual correction is performed on the basis of automatic computer segmentation, so that the brain tissue probability map corresponding to each template in the brain template library is obtained.
  • the tissue probability map includes a white matter probability map, a gray matter probability map, and a cerebrospinal fluid probability map.
  • the brain structure, the brain lobe, and the brain tissue segmentation map of the individual to be tested based on the brain template library may be:
  • the above nonlinear registration can adopt a symmetric nonlinear registration algorithm based on a differential homeomorphic model.
  • the brain tissue probability map of the above-mentioned test subject is: the brain tissue probability map of the template with the highest similarity to the T1-weighted brain magnetic resonance image of the individual to be tested is obtained by spatial mapping.
  • the label fusion algorithm can be used to fuse the brain structure partitions of the individuals to be tested according to each template to obtain the brain structure partition of the final test subject. Similarly, using the label fusion algorithm, the test samples obtained according to each template are to be tested. The individual's brain lobe is fused to obtain the brain lobe of the final individual to be tested.
  • the label fusion algorithm is an algorithm that combines the segmentation results corresponding to each template by using a desired maximum algorithm. The label fusion algorithm can solve the deviation caused by a single template, and the result is more accurate.
  • Commonly used tag fusion algorithms include a reliability-based tag fusion algorithm, a majority consent rule tag fusion algorithm, and a weighted tag fusion algorithm.
  • Step S103 calculating a normalized volume of the brain structure, and calculating a brain atrophy value of the brain lobe based on the brain tissue segmentation map;
  • the volume of the brain structure of the individual to be tested obtained in step S102 can be estimated, and the total volume of the brain is divided to calculate the normalized volume of the brain structure of the individual to be tested.
  • the brain gray matter volume, the white matter volume, and the cerebrospinal fluid volume in the brain lobe may be calculated according to the brain lobe partition and the brain tissue segmentation map obtained in step S102, and the brain atrophy value of the brain lobe is calculated according to the brain atrophy calculation formula, wherein The formula for brain atrophy is:
  • the brain atrophy value obtained in the above formula can directly reflect the degree of atrophy of the brain lobe.
  • Step S104 input the normalized volume and the brain atrophy value into a brain age estimation model, and obtain the brain age of the test subject;
  • the brain age estimation model may be established in advance, firstly, the normalized volume of the training sample individual and the brain atrophy value are collected, and the normalized volume of the training sample individual and the brain atrophy value are used as independent variables, and the training will be performed.
  • the brain age of the sample individuals was used as the dependent variable, and a linear support vector machine was used to establish a brain age estimation model.
  • the establishment of the brain age estimation model is not limited to the linear support vector machine, but also the hidden Markov model, the neural network, etc. can be used to establish the brain age estimation model.
  • the optimization goal is: The specific values of W and b when the optimization target value is minimum can be obtained.
  • the data fitting algorithm such as gradient descent algorithm or genetic algorithm can be used to solve the values of W and b, so as to obtain the brain age estimation model.
  • the optimization target can also be The formula for optimizing the target is not limited here.
  • the T1-weighted brain magnetic resonance image may also be pre-processed after acquiring the T1-weighted brain magnetic resonance image of the individual to be tested, the pre-processing including the following Or more than two: noise reduction, de-fielding, and pixel range normalization.
  • the present invention estimates the Rice noise variance based on the gray-scale distribution skewness noise estimation method, and then uses the non-local mean algorithm to reduce noise according to the estimated noise variance; the present application de-biasing field is used to remove the uneven magnetic field.
  • the phenomenon that the gray level of the brain tissue is inconsistent; this application uses the method of histogram matching to normalize the intensity ranges of different images to a common range.
  • the T1 weighted brain magnetic resonance images of the training samples were from different models of 16 different hospitals in China, including SIEMENS. GE, PHILIPS NMR machine. Firstly, the T1 weighted brain magnetic resonance images of the collected training sample individuals are preprocessed; the 38 normalized volumes corresponding to each training sample individual and 12 brain atrophy values are obtained; since the healthy individual's brain age and physiological age are close Equal, so when training the brain age estimation model, the physiological age of the training sample individual is used as the brain age of the training sample individual to establish a brain age estimation model.
  • the difference between the established brain age estimation model and its physiological age is 5.44 years old, that is, according to the technical solution provided by the first embodiment of the present application, the error of the established brain age estimation model is 5.44 years old.
  • the individuals to be tested were 12 healthy individuals (age range 67.3 ⁇ 9.3 years) and 14 individuals with Alzheimer's disease (age range 61.3 ⁇ 15.6 years).
  • the T1-weighted brain magnetic resonance images of the individuals to be tested were from PHILIPS NMR.
  • the resonance machine first preprocesses the T1-weighted brain magnetic resonance image of the individual to be measured, and extracts 38 normalized brain volumes and 12 brain atrophy values of the individuals to be tested, and estimates according to the brain age estimation model obtained above. The age of the individual to be tested.
  • the test results showed that the average age difference between the predicted brain age and the physiological age of the above 12 healthy individuals was 6.4 years, and the predicted brain age and the physiological age difference of the individuals with Alzheimer's disease were 19.7 years old. It can be seen from the above test results that the brain of the Alzheimer's disease patient has a higher degree of brain aging than the healthy individual.
  • the brain structure size and the brain atrophy value are used as the brain age test parameters for the first time.
  • the user can estimate the brain age of the user to understand the state of the brain health and facilitate the advancement of the brain health.
  • Intervention to delay brain aging, through the technical solutions provided in this application can estimate brain degeneration and improve people's awareness of brain health.
  • the brain age testing method in the embodiment of the present application includes:
  • Step S201 acquiring a T1-weighted brain magnetic resonance image of the training sample individual
  • a method for establishing a brain age estimation model is specifically given.
  • a large number of parameters related to brain development are often required.
  • many parameters related to brain development are not necessary, even if this parameter is added to the brain age estimation model, it cannot be To increase the accuracy of the brain age estimation model too much, it will occupy a large amount of computing resources in the process of calculating the brain age of the individual to be tested. Therefore, this embodiment provides a method for establishing a brain age estimation model, which can eliminate parameters that are not needed when establishing a brain age estimation model, and release certain computing resources.
  • the T1-weighted brain magnetic resonance image of each training sample individual may be pre-processed, and the above pre-processing includes one or more of the following: noise reduction, de-biasing
  • the field and pixel range are normalized.
  • Step S202 acquiring a brain lobe and a brain tissue segmentation map of each training sample individual, and calculating a brain atrophy value of the brain lobe based on the brain tissue segmentation map;
  • the brain lobe and the brain tissue probability map of each training sample individual are obtained, and the brain atrophy value of the brain lobe is calculated, and the specific implementation method can be referred to
  • the first embodiment S102 and S103 are as described above, and are not described herein again.
  • brain atrophy values corresponding to each training sample individual including left frontal lobe, left parietal lobe, left occipital lobe, left temporal lobe, left cingulate gyrus and left insular lobe, and right frontal lobe.
  • Step S203 acquiring an i-th brain structure of each training sample individual, calculating a normalized volume of the i-th brain structure, a normalized volume based on the i-th brain structure, and brain atrophy of the brain of the training sample individual Value, establish an i-th brain age estimation model, and obtain the error of the above i-th brain age estimation model;
  • the brain structure in the brain is very large, but the size of some brain structures changes little with the development of the brain, and some changes in the size of the brain structure do not affect the degree of aging of the brain.
  • the normalized volume is not required in brain age calculations. Therefore, in the embodiment of the present application, we can first select a certain brain structure, calculate the normalized volume of the selected brain structure, and use the normalized volume and the above 12 brain leaf atrophy values as the model for establishing the brain age.
  • the 13 dependent variables can be used to establish the brain age estimation model corresponding to the brain structure by using the linear support vector machine, and the minimum value of the optimization target is used as the error of the brain age estimation model.
  • Step S204 determining whether to traverse all brain structures of each training sample individual
  • step S203 it is determined whether all brain structures have been traversed. If not, step S205 is performed, and if yes, step S206 is performed.
  • Step S205 increasing the value of i by one
  • the normalized volume of the next brain structure and the above 12 brain leaf atrophy values are used as the 13 dependent variables for establishing the brain age estimation model, and the next brain structure is established.
  • the corresponding brain age estimation model is obtained, and the error of the brain age estimation model is obtained.
  • Step S206 selecting a normalized volume of the brain structure corresponding to the brain age estimation model with a small error, and establishing a final brain age estimation model according to the normalized volume of the selected brain structure and the brain blade atrophy value;
  • the normalized volume of the brain structure corresponding to the brain age estimation model with less error is selected.
  • the threshold may be preset, for example, 10 years old, and the error of the brain age estimation model corresponding to each brain structure obtained above is compared with 10 years old, and if it is less than the preset 10 years old, the normalized volume of the brain structure is selected.
  • the normalized volume of the brain structure that may be selected is only the normalized volume of the cerebellum and the normalized volume of the hippocampus.
  • the normalized volume of the selected brain structure ie, the normalized volume of the cerebellum and the hippocampus normalized.
  • the volume, and 12 brain leaf atrophy values were used as the 14 dependent variables of the final brain age estimation model to establish the final brain age estimation model.
  • the determined brain structure of the test subject is corresponding to the brain structure in the final brain age estimation model, for example, if the final brain age estimation model has only the cerebellar normalized volume and the hippocampus If the volume is normalized, then in step S208, only the cerebellum and hippocampus of the individual to be tested need to be determined. In step S209, only the normalized volume of the cerebellum and the normalized volume of the brain structure need to be calculated.
  • the foregoing steps S207-S209 are the same as the implementations of the steps S101-S103 in the first embodiment. For details, refer to the description of the first embodiment, and details are not described herein again.
  • not every brain age test needs to perform steps S201-S206. After the final brain age estimation model is established, steps S201-S206 need not be performed in subsequent brain age test.
  • a method for establishing a brain age estimation model is specifically provided, which can eliminate parameters that are not needed when establishing a brain age estimation model, release certain computing resources, and firstly determine brain structure size and brain atrophy value.
  • the user can estimate the age of the brain by using the solution of the present application, in order to understand the state of the brain health, and facilitate the intervention of the brain health in advance to delay the aging of the brain, through the technical solution provided by the application. Brain degeneration can be estimated and people's awareness of brain health can be improved.
  • a third embodiment of the present invention provides a brain age testing device based on a magnetic resonance image. For convenience of description, only parts related to the present application are shown. As shown in FIG. 3, the brain age testing device 300 includes:
  • the image obtaining unit 301 is configured to acquire a T1-weighted brain magnetic resonance image of the individual to be tested;
  • the image analyzing unit 302 is configured to determine a brain structure, a brain lobe, and a brain tissue segmentation map of the test subject according to the T1 weighted brain magnetic resonance image, where the brain tissue segmentation map includes a gray matter segmentation map, a white matter segmentation map, and a cerebrospinal fluid Split the map;
  • the parameter calculation unit 303 is configured to calculate a normalized volume of the brain structure, and calculate a brain atrophy value of the brain lobe based on the brain tissue segmentation map;
  • the brain age obtaining unit 304 is configured to input the normalized volume and the brain atrophy value into a brain age estimation model to obtain a brain age of the test subject, wherein the brain age estimation model is based on a normalized volume of the training sample individual Brain atrophy and brain age training are obtained.
  • the image analysis unit 302 is specifically configured to:
  • the T1 weighted brain magnetic resonance image is registered with a preset brain template library to obtain a brain structure, a brain lobe, and a brain tissue probability map corresponding to the T1 weighted brain magnetic resonance image; and the T1 weighted brain magnetic resonance image is obtained.
  • Corresponding brain tissue probability map is used as a priori knowledge, and brain tissue segmentation is performed on the T1-weighted brain magnetic resonance image to obtain a brain tissue segmentation map corresponding to the T1-weighted brain magnetic resonance image;
  • the brain template library includes: a T1-weighted brain magnetic resonance image of two or more different brains, and a brain structure, a brain lobe, and a brain tissue probability map corresponding to each T1-weighted brain magnetic resonance image, and the brain tissue probability map. Including the gray matter probability map, the white matter probability map and the cerebrospinal fluid probability map.
  • the parameter calculation unit 303 is specifically configured to:
  • the brain atrophy value of the brain lobe is calculated according to the brain atrophy calculation formula, wherein the above brain atrophy calculation formula is:
  • the brain age testing device 300 further includes:
  • the model building unit is configured to use the normalized volume of the training sample individual and the brain atrophy value as independent variables, and use the brain support age of the training individual as a dependent variable, and apply a linear support vector machine to establish a brain age estimation model.
  • the brain age testing device 300 further includes:
  • the individual preprocessing unit to be tested is configured to preprocess the T1-weighted brain magnetic resonance image of the test subject to be tested before determining the brain structure, the brain lobe and the brain tissue segmentation map of the test subject according to the T1 weighted brain magnetic resonance image.
  • the above preprocessing includes one or more of the following: noise reduction, de-biasing, and pixel range normalization.
  • a fourth embodiment of the present invention provides a brain age testing device based on a magnetic resonance image. For convenience of explanation, only parts related to the present application are shown. As shown in FIG. 4, the brain age testing device 400 includes:
  • a training sample individual image acquiring unit 401 configured to acquire a T1-weighted brain magnetic resonance image of the training sample individual
  • the training sample individual brain atrophy calculation unit 402 is configured to acquire a brain lobe and a brain tissue segmentation map of each training sample individual, and calculate a brain atrophy value of the brain lobe based on the brain tissue segmentation map;
  • the i-th brain age model acquisition unit 403 is configured to acquire an i-th brain structure of each training sample individual, calculate a normalized volume of the i-th brain structure, based on a normalized volume of the i-th brain structure, and Brain atrophy value of the brain, establish an i-th brain age estimation model, and obtain the error of the above i-th brain age estimation model;
  • the determining unit 404 is configured to determine whether to traverse all brain structures of each training sample individual;
  • the final model establishing unit 406 is configured to select a normalized volume of the brain structure corresponding to the brain age estimation model with less error, and establish a final brain according to the normalized volume of the selected brain structure and the brain blade atrophy value.
  • Age estimation model
  • the individual image acquisition unit 407 to be tested is configured to acquire a T1-weighted brain magnetic resonance image of the individual to be tested;
  • the individual image analysis unit 408 to be tested is configured to determine a brain structure, a brain lobe, and a brain tissue segmentation map of the individual to be tested based on the T1 weighted brain magnetic resonance image of the individual to be tested;
  • the individual parameter calculation unit 409 to be tested is configured to calculate a normalized volume of the brain structure of the individual to be tested, and calculate a brain atrophy value of the brain of the individual to be tested based on the brain tissue segmentation map of the individual to be tested;
  • the individual brain age obtaining unit 410 is configured to input the normalized volume of the test subject and the brain atrophy value of the test subject into a final brain age estimation model, and obtain the brain age of the test subject;
  • the determined brain structure of the test subject is corresponding to the brain structure in the final brain age estimation model, for example, if the final brain age estimation model has only the cerebellar normalized volume and the hippocampus
  • the individual image analysis unit 408 to be tested only needs to determine the cerebellum and hippocampus of the individual to be tested
  • the individual parameter calculation unit 409 to be tested only needs to calculate the normalized volume of the cerebellum and the normalization of the brain structure. volume.
  • the foregoing units 407-409 are the same as the implementations of the units 301-303 in the third embodiment. For details, refer to the description of the third embodiment, and details are not described herein again.
  • the brain age testing device 400 further includes:
  • a pre-processing unit configured to pre-process the T1-weighted brain magnetic resonance image of the individual to be tested, and to obtain each training sample individual, before determining the brain structure, the brain lobe, and the brain tissue segmentation map of the individual to be tested Before the brain structure, brain lobe and brain tissue segmentation map, the T1-weighted brain magnetic resonance images of each training sample individual are preprocessed, and the above preprocessing includes one or more of the following: noise reduction, de-biasing field, pixel range Normalized.
  • FIG. 5 is a schematic diagram of an electronic device according to Embodiment 5 of the present application.
  • the electronic device 5 of this embodiment includes a processor 50, a memory 51, and a computer program 52 stored in the above-described memory 51 and operable on the processor 50 described above.
  • the processor 50 executes the computer program 52 described above, the steps in the above various method embodiments are implemented, such as steps S101 to S104 shown in FIG.
  • the processor 50 executes the computer program 52
  • the functions of the modules/units in the above-described respective device embodiments such as the functions of the modules 301 to 304 shown in FIG. 3, are implemented.
  • the above electronic device 5 may be a computing device such as a desktop computer, a notebook, a palmtop computer, and a cloud server.
  • the above electronic device may include, but is not limited to, the processor 50 and the memory 51. It will be understood by those skilled in the art that FIG. 5 is only an example of the electronic device 5, and does not constitute a limitation on the electronic device 5, and may include more or less components than those illustrated, or combine some components, or different components.
  • the electronic device 5 described above may further include an input/output device, a network access device, a bus, and the like.
  • the processor 50 may be a central processing unit (CPU), or may be other general-purpose processors, a digital signal processor (DSP), an application specific integrated circuit (ASIC), Field-Programmable Gate Array (FPGA) or other programmable logic device, discrete gate or transistor logic device, discrete hardware components, etc.
  • the general purpose processor may be a microprocessor or the processor or any conventional processor or the like.
  • the memory 51 may be an internal storage unit of the electronic device 5, such as a hard disk or an internal memory of the electronic device 5.
  • the memory 51 may be an external storage device of the electronic device 5, such as a plug-in hard disk equipped with the above-mentioned electronic device 5, a smart memory card (SMC), a Secure Digital (SD) card, and a flash memory. Flash card, etc.
  • the above-mentioned memory 51 may also include both an internal storage unit of the above-described electronic device 5 and an external storage device.
  • the above memory 51 is used to store the above computer program and other programs and data required for the above electronic device.
  • the above-described memory 51 can also be used to temporarily store data that has been output or is about to be output.
  • each functional unit and module in the foregoing system may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit, and the integrated unit may be hardware.
  • Formal implementation can also be implemented in the form of software functional units.
  • the specific names of the respective functional units and modules are only for the purpose of facilitating mutual differentiation, and are not intended to limit the scope of protection of the present application.
  • the disclosed apparatus/electronic device and method may be implemented in other manners.
  • the device/electronic device embodiments described above are merely illustrative.
  • the division of the above modules or units is only a logical functional division, and may be implemented in another manner, such as multiple units or Components can be combined or integrated into another system, or some features can be ignored or not executed.
  • the mutual coupling or direct coupling or communication connection shown or discussed may be an indirect coupling or communication connection through some interface, device or unit, and may be in electrical, mechanical or other form.
  • the units described above as separate components may or may not be physically separated.
  • the components displayed as units may or may not be physical units, that is, may be located in one place, or may be distributed to multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of the embodiment.
  • each functional unit in each embodiment of the present application may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit.
  • the above integrated unit can be implemented in the form of hardware or in the form of a software functional unit.
  • the above-described integrated modules/units if implemented in the form of software functional units and sold or used as stand-alone products, may be stored in a computer readable storage medium. Based on such understanding, the present application implements all or part of the processes in the above embodiments, and may also be completed by a computer program to instruct related hardware.
  • the computer program may be stored in a computer readable storage medium.
  • the steps of the various method embodiments described above may be implemented when executed by a processor.
  • the above computer program comprises computer program code
  • the computer program code may be in the form of source code, object code form, executable file or some intermediate form.
  • the computer readable medium may include any entity or device capable of carrying the above computer program code, a recording medium, a USB flash drive, a removable hard disk, a magnetic disk, an optical disk, a computer memory, a Read-Only Memory (ROM), a random Access memory (RAM, Random Access Memory), electrical carrier signals, telecommunications signals, and software distribution media. It should be noted that the contents of the above computer readable medium may be appropriately increased or decreased according to the requirements of legislation and patent practice in jurisdictions. For example, in some jurisdictions, according to legislation and patent practice, computer readable media are not Includes electrical carrier signals and telecommunications signals.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Probability & Statistics with Applications (AREA)
  • Software Systems (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

Disclosed are a magnetic resonance image based brain age test method, a magnetic resonance image based brain age test apparatus, an electronic device and a computer-readable storage medium. The magnetic resonance image based brain age test method comprises: acquiring a T1-weighted brain magnetic resonance image of an object to be tested (S101); determining, based on the T1-weighted brain magnetic resonance image, a brain structure, a brain lobe and a brain tissue segmentation map of the object to be tested (S102); calculating a normalized volume of the brain structure and a brain atrophy value of the brain lobe (S103); and inputting the normalized volume and the brain atrophy value into a brain age estimation model to acquire a brain age of the object to be tested (S104). The method can be used to estimate a brain age, and facilitates people in intervening with the health of the brain in advance and improves the health awareness of people for the brain.

Description

一种基于磁共振图像的脑龄测试方法及脑龄测试装置Brain age test method based on magnetic resonance image and brain age test device 技术领域Technical field
本申请属于图像处理技术领域,尤其涉及一种基于磁共振图像的脑龄测试方法、基于磁共振图像的脑龄测试装置、电子设备及计算机可读存储介质。The present application belongs to the field of image processing technologies, and in particular, to a brain age testing method based on magnetic resonance images, a brain age testing device based on magnetic resonance images, an electronic device, and a computer readable storage medium.
背景技术Background technique
人从胎儿到婴幼儿,再到儿童、青少年、中年、老年,大脑和身体同时在发育,常用大脑年龄(简称脑龄)衡量人大脑的发育程度。此外,人还具有生理年龄,用来描述人距离出生日期的时间。通常情况下,生理年龄与大脑年龄并不相等。比如,当人的大脑发生病变时,会导致大脑急剧衰老,出现记忆力衰退以及反应迟钝等现象,使得大脑年龄大于生理年龄;当人经常锻炼身体,长时间保持身心愉悦时,会延缓大脑的衰老,使得大脑年龄小于生理年龄。From fetus to infant, to children, adolescents, middle-aged, old age, brain and body at the same time, the brain age (referred to as brain age) is used to measure the development of the human brain. In addition, people also have a physiological age, which is used to describe the time from the date of birth. Normally, the physiological age is not equal to the brain age. For example, when a person's brain develops a lesion, it causes a sharp aging of the brain, a memory loss and a slow response, which makes the brain older than the physiological age. When people exercise regularly and maintain physical and mental pleasure for a long time, it will delay the aging of the brain. That makes the brain age less than the physiological age.
大脑年龄的确定,可提高人们对大脑健康的意识,提前对大脑健康进行干预,延缓大脑衰老。The determination of brain age can raise people's awareness of brain health, intervene in brain health in advance, and delay brain aging.
技术问题technical problem
有鉴于此,本申请实施例提供了一种基于磁共振图像的脑龄测试方法、基于磁共振图像的脑龄测试装置、电子设备及计算机可读存储介质,可以实现对人体脑龄的测试。In view of this, the embodiments of the present application provide a brain age testing method based on magnetic resonance images, a brain age testing device based on magnetic resonance images, an electronic device, and a computer readable storage medium, which can implement testing of human brain age.
技术解决方案Technical solution
本申请第一方面提供了一种基于磁共振图像的脑龄测试方法,包括:A first aspect of the present application provides a brain age testing method based on a magnetic resonance image, comprising:
获取待测个体的T1加权大脑磁共振图像;Obtaining a T1-weighted brain magnetic resonance image of the individual to be tested;
基于所述T1加权大脑磁共振图像确定所述待测个体的脑结构、脑叶以及脑组织分割图谱;Determining a brain structure, a brain lobe, and a brain tissue segmentation map of the test subject based on the T1-weighted brain magnetic resonance image;
计算所述脑结构的归一化体积,并基于所述脑组织分割图谱计算所述脑叶的脑萎缩值;Calculating a normalized volume of the brain structure, and calculating a brain atrophy value of the brain lobe based on the brain tissue segmentation map;
将所述归一化体积以及所述脑萎缩值输入脑龄估计模型,获取所述待测个体的脑龄,其中,所述脑龄估计模型根据训练样本个体的归一化体积、脑萎缩值以及脑龄训练得到。Entering the normalized volume and the brain atrophy value into a brain age estimation model to obtain a brain age of the test subject, wherein the brain age estimation model is based on a normalized volume of the training sample individual and a brain atrophy value And brain age training is available.
本申请第二方面提供了一种基于磁共振图像的脑龄测试装置,包括:A second aspect of the present application provides a brain age testing device based on a magnetic resonance image, comprising:
图像获取单元,用于获取待测个体的T1加权大脑磁共振图像;An image acquisition unit, configured to acquire a T1-weighted brain magnetic resonance image of the individual to be tested;
图像分析单元,用于基于所述T1加权大脑磁共振图像确定所述待测个体的脑结构、脑叶以及脑组织分割图谱;An image analyzing unit, configured to determine a brain structure, a brain lobe, and a brain tissue segmentation map of the test subject based on the T1-weighted brain magnetic resonance image;
参数计算单元,用于计算所述脑结构的归一化体积,并基于所述脑组织分割图谱计算所述脑叶的脑萎缩值;a parameter calculation unit, configured to calculate a normalized volume of the brain structure, and calculate a brain atrophy value of the brain lobe based on the brain tissue segmentation map;
脑龄获取单元,用于将所述归一化体积以及所述脑萎缩值输入脑龄估计模型,获取所 述待测个体的脑龄,其中,所述脑龄估计模型根据训练样本个体的归一化体积、脑萎缩值以及脑龄训练得到。a brain age acquisition unit, configured to input the normalized volume and the brain atrophy value into a brain age estimation model, and obtain a brain age of the test subject, wherein the brain age estimation model is based on the individuality of the training sample individual A volume, brain atrophy value and brain age training are obtained.
本申请第三方面提供了一种电子设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现如上所述方法的步骤。A third aspect of the present application provides an electronic device including a memory, a processor, and a computer program stored in the memory and operable on the processor, wherein the processor executes the computer program The steps of the method as described above are achieved.
本申请第四方面提供了一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如上所述方法的步骤。A fourth aspect of the present application provides a computer readable storage medium storing a computer program, wherein the computer program is executed by a processor to implement the steps of the method as described above.
有益效果Beneficial effect
由上可见,本申请提供了一种基于磁共振图像的脑龄测试方法,首先获取待测个体的T1加权大脑磁共振图像,其次,对该T1加权大脑磁共振图像进行分析,获得脑结构、脑叶以及脑组织分割图谱,并根据获得的上述脑结构、脑叶以及脑组织分割图谱计算所述待测个体的归一化体积以及脑萎缩值,将所述归一化体积以及所述脑萎缩值输入脑龄估计模型,获取所述待测个体的脑龄,从而实现对待测个体脑龄的测试。本申请首次将脑结构大小和脑萎缩值作为脑龄测试参数,通过本申请方案,用户可以估算自身的脑龄,以便了解自身大脑健康的状态,方便对自身的大脑健康提前进行干预,以延缓大脑衰老,通过本申请所提供的技术方案可以估计脑退化状况,提高人们对大脑健康的意识。As can be seen from the above, the present application provides a brain age testing method based on magnetic resonance images, first obtaining a T1-weighted brain magnetic resonance image of an individual to be tested, and secondly, analyzing the T1-weighted brain magnetic resonance image to obtain a brain structure, Brain lobe and brain tissue segmentation map, and calculating a normalized volume of the test subject and a brain atrophy value according to the obtained brain structure, brain lobe and brain tissue segmentation map, the normalized volume and the brain The atrophy value is input into the brain age estimation model, and the brain age of the individual to be tested is obtained, thereby realizing the test of the brain age of the individual to be tested. For the first time, the brain size and brain atrophy value are used as brain age test parameters. Through the application of the present application, the user can estimate his or her brain age, so as to understand the state of his brain health, and facilitate the intervention of his brain health in advance to delay. Brain aging, through the technical solutions provided by this application can estimate brain degeneration and improve people's awareness of brain health.
附图说明DRAWINGS
为了更清楚地说明本申请实施例中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其它的附图。In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings used in the embodiments or the prior art description will be briefly described below. Obviously, the drawings in the following description are only the present application. For some embodiments, other drawings may be obtained from those skilled in the art without departing from the drawings.
图1是本申请实施例一提供的基于磁共振图像的脑龄测试方法的实现流程示意图;1 is a schematic flowchart showing an implementation process of a brain age testing method based on magnetic resonance images according to Embodiment 1 of the present application;
图2是本申请实施例二提供的基于磁共振图像的脑龄测试方法的实现流程示意图;2 is a schematic diagram showing an implementation flow of a brain age testing method based on magnetic resonance images provided in Embodiment 2 of the present application;
图3是本申请实施例三提供的基于磁共振图像的脑龄测试装置的结构示意图;3 is a schematic structural diagram of a brain age testing device based on magnetic resonance images provided in Embodiment 3 of the present application;
图4是本申请实施例四提供的基于磁共振图像的脑龄测试装置的结构示意图;4 is a schematic structural diagram of a brain age testing device based on magnetic resonance images provided in Embodiment 4 of the present application;
图5是本申请实施例五提供的电子设备的示意图。FIG. 5 is a schematic diagram of an electronic device according to Embodiment 5 of the present application.
本发明的实施方式Embodiments of the invention
以下描述中,为了说明而不是为了限定,提出了诸如特定系统结构、技术之类的具体细节,以便透彻理解本申请实施例。然而,本领域的技术人员应当清楚,在没有这些具体细节的其它实施例中也可以实现本申请。在其它情况中,省略对众所周知的系统、装置、电路以及方法的详细说明,以免不必要的细节妨碍本申请的描述。In the following description, for purposes of illustration and description However, it will be apparent to those skilled in the art that the present invention may be practiced in other embodiments without these specific details. In other instances, detailed descriptions of well-known systems, devices, circuits, and methods are omitted so as not to obscure the description of the application.
本申请实施例提供的基于磁共振图像的脑龄测试方法适用于电子设备,示例性地,上 述电子设备包括但不限于:台式电脑、平板电脑、云端服务器、手机终端等。The brain age testing method based on the magnetic resonance image provided by the embodiment of the present application is applicable to an electronic device. Illustratively, the above electronic device includes, but is not limited to, a desktop computer, a tablet computer, a cloud server, a mobile phone terminal, and the like.
为了说明本申请上述的技术方案,下面通过具体实施例来进行说明。In order to explain the above technical solutions of the present application, the following description will be made by way of specific embodiments.
实施例1Example 1
下面对本申请实施例一提供的基于磁共振图像的脑龄测试方法进行描述,请参阅附图1,本申请实施例中的脑龄测试方法包括:The method for testing the brain age based on the magnetic resonance image provided in the first embodiment of the present application is described below. Referring to FIG. 1, the brain age testing method in the embodiment of the present application includes:
步骤S101,获取待测个体的T1加权大脑磁共振图像;Step S101, acquiring a T1-weighted brain magnetic resonance image of the individual to be tested;
在本申请实施例中,需要首先获取待测个体的T1加权大脑磁共振图像,以便于后续可以从上述T1加权大脑磁共振图像中获取该待测个体的关于脑发育的相关参数,从而利用该参数估算待测个体的脑龄,其中,所述待测个体的T1加权大脑磁共振图像为三维图像。In the embodiment of the present application, the T1 weighted brain magnetic resonance image of the individual to be tested is first acquired, so that the related parameters about brain development of the test subject can be obtained from the T1 weighted brain magnetic resonance image, thereby utilizing the The parameter estimates the brain age of the individual to be tested, wherein the T1-weighted brain magnetic resonance image of the individual to be tested is a three-dimensional image.
步骤S102,基于上述T1加权大脑磁共振图像确定上述待测个体的脑结构、脑叶以及脑组织分割图谱;Step S102, determining a brain structure, a brain lobe, and a brain tissue segmentation map of the test subject according to the T1 weighted brain magnetic resonance image;
在本申请实施例中,上述脑组织分割图谱包括脑灰质分割图谱、脑白质分割图谱以及脑脊液分割图谱,为了保证可以更精确的确定待测个体的脑结构、脑叶以及脑组织分割图谱,优选地,可以将上述T1加权大脑磁共振图像与预设的脑模板库进行配准。其中,上述脑模板库是预先存储的,该脑模板库中包括两个以上不同大脑的T1加权大脑磁共振图像,即包括两个以上不同的模板,优选地,为保证待测个体脑结构、脑叶以及脑组织分割图谱的准确性,该脑模板库中可包括不同大脑健康程度、不同性别、年龄在10岁~90岁的人的T1加权大脑磁共振图像,且图像个数多于5个。In the embodiment of the present application, the brain tissue segmentation map includes a gray matter segmentation map, a white matter segmentation map, and a cerebrospinal fluid segmentation map. In order to ensure more accurate determination of the brain structure, brain lobe, and brain tissue segmentation map of the individual to be tested, The above T1-weighted brain magnetic resonance image can be registered with a preset brain template library. Wherein, the brain template library is pre-stored, and the brain template library includes T1 weighted brain magnetic resonance images of two or more different brains, that is, including two or more different templates, preferably, to ensure the brain structure of the individual to be tested, The accuracy of the brain leaf and brain tissue segmentation map, which may include T1 weighted brain magnetic resonance images of people with different brain health, gender, age and age between 10 and 90 years old, and the number of images is more than 5 One.
在本申请实施例中,可以事先对脑模板库中的各个T1加权大脑磁共振图像(各个模板)进行脑结构以及脑叶的手动分割,得到脑模板库中各个模板所对应的脑结构以及脑叶,上述脑结构可以包括脑实质、小脑、海马体、杏仁核、脑室、侧脑室、丘脑、尾状核、壳核、苍白球、伏隔核、中脑、桥脑、延脑等和大脑发育和老化相关的重要脑结构,上述脑叶包括左额叶、左顶叶、左枕叶、左颞叶、左扣带回和左岛叶以及右额叶、右顶叶、右枕叶、右颞叶、右扣带回和右岛叶。同时对脑模板库中各个T1加权大脑磁共振图像进行脑组织自动分割,并在电脑自动分割的基础上进行人工手动修正,从而得到脑模板库中各个模板所对应的脑组织概率图谱,上述脑组织概率图谱包括脑白质概率图谱、脑灰质概率图谱以及脑脊液概率图谱。In the embodiment of the present application, the brain structure and the brain leaf can be manually segmented in each T1 weighted brain magnetic resonance image (each template) in the brain template library, and the brain structure and brain corresponding to each template in the brain template library are obtained. Leaf, the above brain structure may include brain parenchyma, cerebellum, hippocampus, amygdala, ventricle, lateral ventricle, thalamus, caudate nucleus, putamen, globus pallidus, nucleus accumbens, midbrain, pons, cerebral ventricle, etc. An important brain structure associated with development and aging, including the left frontal lobe, left parietal lobe, left occipital lobe, left temporal lobe, left cingulate gyrus and left insular lobe, and right frontal lobe, right parietal lobe, right occipital lobe, Right temporal lobe, right cingulate and right island leaves. At the same time, the brain tissue is automatically segmented in each T1 weighted brain magnetic resonance image in the brain template library, and manual manual correction is performed on the basis of automatic computer segmentation, so that the brain tissue probability map corresponding to each template in the brain template library is obtained. The tissue probability map includes a white matter probability map, a gray matter probability map, and a cerebrospinal fluid probability map.
具体的,基于脑模板库获得待测个体的脑结构、脑叶以及脑组织分割图谱可以为:Specifically, the brain structure, the brain lobe, and the brain tissue segmentation map of the individual to be tested based on the brain template library may be:
S1021,利用非线性配准将上述脑模板库中各个模板映射到待测个体的T1加权大脑磁共振图像上,得到该脑模板库中各个模板与待测个体的T1加权大脑磁共振图像之间的空间映射关系;S1021, using a non-linear registration to map each template in the brain template library to a T1-weighted brain magnetic resonance image of the individual to be tested, to obtain a T1-weighted brain magnetic resonance image between each template in the brain template library and the individual to be tested. Spatial mapping relationship;
S1022,利用该空间映射关系,将该脑模板库中各个模板所对应的各个脑结构分区、 脑叶分区以及脑组织概率图谱映射到待测个体的T1加权大脑磁共振图像上,得到各个模板所对应的待测个体的脑结构分区以及脑叶分区,以及待测个体的脑组织概率图谱;S1022, by using the spatial mapping relationship, mapping each brain structure partition, brain leaf partition, and brain tissue probability map corresponding to each template in the brain template library to a T1-weighted brain magnetic resonance image of the individual to be tested, and obtaining each template Corresponding individual brain tissue division and brain leaf division, and brain tissue probability map of the individual to be tested;
其中,上述非线性配准可以采用基于微分同胚模型的对称性非线性配准算法。上述待测个体的脑组织概率图谱为:与待测个体的T1加权大脑磁共振图像相似度最高的模板的脑组织概率图谱经过空间映射得到。Wherein, the above nonlinear registration can adopt a symmetric nonlinear registration algorithm based on a differential homeomorphic model. The brain tissue probability map of the above-mentioned test subject is: the brain tissue probability map of the template with the highest similarity to the T1-weighted brain magnetic resonance image of the individual to be tested is obtained by spatial mapping.
S1023,利用标签融合算法,融合得到最终的上述待测个体的脑结构分区和脑叶分区;S1023, using a label fusion algorithm to fuse the brain structure partition and the brain leaf partition of the final test object;
可以利用标签融合算法,将根据各个模板所得到的待测个体的脑结构分区进行融合得到最终的待测个体的脑结构分区,同理,利用标签融合算法,将根据各个模板所得到的待测个体的脑叶分区进行融合得到最终的待测个体的脑叶分区。其中,标签融合算法是利用期望最大算法将各个模板所对应的分割结果进行组合的算法,通过标签融合算法可以解决单一模板所带来的偏差,使结果更准确。常用的标签融合算法有基于信度的标签融合算法、多数同意规则标签融合算法、带权重的标签融合算法等。The label fusion algorithm can be used to fuse the brain structure partitions of the individuals to be tested according to each template to obtain the brain structure partition of the final test subject. Similarly, using the label fusion algorithm, the test samples obtained according to each template are to be tested. The individual's brain lobe is fused to obtain the brain lobe of the final individual to be tested. The label fusion algorithm is an algorithm that combines the segmentation results corresponding to each template by using a desired maximum algorithm. The label fusion algorithm can solve the deviation caused by a single template, and the result is more accurate. Commonly used tag fusion algorithms include a reliability-based tag fusion algorithm, a majority consent rule tag fusion algorithm, and a weighted tag fusion algorithm.
S1024,将上述待测个体的脑组织概率图谱作为先验知识,利用贝叶斯网络的脑组织分割方法,对待测个体的T1加权大脑磁共振图像进行脑组织分割,得到待测个体的脑组织分割图谱。S1024, using the brain tissue probability map of the above-mentioned test subject as a priori knowledge, using a brain tissue segmentation method of the Bayesian network, and performing brain tissue segmentation on the T1-weighted brain magnetic resonance image of the individual to be measured, and obtaining brain tissue of the individual to be tested Split the map.
步骤S103,计算上述脑结构的归一化体积,并基于上述脑组织分割图谱计算上述脑叶的脑萎缩值;Step S103, calculating a normalized volume of the brain structure, and calculating a brain atrophy value of the brain lobe based on the brain tissue segmentation map;
在本申请实施例中,可以估算步骤S102中获得的待测个体的脑结构的体积,并除以大脑总体积来计算待测个体的脑结构的归一化体积。In the embodiment of the present application, the volume of the brain structure of the individual to be tested obtained in step S102 can be estimated, and the total volume of the brain is divided to calculate the normalized volume of the brain structure of the individual to be tested.
可以根据步骤S102得到的脑叶分区以及脑组织分割图谱,计算该脑叶中的脑灰质体积、脑白质体积以及脑脊液体积,并根据脑萎缩计算公式计算该脑叶的脑萎缩值,其中,上述脑萎缩计算公式为:The brain gray matter volume, the white matter volume, and the cerebrospinal fluid volume in the brain lobe may be calculated according to the brain lobe partition and the brain tissue segmentation map obtained in step S102, and the brain atrophy value of the brain lobe is calculated according to the brain atrophy calculation formula, wherein The formula for brain atrophy is:
Figure PCTCN2018072603-appb-000001
Figure PCTCN2018072603-appb-000001
上述计算公式中得到的脑萎缩值可以直接反映出脑叶的萎缩程度,数值越大,表面该脑叶的脑萎缩程度越大。The brain atrophy value obtained in the above formula can directly reflect the degree of atrophy of the brain lobe. The larger the value, the greater the degree of brain atrophy on the surface of the brain.
步骤S104,将上述归一化体积以及上述脑萎缩值输入脑龄估计模型,获取上述待测个体的脑龄;Step S104, input the normalized volume and the brain atrophy value into a brain age estimation model, and obtain the brain age of the test subject;
在本申请实施例中,可以事先建立脑龄估计模型,首先采集训练样本个体的归一化体积以及脑萎缩值,将上述训练样本个体的归一化体积以及脑萎缩值作为自变量,将训练样本个体的脑龄作为因变量,应用线性支持向量机,建立脑龄估计模型。此外,脑龄估计模型的建立不仅限于线性支持向量机,还可以采用隐藏式马尔科夫模型、类神经网络等来建 立脑龄估计模型。In the embodiment of the present application, the brain age estimation model may be established in advance, firstly, the normalized volume of the training sample individual and the brain atrophy value are collected, and the normalized volume of the training sample individual and the brain atrophy value are used as independent variables, and the training will be performed. The brain age of the sample individuals was used as the dependent variable, and a linear support vector machine was used to establish a brain age estimation model. In addition, the establishment of the brain age estimation model is not limited to the linear support vector machine, but also the hidden Markov model, the neural network, etc. can be used to establish the brain age estimation model.
下面详细论述如何应用线性支持向量机,建立脑龄估计模型。The following is a detailed discussion of how to apply a linear support vector machine to establish a brain age estimation model.
假设训练样本个体为N个,训练数据依次为:(X 1,y 1)、(X 2,y 2)……(X N,y N),其中X i=(x i,1,x i,2…x i,K) Τ,i=1,2…N,X i包括K个与脑发育相关的参数,在本申请实施例中,可以为多个脑结构的归一化体积以及多个脑叶的脑萎缩值,y i为第i个训练样本个体所对应的脑龄。 Assuming that there are N training samples, the training data are: (X 1 , y 1 ), (X 2 , y 2 ), ... (X N , y N ), where X i = (x i,1 ,x i , 2 ... x i, K ) Τ , i = 1, 2 ... N, X i includes K parameters related to brain development, and in the embodiment of the present application, may be a normalized volume of multiple brain structures and The brain atrophy value of the brain lobe, y i is the brain age corresponding to the i-th training sample individual.
假设脑龄估计模型的表达式为:f(X)=W·X+b,W∈R 1×K,b∈R Assume that the expression of the brain age estimation model is: f(X)=W·X+b, W∈R 1×K , b∈R
优化目标为:
Figure PCTCN2018072603-appb-000002
求取优化目标值最小时W与b的具体数值,可利用梯度下降算法或遗传算法等数据拟合算法,求解出W与b的数值,从而得到脑龄估计模型,此外,优化目标还可以为
Figure PCTCN2018072603-appb-000003
此处对优化目标的公式不作限定。
The optimization goal is:
Figure PCTCN2018072603-appb-000002
The specific values of W and b when the optimization target value is minimum can be obtained. The data fitting algorithm such as gradient descent algorithm or genetic algorithm can be used to solve the values of W and b, so as to obtain the brain age estimation model. In addition, the optimization target can also be
Figure PCTCN2018072603-appb-000003
The formula for optimizing the target is not limited here.
优选地,为了更精确地估计待测个体的脑龄,还可以在获取待测个体的T1加权大脑磁共振图像之后,对该T1加权大脑磁共振图像进行预处理,上述预处理包括如下一项或两项以上:降噪、去偏场以及像素范围归一化。本申请基于灰度分布偏度的噪声估计方法来估计莱斯噪声方差,然后根据估计的噪声方差,利用非局部均值算法降噪;本申请预先去偏场用以去除不均匀磁场所造成的同种脑组织灰度不一致的现象;本申请利用柱状图匹配的方法,把不同图像的强度范围归一化到一个共同的范围。Preferably, in order to more accurately estimate the brain age of the individual to be tested, the T1-weighted brain magnetic resonance image may also be pre-processed after acquiring the T1-weighted brain magnetic resonance image of the individual to be tested, the pre-processing including the following Or more than two: noise reduction, de-fielding, and pixel range normalization. The present invention estimates the Rice noise variance based on the gray-scale distribution skewness noise estimation method, and then uses the non-local mean algorithm to reduce noise according to the estimated noise variance; the present application de-biasing field is used to remove the uneven magnetic field. The phenomenon that the gray level of the brain tissue is inconsistent; this application uses the method of histogram matching to normalize the intensity ranges of different images to a common range.
收集1936个健康训练样本个体,男女数量比例为1:1,年龄跨度从40岁到90岁,其中训练样本个体的T1加权大脑磁共振图像来自中国16家不同医院的不同机型,包括SIEMENS,GE,PHILIPS核磁共振机器。首先对收集的训练样本个体的T1加权大脑磁共振图像进行预处理;获取每个训练样本个体所对应的38个归一化体积以及12个脑萎缩值;由于健康个体的脑龄与生理年龄近乎相等,因此在训练脑龄估计模型时,将训练样本个体的生理年龄作为训练样本个体的脑龄来建立脑龄估计模型。经训练,所建立的脑龄估计模型与其生理年龄的差异为5.44岁,也即按照本申请实施例一所提供的技术方案,所建立的脑龄估计模型的误差为5.44岁。待测个体为12个健康个体(年龄范围为67.3±9.3岁)和14个患有老年痴呆症的个体(年龄范围为61.3±15.6岁),待测个体的T1加权大脑磁共振图像来自PHILIPS核磁共振机器,首先对待测个体的T1加权大脑磁共振图像进行预处理,并提取待测个体的38个归一化脑体积以及12个脑萎缩值,根据我们上述获得的脑龄估计模型,来估计待测个体的脑龄。测试结果为:上述12个健康个体的预测脑龄与其生理年龄平均差异为6.4岁,患有老年痴呆症的个体的预测脑龄和其生理年龄差异为19.7岁。从上 述测试结果中可以看出老年痴呆症患者的大脑比健康个体的大脑老化程度高。A total of 1936 healthy training samples were collected. The ratio of male to female was 1:1, and the age span ranged from 40 to 90 years. The T1 weighted brain magnetic resonance images of the training samples were from different models of 16 different hospitals in China, including SIEMENS. GE, PHILIPS NMR machine. Firstly, the T1 weighted brain magnetic resonance images of the collected training sample individuals are preprocessed; the 38 normalized volumes corresponding to each training sample individual and 12 brain atrophy values are obtained; since the healthy individual's brain age and physiological age are close Equal, so when training the brain age estimation model, the physiological age of the training sample individual is used as the brain age of the training sample individual to establish a brain age estimation model. After training, the difference between the established brain age estimation model and its physiological age is 5.44 years old, that is, according to the technical solution provided by the first embodiment of the present application, the error of the established brain age estimation model is 5.44 years old. The individuals to be tested were 12 healthy individuals (age range 67.3±9.3 years) and 14 individuals with Alzheimer's disease (age range 61.3±15.6 years). The T1-weighted brain magnetic resonance images of the individuals to be tested were from PHILIPS NMR. The resonance machine first preprocesses the T1-weighted brain magnetic resonance image of the individual to be measured, and extracts 38 normalized brain volumes and 12 brain atrophy values of the individuals to be tested, and estimates according to the brain age estimation model obtained above. The age of the individual to be tested. The test results showed that the average age difference between the predicted brain age and the physiological age of the above 12 healthy individuals was 6.4 years, and the predicted brain age and the physiological age difference of the individuals with Alzheimer's disease were 19.7 years old. It can be seen from the above test results that the brain of the Alzheimer's disease patient has a higher degree of brain aging than the healthy individual.
在本申请实施例中,首次将脑结构大小和脑萎缩值作为脑龄测试参数,通过本申请方案,用户可以估算自身的脑龄,以便了解自身大脑健康的状态,方便对自身的大脑健康提前进行干预,以延缓大脑衰老,通过本申请所提供的技术方案可以估计脑退化状况,提高人们对大脑健康的意识。In the embodiment of the present application, the brain structure size and the brain atrophy value are used as the brain age test parameters for the first time. Through the application scheme, the user can estimate the brain age of the user to understand the state of the brain health and facilitate the advancement of the brain health. Intervention to delay brain aging, through the technical solutions provided in this application can estimate brain degeneration and improve people's awareness of brain health.
实施例2Example 2
下面对本申请实施例二提供的另一种基于磁共振图像的脑龄测试方法进行描述,请参阅附图2。本申请实施例中的脑龄测试方法包括:Another method for testing a brain age based on magnetic resonance images provided in the second embodiment of the present application is described below. Please refer to FIG. 2 . The brain age testing method in the embodiment of the present application includes:
步骤S201,获取训练样本个体的T1加权大脑磁共振图像;Step S201, acquiring a T1-weighted brain magnetic resonance image of the training sample individual;
在本申请实施例中,具体给出了一种脑龄估计模型的建立方法。为了确定一种比较精确的估计获取模型,常常需要大量的与脑发育相关的参数,然而,很多与脑发育相关的参数并非是必要的,即便在脑龄估计模型中增加该参数,也并不能使脑龄估计模型的精确度增加太多,反而会在后续计算待测个体脑龄的过程中占用大量的计算资源。因此,该实施例给出一种脑龄估计模型的建立方法,可以剔除建立脑龄估计模型时不需要的参数,释放一定的计算资源。In the embodiment of the present application, a method for establishing a brain age estimation model is specifically given. In order to determine a more accurate estimation acquisition model, a large number of parameters related to brain development are often required. However, many parameters related to brain development are not necessary, even if this parameter is added to the brain age estimation model, it cannot be To increase the accuracy of the brain age estimation model too much, it will occupy a large amount of computing resources in the process of calculating the brain age of the individual to be tested. Therefore, this embodiment provides a method for establishing a brain age estimation model, which can eliminate parameters that are not needed when establishing a brain age estimation model, and release certain computing resources.
为建立脑龄估计模型,首先需要获取训练样本个体的T1加权大脑磁共振图像,以便于后续提取训练样本个体的关于脑发育的相关参数,从而利用该参数建立脑龄估计模型。In order to establish a brain age estimation model, it is first necessary to obtain a T1-weighted brain magnetic resonance image of the training sample individual, so as to subsequently extract relevant parameters of brain development of the training sample individual, thereby using the parameter to establish a brain age estimation model.
为了更精确地提取训练样本个体的关于脑发育的相关参数,可以对各个训练样本个体的T1加权大脑磁共振图像进行预处理,上述预处理包括如下一项或两项以上:降噪、去偏场以及像素范围归一化。In order to more accurately extract the relevant parameters of brain development of the training sample individual, the T1-weighted brain magnetic resonance image of each training sample individual may be pre-processed, and the above pre-processing includes one or more of the following: noise reduction, de-biasing The field and pixel range are normalized.
步骤S202,获取各个训练样本个体的脑叶以及脑组织分割图谱,基于上述脑组织分割图谱,计算该脑叶的脑萎缩值;Step S202, acquiring a brain lobe and a brain tissue segmentation map of each training sample individual, and calculating a brain atrophy value of the brain lobe based on the brain tissue segmentation map;
在本申请实施例中,根据各个训练样本个体的T1加权大脑磁共振图像,获取各个训练样本个体的脑叶以及脑组织概率图谱,并计算该脑叶的脑萎缩值,具体的实施方法可参见实施例一S102以及S103上述,此处不再赘述。In the embodiment of the present application, according to the T1-weighted brain magnetic resonance image of each training sample individual, the brain lobe and the brain tissue probability map of each training sample individual are obtained, and the brain atrophy value of the brain lobe is calculated, and the specific implementation method can be referred to The first embodiment S102 and S103 are as described above, and are not described herein again.
在该实施例中,优选获取各个训练样本个体所对应的12个脑萎缩值,包括左额叶、左顶叶、左枕叶、左颞叶、左扣带回和左岛叶以及右额叶、右顶叶、右枕叶、右颞叶、右扣带回和右岛叶的脑萎缩值。In this embodiment, it is preferred to obtain 12 brain atrophy values corresponding to each training sample individual, including left frontal lobe, left parietal lobe, left occipital lobe, left temporal lobe, left cingulate gyrus and left insular lobe, and right frontal lobe. Brain atrophy values of right parietal lobe, right occipital lobe, right temporal lobe, right cingulate gyrus, and right island lobes.
步骤S203,获取各个训练样本个体的第i个脑结构,计算上述第i个脑结构的归一化体积,基于第i个脑结构的归一化体积,以及训练样本个体的脑叶的脑萎缩值,建立第i个脑龄估计模型,并获取上述第i个脑龄估计模型的误差;Step S203, acquiring an i-th brain structure of each training sample individual, calculating a normalized volume of the i-th brain structure, a normalized volume based on the i-th brain structure, and brain atrophy of the brain of the training sample individual Value, establish an i-th brain age estimation model, and obtain the error of the above i-th brain age estimation model;
在本申请实施例中,大脑中的脑结构非常多,但是有些脑结构的大小随着大脑发育变 化很小,也有些脑结构的大小变化并不影响大脑的老化程度,显然,这些脑结构的归一化体积在脑龄计算中并不需要。因此,在本申请实施例中,我们可以首先选取某一个脑结构,计算所选的脑结构的归一化体积,将该归一化体积以及上述12个脑叶萎缩值作为建立脑龄估计模型的13个因变量,可以应用线性支持向量机建立该脑结构所对应的脑龄估计模型,并将优化目标的最小值作为该脑龄估计模型的误差。In the embodiment of the present application, the brain structure in the brain is very large, but the size of some brain structures changes little with the development of the brain, and some changes in the size of the brain structure do not affect the degree of aging of the brain. Obviously, the structure of these brains The normalized volume is not required in brain age calculations. Therefore, in the embodiment of the present application, we can first select a certain brain structure, calculate the normalized volume of the selected brain structure, and use the normalized volume and the above 12 brain leaf atrophy values as the model for establishing the brain age. The 13 dependent variables can be used to establish the brain age estimation model corresponding to the brain structure by using the linear support vector machine, and the minimum value of the optimization target is used as the error of the brain age estimation model.
步骤S204,判断是否遍历完每一个训练样本个体的所有脑结构;Step S204, determining whether to traverse all brain structures of each training sample individual;
在本申请实施例中,步骤S203之后,判断是否遍历完所有的脑结构,若否,则执行步骤S205,若是,执行步骤S206。In the embodiment of the present application, after step S203, it is determined whether all brain structures have been traversed. If not, step S205 is performed, and if yes, step S206 is performed.
步骤S205,使i值增加1;Step S205, increasing the value of i by one;
在该步骤中,若没有遍历完所有的脑结构,则将下一个脑结构的归一化体积以及上述12个脑叶萎缩值作为建立脑龄估计模型的13个因变量,建立下一个脑结构所对应的脑龄估计模型,并获取该脑龄估计模型的误差。In this step, if all the brain structures have not been traversed, the normalized volume of the next brain structure and the above 12 brain leaf atrophy values are used as the 13 dependent variables for establishing the brain age estimation model, and the next brain structure is established. The corresponding brain age estimation model is obtained, and the error of the brain age estimation model is obtained.
步骤S206,选取误差较小的脑龄估计模型所对应的脑结构的归一化体积,并根据选取后的脑结构的归一化体积以及脑叶萎缩值,建立最终的脑龄估计模型;Step S206, selecting a normalized volume of the brain structure corresponding to the brain age estimation model with a small error, and establishing a final brain age estimation model according to the normalized volume of the selected brain structure and the brain blade atrophy value;
在该步骤中,若遍历完所有的脑结构,则选取误差较小的脑龄估计模型所对应的脑结构的归一化体积。可以预先设置阈值,比如为10岁,将上述获得的各个脑结构所对应的脑龄估计模型的误差与10岁相比较,若小于预设的10岁,则选取该脑结构的归一化体积,比如,可能选取出的脑结构归一化体积仅仅只有小脑归一化体积以及海马体归一化体积,将选取后的脑结构的归一化体积即小脑归一化体积以及海马体归一化体积,以及12个脑叶萎缩值作为最终的脑龄估计模型的14个因变量,建立最终的脑龄估计模型。In this step, if all brain structures are traversed, the normalized volume of the brain structure corresponding to the brain age estimation model with less error is selected. The threshold may be preset, for example, 10 years old, and the error of the brain age estimation model corresponding to each brain structure obtained above is compared with 10 years old, and if it is less than the preset 10 years old, the normalized volume of the brain structure is selected. For example, the normalized volume of the brain structure that may be selected is only the normalized volume of the cerebellum and the normalized volume of the hippocampus. The normalized volume of the selected brain structure, ie, the normalized volume of the cerebellum and the hippocampus normalized. The volume, and 12 brain leaf atrophy values were used as the 14 dependent variables of the final brain age estimation model to establish the final brain age estimation model.
S207,获取待测个体的T1加权大脑磁共振图像;S207. Acquire a T1-weighted brain magnetic resonance image of the individual to be tested;
S208,基于上述待测个体的T1加权大脑磁共振图像确定上述待测个体的脑结构、脑叶以及脑组织分割图谱;S208, determining a brain structure, a brain lobe, and a brain tissue segmentation map of the test subject according to the T1-weighted brain magnetic resonance image of the test subject;
S209,计算上述待测个体的脑结构的归一化体积,并基于上述待测个体的脑组织分割图谱计算上述待测个体脑叶的脑萎缩值;S209, calculating a normalized volume of the brain structure of the test subject, and calculating a brain atrophy value of the brain of the test subject according to the brain tissue segmentation map of the test subject;
S210,将上述待测个体的归一化体积以及上述待测个体的脑萎缩值输入最终的脑龄估计模型,获取上述待测个体的脑龄;S210, inputting the normalized volume of the individual to be tested and the brain atrophy value of the test subject to the final brain age estimation model, and obtaining the brain age of the test subject;
在本申请实施例中,确定的上述待测个体的脑结构是与上述最终的脑龄估计模型中的脑结构对应的,比如,若最终的脑龄估计模型中只有小脑归一化体积与海马体归一化体积,则步骤S208中,只需要确定上述待测个体的小脑以及海马体,步骤S209中,只需要计算小脑归一化体积与脑结构归一化体积。此外,上述步骤S207-S209与实施例一中的步骤S101-S103实施方式相同,具体可参见实施例一的描述,此处不再赘述。此外,在本申请实 施例中,并非每一次脑龄测试都需要执行步骤S201-S206,在最终脑龄估计模型建立之后,后续进行脑龄测试时无需再执行步骤S201-S206。In the embodiment of the present application, the determined brain structure of the test subject is corresponding to the brain structure in the final brain age estimation model, for example, if the final brain age estimation model has only the cerebellar normalized volume and the hippocampus If the volume is normalized, then in step S208, only the cerebellum and hippocampus of the individual to be tested need to be determined. In step S209, only the normalized volume of the cerebellum and the normalized volume of the brain structure need to be calculated. In addition, the foregoing steps S207-S209 are the same as the implementations of the steps S101-S103 in the first embodiment. For details, refer to the description of the first embodiment, and details are not described herein again. In addition, in the embodiment of the present application, not every brain age test needs to perform steps S201-S206. After the final brain age estimation model is established, steps S201-S206 need not be performed in subsequent brain age test.
在本申请实施例中,具体给出了一种脑龄估计模型的建立方法,可以剔除建立脑龄估计模型时不需要的参数,释放一定的计算资源;并且首次将脑结构大小和脑萎缩值作为脑龄测试参数,通过本申请方案,用户可以估算自身的脑龄,以便了解自身大脑健康的状态,方便对自身的大脑健康提前进行干预,以延缓大脑衰老,通过本申请所提供的技术方案可以估计脑退化状况,提高人们对大脑健康的意识。In the embodiment of the present application, a method for establishing a brain age estimation model is specifically provided, which can eliminate parameters that are not needed when establishing a brain age estimation model, release certain computing resources, and firstly determine brain structure size and brain atrophy value. As a brain age test parameter, the user can estimate the age of the brain by using the solution of the present application, in order to understand the state of the brain health, and facilitate the intervention of the brain health in advance to delay the aging of the brain, through the technical solution provided by the application. Brain degeneration can be estimated and people's awareness of brain health can be improved.
应理解,上述实施例中各步骤的序号的大小并不意味着执行顺序的先后,各过程的执行顺序应以其功能和内在逻辑确定,而不应对本申请实施例的实施过程构成任何限定。It should be understood that the size of the sequence of the steps in the above embodiments does not mean that the order of execution is performed. The order of execution of each process should be determined by its function and internal logic, and should not be construed as limiting the implementation process of the embodiments of the present application.
实施例3Example 3
本申请实施例三提供一种基于磁共振图像的脑龄测试装置,为了便于说明,仅示出与本申请相关的部分,如图3所示,上述脑龄测试装置300包括:A third embodiment of the present invention provides a brain age testing device based on a magnetic resonance image. For convenience of description, only parts related to the present application are shown. As shown in FIG. 3, the brain age testing device 300 includes:
图像获取单元301,用于获取待测个体的T1加权大脑磁共振图像;The image obtaining unit 301 is configured to acquire a T1-weighted brain magnetic resonance image of the individual to be tested;
图像分析单元302,用于基于上述T1加权大脑磁共振图像确定上述待测个体的脑结构、脑叶以及脑组织分割图谱,所述脑组织分割图谱包括脑灰质分割图谱、脑白质分割图谱以及脑脊液分割图谱;The image analyzing unit 302 is configured to determine a brain structure, a brain lobe, and a brain tissue segmentation map of the test subject according to the T1 weighted brain magnetic resonance image, where the brain tissue segmentation map includes a gray matter segmentation map, a white matter segmentation map, and a cerebrospinal fluid Split the map;
参数计算单元303,用于计算上述脑结构的归一化体积,并基于上述脑组织分割图谱计算上述脑叶的脑萎缩值;The parameter calculation unit 303 is configured to calculate a normalized volume of the brain structure, and calculate a brain atrophy value of the brain lobe based on the brain tissue segmentation map;
脑龄获取单元304,用于将上述归一化体积以及上述脑萎缩值输入脑龄估计模型,获取上述待测个体的脑龄,其中,上述脑龄估计模型根据训练样本个体的归一化体积、脑萎缩值以及脑龄训练得到。The brain age obtaining unit 304 is configured to input the normalized volume and the brain atrophy value into a brain age estimation model to obtain a brain age of the test subject, wherein the brain age estimation model is based on a normalized volume of the training sample individual Brain atrophy and brain age training are obtained.
优选地,上述图像分析单元302具体用于:Preferably, the image analysis unit 302 is specifically configured to:
将上述T1加权大脑磁共振图像与预设的脑模板库进行配准,获得上述T1加权大脑磁共振图像所对应的脑结构、脑叶以及脑组织概率图谱;将所述T1加权大脑磁共振图像所对应的脑组织概率图谱作为先验知识,对所述T1加权大脑磁共振图像进行脑组织分割,得到所述T1加权大脑磁共振图像所对应的脑组织分割图谱;The T1 weighted brain magnetic resonance image is registered with a preset brain template library to obtain a brain structure, a brain lobe, and a brain tissue probability map corresponding to the T1 weighted brain magnetic resonance image; and the T1 weighted brain magnetic resonance image is obtained. Corresponding brain tissue probability map is used as a priori knowledge, and brain tissue segmentation is performed on the T1-weighted brain magnetic resonance image to obtain a brain tissue segmentation map corresponding to the T1-weighted brain magnetic resonance image;
其中,上述脑模板库中包括:两个以上不同大脑的T1加权大脑磁共振图像,以及分别与各个T1加权大脑磁共振图像对应的脑结构、脑叶以及脑组织概率图谱,上述脑组织概率图谱包括脑灰质概率图谱、脑白质概率图谱和脑脊液概率图谱。The brain template library includes: a T1-weighted brain magnetic resonance image of two or more different brains, and a brain structure, a brain lobe, and a brain tissue probability map corresponding to each T1-weighted brain magnetic resonance image, and the brain tissue probability map. Including the gray matter probability map, the white matter probability map and the cerebrospinal fluid probability map.
优选地,上述参数计算单元303具体用于:Preferably, the parameter calculation unit 303 is specifically configured to:
计算上述脑结构的归一化体积,并基于上述待测个体的T1加权大脑磁共振图像的脑组织分割图谱,确定上述待测个体的T1加权大脑磁共振图像的上述脑叶中的脑灰质体积、 脑白质体积以及脑脊液体积;Calculating a normalized volume of the brain structure, and determining a gray matter volume in the brain lobe of the T1-weighted brain magnetic resonance image of the individual to be tested based on a brain tissue segmentation map of the T1-weighted brain magnetic resonance image of the individual to be tested , white matter volume and cerebrospinal fluid volume;
根据脑萎缩计算公式计算上述脑叶的脑萎缩值,其中,上述脑萎缩计算公式为:The brain atrophy value of the brain lobe is calculated according to the brain atrophy calculation formula, wherein the above brain atrophy calculation formula is:
Figure PCTCN2018072603-appb-000004
Figure PCTCN2018072603-appb-000004
优选地,上述脑龄测试装置300还包括:Preferably, the brain age testing device 300 further includes:
模型建立单元,用于将训练样本个体的归一化体积以及脑萎缩值作为自变量,以训练样本个体的脑龄作为因变量,应用线性支持向量机,建立脑龄估计模型。The model building unit is configured to use the normalized volume of the training sample individual and the brain atrophy value as independent variables, and use the brain support age of the training individual as a dependent variable, and apply a linear support vector machine to establish a brain age estimation model.
优选地,上述脑龄测试装置300还包括:Preferably, the brain age testing device 300 further includes:
待测个体预处理单元,用于在上述T1加权大脑磁共振图像确定上述待测个体的脑结构、脑叶以及脑组织分割图谱之前,对上述待测个体的T1加权大脑磁共振图像进行预处理,上述预处理包括如下一项或两项以上:降噪、去偏场、像素范围归一化。The individual preprocessing unit to be tested is configured to preprocess the T1-weighted brain magnetic resonance image of the test subject to be tested before determining the brain structure, the brain lobe and the brain tissue segmentation map of the test subject according to the T1 weighted brain magnetic resonance image. The above preprocessing includes one or more of the following: noise reduction, de-biasing, and pixel range normalization.
需要说明的是,上述装置各模块/单元之间的信息交互、执行过程等内容,由于与本申请方法实施例一基于同一构思,其具体功能及带来的技术效果,具体可参见方法实施例一部分,此处不再赘述。It should be noted that the information interaction, the execution process, and the like between the modules/units of the foregoing device are based on the same concept, the specific functions and the technical effects brought about by the method embodiment of the present application. Part of it will not be described here.
实施例4Example 4
本申请实施例四提供一种基于磁共振图像的脑龄测试装置,为了便于说明,仅示出与本申请相关的部分,如图4所示,上述脑龄测试装置400包括:A fourth embodiment of the present invention provides a brain age testing device based on a magnetic resonance image. For convenience of explanation, only parts related to the present application are shown. As shown in FIG. 4, the brain age testing device 400 includes:
训练样本个体图像获取单元401,用于获取训练样本个体的T1加权大脑磁共振图像;a training sample individual image acquiring unit 401, configured to acquire a T1-weighted brain magnetic resonance image of the training sample individual;
训练样本个体脑萎缩计算单元402,用于获取各个训练样本个体的脑叶以及脑组织分割图谱,基于上述脑组织分割图谱,计算脑叶的脑萎缩值;The training sample individual brain atrophy calculation unit 402 is configured to acquire a brain lobe and a brain tissue segmentation map of each training sample individual, and calculate a brain atrophy value of the brain lobe based on the brain tissue segmentation map;
第i个脑龄模型获取单元403,用于获取各个训练样本个体的第i个脑结构,计算上述第i个脑结构的归一化体积,基于第i个脑结构的归一化体积,以及脑叶的脑萎缩值,建立第i个脑龄估计模型,并获取上述第i个脑龄估计模型的误差;The i-th brain age model acquisition unit 403 is configured to acquire an i-th brain structure of each training sample individual, calculate a normalized volume of the i-th brain structure, based on a normalized volume of the i-th brain structure, and Brain atrophy value of the brain, establish an i-th brain age estimation model, and obtain the error of the above i-th brain age estimation model;
判断单元404,用于判断是否遍历完每一个训练样本个体的所有脑结构;The determining unit 404 is configured to determine whether to traverse all brain structures of each training sample individual;
计算序号增加单元405,用于在没有遍历完所有脑结构的情况下,使i值增加1;Calculating a sequence number increasing unit 405 for increasing the value of i by 1 without traversing all brain structures;
最终模型建立单元406,用于选取误差较小的脑龄估计模型所对应的脑结构的归一化体积,并根据选取后的脑结构的归一化体积以及脑叶萎缩值,建立最终的脑龄估计模型;The final model establishing unit 406 is configured to select a normalized volume of the brain structure corresponding to the brain age estimation model with less error, and establish a final brain according to the normalized volume of the selected brain structure and the brain blade atrophy value. Age estimation model;
待测个体图像获取单元407,用于获取待测个体的T1加权大脑磁共振图像;The individual image acquisition unit 407 to be tested is configured to acquire a T1-weighted brain magnetic resonance image of the individual to be tested;
待测个体图像分析单元408,用于基于上述待测个体的T1加权大脑磁共振图像确定上述待测个体的脑结构、脑叶以及脑组织分割图谱;The individual image analysis unit 408 to be tested is configured to determine a brain structure, a brain lobe, and a brain tissue segmentation map of the individual to be tested based on the T1 weighted brain magnetic resonance image of the individual to be tested;
待测个体参数计算单元409,用于计算上述待测个体的脑结构的归一化体积,并基于上述待测个体的脑组织分割图谱计算上述待测个体脑叶的脑萎缩值;The individual parameter calculation unit 409 to be tested is configured to calculate a normalized volume of the brain structure of the individual to be tested, and calculate a brain atrophy value of the brain of the individual to be tested based on the brain tissue segmentation map of the individual to be tested;
待测个体脑龄获取单元410,用于将上述待测个体的归一化体积以及上述待测个体的脑萎缩值输入最终的脑龄估计模型,获取上述待测个体的脑龄;The individual brain age obtaining unit 410 is configured to input the normalized volume of the test subject and the brain atrophy value of the test subject into a final brain age estimation model, and obtain the brain age of the test subject;
在本申请实施例中,确定的上述待测个体的脑结构是与上述最终的脑龄估计模型中的脑结构对应的,比如,若最终的脑龄估计模型中只有小脑归一化体积与海马体归一化体积,则待测个体图像分析单元408,只需要确定上述待测个体的小脑以及海马体,待测个体参数计算单元409,只需要计算小脑归一化体积与脑结构归一化体积。此外,上述单元407-409与实施例三中的单元301-303实施方式相同,具体可参见实施例三的描述,此处不再赘述。In the embodiment of the present application, the determined brain structure of the test subject is corresponding to the brain structure in the final brain age estimation model, for example, if the final brain age estimation model has only the cerebellar normalized volume and the hippocampus For the normalized volume, the individual image analysis unit 408 to be tested only needs to determine the cerebellum and hippocampus of the individual to be tested, and the individual parameter calculation unit 409 to be tested only needs to calculate the normalized volume of the cerebellum and the normalization of the brain structure. volume. In addition, the foregoing units 407-409 are the same as the implementations of the units 301-303 in the third embodiment. For details, refer to the description of the third embodiment, and details are not described herein again.
优选地,上述脑龄测试装置400还包括:Preferably, the brain age testing device 400 further includes:
预处理单元,用于在确定上述待测个体的脑结构、脑叶以及脑组织分割图谱之前,对上述待测个体的T1加权大脑磁共振图像进行预处理,以及用于在获取各个训练样本个体的脑结构、脑叶以及脑组织分割图谱之前,对各个训练样本个体的T1加权大脑磁共振图像进行预处理,上述预处理包括如下一项或两项以上:降噪、去偏场、像素范围归一化。a pre-processing unit, configured to pre-process the T1-weighted brain magnetic resonance image of the individual to be tested, and to obtain each training sample individual, before determining the brain structure, the brain lobe, and the brain tissue segmentation map of the individual to be tested Before the brain structure, brain lobe and brain tissue segmentation map, the T1-weighted brain magnetic resonance images of each training sample individual are preprocessed, and the above preprocessing includes one or more of the following: noise reduction, de-biasing field, pixel range Normalized.
需要说明的是,上述装置各模块/单元之间的信息交互、执行过程等内容,由于与本申请方法实施例二基于同一构思,其具体功能及带来的技术效果,具体可参见方法实施例二部分,此处不再赘述。It should be noted that the content of the information exchange, the execution process, and the like between the modules/units of the foregoing device are based on the same concept, the specific functions and the technical effects brought by the second embodiment of the present application. The second part is not repeated here.
实施例5Example 5
图5是本申请实施例五提供的电子设备的示意图。如图5所示,该实施例的电子设备5包括:处理器50、存储器51以及存储在上述存储器51中并可在上述处理器50上运行的计算机程序52。上述处理器50执行上述计算机程序52时实现上述各个方法实施例中的步骤,例如图1所示的步骤S101至S104。或者,上述处理器50执行上述计算机程序52时实现上述各装置实施例中各模块/单元的功能,例如图3所示模块301至304的功能。FIG. 5 is a schematic diagram of an electronic device according to Embodiment 5 of the present application. As shown in FIG. 5, the electronic device 5 of this embodiment includes a processor 50, a memory 51, and a computer program 52 stored in the above-described memory 51 and operable on the processor 50 described above. When the processor 50 executes the computer program 52 described above, the steps in the above various method embodiments are implemented, such as steps S101 to S104 shown in FIG. Alternatively, when the processor 50 executes the computer program 52, the functions of the modules/units in the above-described respective device embodiments, such as the functions of the modules 301 to 304 shown in FIG. 3, are implemented.
上述电子设备5可以是桌上型计算机、笔记本、掌上电脑及云端服务器等计算设备。上述电子设备可包括,但不仅限于,处理器50、存储器51。本领域技术人员可以理解,图5仅仅是电子设备5的示例,并不构成对电子设备5的限定,可以包括比图示更多或更少的部件,或者组合某些部件,或者不同的部件,例如上述电子设备5还可以包括输入输出设备、网络接入设备、总线等。The above electronic device 5 may be a computing device such as a desktop computer, a notebook, a palmtop computer, and a cloud server. The above electronic device may include, but is not limited to, the processor 50 and the memory 51. It will be understood by those skilled in the art that FIG. 5 is only an example of the electronic device 5, and does not constitute a limitation on the electronic device 5, and may include more or less components than those illustrated, or combine some components, or different components. For example, the electronic device 5 described above may further include an input/output device, a network access device, a bus, and the like.
所称处理器50可以是中央处理单元(Central Processing Unit,CPU),还可以是其它通用处理器、数字信号处理器(Digital Signal Processor,DSP)、专用集成电路(Application Specific Integrated Circuit,ASIC)、现场可编程门阵列(Field-Programmable Gate Array,FPGA)或者其它可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。The processor 50 may be a central processing unit (CPU), or may be other general-purpose processors, a digital signal processor (DSP), an application specific integrated circuit (ASIC), Field-Programmable Gate Array (FPGA) or other programmable logic device, discrete gate or transistor logic device, discrete hardware components, etc. The general purpose processor may be a microprocessor or the processor or any conventional processor or the like.
上述存储器51可以是上述电子设备5的内部存储单元,例如电子设备5的硬盘或内 存。上述存储器51也可以是上述电子设备5的外部存储设备,例如上述电子设备5上配备的插接式硬盘,智能存储卡(Smart Media Card,SMC),安全数字(Secure Digital,SD)卡,闪存卡(Flash Card)等。进一步地,上述存储器51还可以既包括上述电子设备5的内部存储单元也包括外部存储设备。上述存储器51用于存储上述计算机程序以及上述电子设备所需的其它程序和数据。上述存储器51还可以用于暂时地存储已经输出或者将要输出的数据。The memory 51 may be an internal storage unit of the electronic device 5, such as a hard disk or an internal memory of the electronic device 5. The memory 51 may be an external storage device of the electronic device 5, such as a plug-in hard disk equipped with the above-mentioned electronic device 5, a smart memory card (SMC), a Secure Digital (SD) card, and a flash memory. Flash card, etc. Further, the above-mentioned memory 51 may also include both an internal storage unit of the above-described electronic device 5 and an external storage device. The above memory 51 is used to store the above computer program and other programs and data required for the above electronic device. The above-described memory 51 can also be used to temporarily store data that has been output or is about to be output.
所属领域的技术人员可以清楚地了解到,为了描述的方便和简洁,仅以上述各功能单元、模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能单元、模块完成,即将上述装置的内部结构划分成不同的功能单元或模块,以完成以上描述的全部或者部分功能。实施例中的各功能单元、模块可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中,上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。另外,各功能单元、模块的具体名称也只是为了便于相互区分,并不用于限制本申请的保护范围。上述系统中单元、模块的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。It will be apparent to those skilled in the art that, for convenience and brevity of description, only the division of each functional unit and module described above is exemplified. In practical applications, the above functions may be assigned to different functional units as needed. The module is completed by dividing the internal structure of the above device into different functional units or modules to perform all or part of the functions described above. Each functional unit and module in the embodiment may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit, and the integrated unit may be hardware. Formal implementation can also be implemented in the form of software functional units. In addition, the specific names of the respective functional units and modules are only for the purpose of facilitating mutual differentiation, and are not intended to limit the scope of protection of the present application. For the specific working process of the unit and the module in the foregoing system, reference may be made to the corresponding process in the foregoing method embodiment, and details are not described herein again.
在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述或记载的部分,可以参见其它实施例的相关描述。In the above embodiments, the descriptions of the various embodiments are different, and the parts that are not detailed or described in a certain embodiment can be referred to the related descriptions of other embodiments.
本领域普通技术人员可以意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本申请的范围。Those of ordinary skill in the art will appreciate that the elements and algorithm steps of the various examples described in connection with the embodiments disclosed herein can be implemented in electronic hardware or a combination of computer software and electronic hardware. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the solution. A person skilled in the art can use different methods to implement the described functions for each particular application, but such implementation should not be considered to be beyond the scope of the present application.
在本申请所提供的实施例中,应该理解到,所揭露的装置/电子设备和方法,可以通过其它的方式实现。例如,以上所描述的装置/电子设备实施例仅仅是示意性的,例如,上述模块或单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通讯连接可以是通过一些接口,装置或单元的间接耦合或通讯连接,可以是电性,机械或其它的形式。In the embodiments provided by the present application, it should be understood that the disclosed apparatus/electronic device and method may be implemented in other manners. For example, the device/electronic device embodiments described above are merely illustrative. For example, the division of the above modules or units is only a logical functional division, and may be implemented in another manner, such as multiple units or Components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the mutual coupling or direct coupling or communication connection shown or discussed may be an indirect coupling or communication connection through some interface, device or unit, and may be in electrical, mechanical or other form.
上述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。The units described above as separate components may or may not be physically separated. The components displayed as units may or may not be physical units, that is, may be located in one place, or may be distributed to multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of the embodiment.
另外,在本申请各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各 个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。In addition, each functional unit in each embodiment of the present application may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit. The above integrated unit can be implemented in the form of hardware or in the form of a software functional unit.
上述集成的模块/单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本申请实现上述实施例方法中的全部或部分流程,也可以通过计算机程序来指令相关的硬件来完成,上述的计算机程序可存储于一计算机可读存储介质中,该计算机程序在被处理器执行时,可实现上述各个方法实施例的步骤。其中,上述计算机程序包括计算机程序代码,上述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。上述计算机可读介质可以包括:能够携带上述计算机程序代码的任何实体或装置、记录介质、U盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、电载波信号、电信信号以及软件分发介质等。需要说明的是,上述计算机可读介质包含的内容可以根据司法管辖区内立法和专利实践的要求进行适当的增减,例如在某些司法管辖区,根据立法和专利实践,计算机可读介质不包括电载波信号和电信信号。The above-described integrated modules/units, if implemented in the form of software functional units and sold or used as stand-alone products, may be stored in a computer readable storage medium. Based on such understanding, the present application implements all or part of the processes in the above embodiments, and may also be completed by a computer program to instruct related hardware. The computer program may be stored in a computer readable storage medium. The steps of the various method embodiments described above may be implemented when executed by a processor. Wherein, the above computer program comprises computer program code, and the computer program code may be in the form of source code, object code form, executable file or some intermediate form. The computer readable medium may include any entity or device capable of carrying the above computer program code, a recording medium, a USB flash drive, a removable hard disk, a magnetic disk, an optical disk, a computer memory, a Read-Only Memory (ROM), a random Access memory (RAM, Random Access Memory), electrical carrier signals, telecommunications signals, and software distribution media. It should be noted that the contents of the above computer readable medium may be appropriately increased or decreased according to the requirements of legislation and patent practice in jurisdictions. For example, in some jurisdictions, according to legislation and patent practice, computer readable media are not Includes electrical carrier signals and telecommunications signals.
以上上述实施例仅用以说明本申请的技术方案,而非对其限制;尽管参照前述实施例对本申请进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围,均应包含在本申请的保护范围之内。The above embodiments are only used to explain the technical solutions of the present application, and are not limited thereto; although the present application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that the foregoing embodiments can still be The technical solutions are modified, or some of the technical features are replaced by equivalents; and the modifications or substitutions do not deviate from the spirit and scope of the technical solutions of the embodiments of the present application, and should be included in the present disclosure. Within the scope of protection of the application.

Claims (10)

  1. 一种基于磁共振图像的脑龄测试方法,其特征在于,包括:A brain age testing method based on magnetic resonance images, characterized in that it comprises:
    获取待测个体的T1加权大脑磁共振图像;Obtaining a T1-weighted brain magnetic resonance image of the individual to be tested;
    基于所述T1加权大脑磁共振图像确定所述待测个体的脑结构、脑叶以及脑组织分割图谱;Determining a brain structure, a brain lobe, and a brain tissue segmentation map of the test subject based on the T1-weighted brain magnetic resonance image;
    计算所述脑结构的归一化体积,并基于所述脑组织分割图谱计算所述脑叶的脑萎缩值;Calculating a normalized volume of the brain structure, and calculating a brain atrophy value of the brain lobe based on the brain tissue segmentation map;
    将所述归一化体积以及所述脑萎缩值输入脑龄估计模型,获取所述待测个体的脑龄,其中,所述脑龄估计模型根据训练样本个体的归一化体积、脑萎缩值以及脑龄训练得到。Entering the normalized volume and the brain atrophy value into a brain age estimation model to obtain a brain age of the test subject, wherein the brain age estimation model is based on a normalized volume of the training sample individual and a brain atrophy value And brain age training is available.
  2. 如权利要求1所述的脑龄测试方法,其特征在于,所述基于所述T1加权大脑磁共振图像确定所述待测个体的脑结构、脑叶以及脑组织分割图谱,包括:The cerebral age testing method according to claim 1, wherein the determining a brain structure, a brain lobe, and a brain tissue segmentation map of the test subject based on the T1-weighted brain magnetic resonance image comprises:
    将所述T1加权大脑磁共振图像与预设的脑模板库进行配准,获得所述T1加权大脑磁共振图像所对应的脑结构、脑叶以及脑组织概率图谱;Registering the T1-weighted brain magnetic resonance image with a preset brain template library to obtain a brain structure, a brain lobe, and a brain tissue probability map corresponding to the T1-weighted brain magnetic resonance image;
    将所述T1加权大脑磁共振图像所对应的脑组织概率图谱作为先验知识,对所述T1加权大脑磁共振图像进行脑组织分割,得到所述T1加权大脑磁共振图像所对应的脑组织分割图谱;Using the brain tissue probability map corresponding to the T1-weighted brain magnetic resonance image as a priori knowledge, the brain tissue segmentation is performed on the T1-weighted brain magnetic resonance image, and the brain tissue segmentation corresponding to the T1-weighted brain magnetic resonance image is obtained. Map
    其中,所述脑模板库中包括:两个以上不同大脑的T1加权大脑磁共振图像,以及分别与各个T1加权大脑磁共振图像对应的脑结构、脑叶以及脑组织概率图谱。The brain template library includes: a T1-weighted brain magnetic resonance image of two or more different brains, and a brain structure, a brain lobe, and a brain tissue probability map corresponding to each T1-weighted brain magnetic resonance image, respectively.
  3. 如权利要求1所述的脑龄测试方法,其特征在于,所述基于所述脑组织分割图谱计算所述脑叶的脑萎缩值,包括:The brain age testing method according to claim 1, wherein the calculating a brain atrophy value of the brain lobe based on the brain tissue segmentation map comprises:
    基于所述脑组织分割图谱,确定所述脑叶中的脑灰质体积、脑白质体积以及脑脊液体积;Determining a gray matter volume, a white matter volume, and a cerebrospinal fluid volume in the brain lobe based on the brain tissue segmentation map;
    根据脑萎缩计算公式计算所述脑叶的脑萎缩值,其中,所述脑萎缩计算公式为:Calculating a brain atrophy value of the brain lobe according to a brain atrophy calculation formula, wherein the brain atrophy calculation formula is:
    Figure PCTCN2018072603-appb-100001
    Figure PCTCN2018072603-appb-100001
  4. 如权利要求1所述的脑龄测试方法,其特征在于,所述脑龄测试方法还包括:The brain age testing method according to claim 1, wherein the brain age testing method further comprises:
    将训练样本个体的归一化体积以及脑萎缩值作为自变量,以训练样本个体的脑龄作为因变量,应用线性支持向量机,建立所述脑龄估计模型。The normalized volume of the training sample individual and the brain atrophy value are used as independent variables, and the brain age of the training sample individual is used as the dependent variable, and the brain age estimation model is established by applying a linear support vector machine.
  5. 如权利要求1至4中任一项所述的脑龄测试方法,其特征在于,所述基于所述T1加权大脑磁共振图像确定所述待测个体的脑结构、脑叶以及脑组织分割图谱之前,还包括:The brain age testing method according to any one of claims 1 to 4, wherein the determining a brain structure, a brain lobe, and a brain tissue segmentation map of the test subject based on the T1-weighted brain magnetic resonance image Previously, it also included:
    对所述待测个体的T1加权大脑磁共振图像进行预处理,所述预处理包括如下一项或两项以上:降噪、去偏场、像素范围归一化;Performing preprocessing on the T1-weighted brain magnetic resonance image of the individual to be tested, the pre-processing including one or more of the following: noise reduction, de-biasing field, pixel range normalization;
    所述基于所述T1加权大脑磁共振图像确定所述待测个体的脑结构、脑叶以及脑组织分割图谱为:Determining, according to the T1-weighted brain magnetic resonance image, a brain structure, a brain lobe, and a brain tissue segmentation map of the test subject:
    基于所述预处理后得到的T1加权大脑磁共振图像确定所述待测个体的脑结构、脑叶以及脑组织分割图谱。A brain structure, a brain lobe, and a brain tissue segmentation map of the test subject are determined based on the T1-weighted brain magnetic resonance image obtained after the pre-processing.
  6. 一种基于磁共振图像的脑龄测试装置,其特征在于,包括:A brain age testing device based on magnetic resonance images, characterized in that it comprises:
    图像获取单元,用于获取待测个体的T1加权大脑磁共振图像;An image acquisition unit, configured to acquire a T1-weighted brain magnetic resonance image of the individual to be tested;
    图像分析单元,用于基于所述T1加权大脑磁共振图像确定所述待测个体的脑结构、脑叶以及脑组织分割图谱;An image analyzing unit, configured to determine a brain structure, a brain lobe, and a brain tissue segmentation map of the test subject based on the T1-weighted brain magnetic resonance image;
    参数计算单元,用于计算所述脑结构的归一化体积,并基于所述脑组织分割图谱计算所述脑叶的脑萎缩值;a parameter calculation unit, configured to calculate a normalized volume of the brain structure, and calculate a brain atrophy value of the brain lobe based on the brain tissue segmentation map;
    脑龄获取单元,用于将所述归一化体积以及所述脑萎缩值输入脑龄估计模型,获取所述待测个体的脑龄,其中,所述脑龄估计模型根据训练样本个体的归一化体积、脑萎缩值以及脑龄训练得到。a brain age acquisition unit, configured to input the normalized volume and the brain atrophy value into a brain age estimation model, and obtain a brain age of the test subject, wherein the brain age estimation model is based on the individuality of the training sample individual A volume, brain atrophy value and brain age training are obtained.
  7. 如权利要求6所述的脑龄测试装置,其特征在于,所述图像分析单元具体用于:The cerebral age testing device according to claim 6, wherein the image analyzing unit is specifically configured to:
    将所述T1加权大脑磁共振图像与预设的脑模板库进行配准,获得所述T1加权大脑磁共振图像所对应的脑结构、脑叶以及脑组织概率图谱;Registering the T1-weighted brain magnetic resonance image with a preset brain template library to obtain a brain structure, a brain lobe, and a brain tissue probability map corresponding to the T1-weighted brain magnetic resonance image;
    将所述T1加权大脑磁共振图像所对应的脑组织概率图谱作为先验知识,对所述T1加权大脑磁共振图像进行脑组织分割,得到所述T1加权大脑磁共振图像所对应的脑组织分割图谱;Using the brain tissue probability map corresponding to the T1-weighted brain magnetic resonance image as a priori knowledge, the brain tissue segmentation is performed on the T1-weighted brain magnetic resonance image, and the brain tissue segmentation corresponding to the T1-weighted brain magnetic resonance image is obtained. Map
    其中,所述脑模板库中包括:两个以上不同大脑的T1加权大脑磁共振图像,以及分别与各个T1加权大脑磁共振图像对应的脑结构、脑叶以及脑组织概率图谱。The brain template library includes: a T1-weighted brain magnetic resonance image of two or more different brains, and a brain structure, a brain lobe, and a brain tissue probability map corresponding to each T1-weighted brain magnetic resonance image, respectively.
  8. 如权利要求6所述的脑龄测试装置,其特征在于,所述参数计算单元具体用于:The cerebral age testing device according to claim 6, wherein the parameter calculating unit is specifically configured to:
    计算所述脑结构的归一化体积,并基于所述脑组织分割图谱,确定所述脑叶中的脑灰质体积、脑白质体积以及脑脊液体积;Calculating a normalized volume of the brain structure, and determining a gray matter volume, a white matter volume, and a cerebrospinal fluid volume in the brain lobe based on the brain tissue segmentation map;
    根据脑萎缩计算公式计算所述脑叶的脑萎缩值,其中,所述脑萎缩计算公式为:Calculating a brain atrophy value of the brain lobe according to a brain atrophy calculation formula, wherein the brain atrophy calculation formula is:
    Figure PCTCN2018072603-appb-100002
    Figure PCTCN2018072603-appb-100002
  9. 一种电子设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现如权利要求1至5任一项所述方法的步骤。An electronic device comprising a memory, a processor, and a computer program stored in the memory and operable on the processor, wherein the processor executes the computer program as claimed in claim 1 5 The steps of any of the methods described.
  10. 一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1至5任一项所述方法的步骤。A computer readable storage medium storing a computer program, wherein the computer program, when executed by a processor, implements the steps of the method of any one of claims 1 to 5.
PCT/CN2018/072603 2018-01-15 2018-01-15 Magnetic resonance image based brain age test method and apparatus WO2019136745A1 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
PCT/CN2018/072603 WO2019136745A1 (en) 2018-01-15 2018-01-15 Magnetic resonance image based brain age test method and apparatus
CN201880000008.4A CN110337670B (en) 2018-01-15 2018-01-15 Brain age testing method and brain age testing device based on magnetic resonance image

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2018/072603 WO2019136745A1 (en) 2018-01-15 2018-01-15 Magnetic resonance image based brain age test method and apparatus

Publications (1)

Publication Number Publication Date
WO2019136745A1 true WO2019136745A1 (en) 2019-07-18

Family

ID=67219253

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2018/072603 WO2019136745A1 (en) 2018-01-15 2018-01-15 Magnetic resonance image based brain age test method and apparatus

Country Status (2)

Country Link
CN (1) CN110337670B (en)
WO (1) WO2019136745A1 (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111134677A (en) * 2020-02-14 2020-05-12 北京航空航天大学 Magnetic resonance image-based stroke-induced disability prediction method and system
CN111415361A (en) * 2020-03-31 2020-07-14 浙江大学 Method and device for estimating brain age of fetus and detecting abnormality based on deep learning
EP3848941A1 (en) * 2020-01-10 2021-07-14 Acer Incorporated Model training method and electronic device
CN113221927A (en) * 2020-01-21 2021-08-06 宏碁股份有限公司 Model training method and electronic device
CN117036793A (en) * 2023-07-31 2023-11-10 复旦大学 Brain age assessment method and device based on multi-scale features of PET (positron emission tomography) images
EP4290453A1 (en) * 2022-06-07 2023-12-13 Gil Medical Center System and method for segmentating brain in mri images

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130202173A1 (en) * 2008-02-19 2013-08-08 vascuVis Inc. Classification of biological tissue by multi-mode data registration, segmentation and characterization
CN105046709A (en) * 2015-07-14 2015-11-11 华南理工大学 Nuclear magnetic resonance imaging based brain age analysis method
CN105212936A (en) * 2015-09-09 2016-01-06 首都医科大学宣武医院 A kind of generation method of brain template
CN105816192A (en) * 2016-03-03 2016-08-03 王雪原 Method for three-dimensional registration and brain tissue extraction of individual human brain multimodality medical images
CN107103612A (en) * 2017-03-28 2017-08-29 深圳博脑医疗科技有限公司 Automate the quantitative calculation method of subregion brain atrophy

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007114238A1 (en) * 2006-03-30 2007-10-11 National University Corporation Shizuoka University Apparatus for determining brain atrophy, method of determining brain atrophy and program for determining brain atrophy
CN107507212B (en) * 2017-08-18 2020-08-18 中国科学院深圳先进技术研究院 Digital brain visualization method and device, computing equipment and storage medium

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130202173A1 (en) * 2008-02-19 2013-08-08 vascuVis Inc. Classification of biological tissue by multi-mode data registration, segmentation and characterization
CN105046709A (en) * 2015-07-14 2015-11-11 华南理工大学 Nuclear magnetic resonance imaging based brain age analysis method
CN105212936A (en) * 2015-09-09 2016-01-06 首都医科大学宣武医院 A kind of generation method of brain template
CN105816192A (en) * 2016-03-03 2016-08-03 王雪原 Method for three-dimensional registration and brain tissue extraction of individual human brain multimodality medical images
CN107103612A (en) * 2017-03-28 2017-08-29 深圳博脑医疗科技有限公司 Automate the quantitative calculation method of subregion brain atrophy

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3848941A1 (en) * 2020-01-10 2021-07-14 Acer Incorporated Model training method and electronic device
US11455498B2 (en) 2020-01-10 2022-09-27 Acer Incorporated Model training method and electronic device
CN113221927A (en) * 2020-01-21 2021-08-06 宏碁股份有限公司 Model training method and electronic device
CN113221927B (en) * 2020-01-21 2024-01-23 宏碁股份有限公司 Model training method and electronic device
CN111134677A (en) * 2020-02-14 2020-05-12 北京航空航天大学 Magnetic resonance image-based stroke-induced disability prediction method and system
CN111134677B (en) * 2020-02-14 2021-05-04 北京航空航天大学 Magnetic resonance image-based stroke-induced disability prediction method and system
CN111415361A (en) * 2020-03-31 2020-07-14 浙江大学 Method and device for estimating brain age of fetus and detecting abnormality based on deep learning
EP4290453A1 (en) * 2022-06-07 2023-12-13 Gil Medical Center System and method for segmentating brain in mri images
CN117036793A (en) * 2023-07-31 2023-11-10 复旦大学 Brain age assessment method and device based on multi-scale features of PET (positron emission tomography) images
CN117036793B (en) * 2023-07-31 2024-04-19 复旦大学 Brain age assessment method and device based on multi-scale features of PET (positron emission tomography) images

Also Published As

Publication number Publication date
CN110337670A (en) 2019-10-15
CN110337670B (en) 2023-09-29

Similar Documents

Publication Publication Date Title
WO2019136745A1 (en) Magnetic resonance image based brain age test method and apparatus
Manjón et al. volBrain: an online MRI brain volumetry system
Giraud et al. An optimized patchmatch for multi-scale and multi-feature label fusion
Cardoso et al. STEPS: Similarity and Truth Estimation for Propagated Segmentations and its application to hippocampal segmentation and brain parcelation
Wolz et al. LEAP: learning embeddings for atlas propagation
Hao et al. Local label learning (LLL) for subcortical structure segmentation: application to hippocampus segmentation
AU2016308097B2 (en) Three-dimensional (3D) convolution with 3D batch normalization
Das et al. Registration based cortical thickness measurement
Mendrik et al. MRBrainS challenge: online evaluation framework for brain image segmentation in 3T MRI scans
Ji et al. Robust spatially constrained fuzzy c-means algorithm for brain MR image segmentation
WO2018176985A1 (en) Quantitative calculation method for level of brain atrophy based on automatic segmentation
Huang et al. Brain extraction based on locally linear representation-based classification
Zhou et al. Segmentation of subcortical brain structures using fuzzy templates
Han et al. Brain extraction from normal and pathological images: a joint PCA/image-reconstruction approach
Nigro et al. Fully automated segmentation of the pons and midbrain using human T1 MR brain images
Bilgel et al. Temporal trajectory and progression score estimation from voxelwise longitudinal imaging measures: Application to amyloid imaging
Ye et al. Segmentation of the cerebellar peduncles using a random forest classifier and a multi-object geometric deformable model: application to spinocerebellar ataxia type 6
Amiri et al. 3D cerebral MR image segmentation using multiple-classifier system
Hu et al. Nonlocal regularization for active appearance model: Application to medial temporal lobe segmentation
CN114708263B (en) Individual brain functional region positioning method, device, equipment and storage medium
Gaweł et al. Automatic spine tissue segmentation from MRI data based on cascade of boosted classifiers and active appearance model
Safavian et al. An automatic level set method for hippocampus segmentation in MR images
Pappas et al. Improved normalization of lesioned brains via cohort‐specific templates
Pujadas et al. Shape-based normalized cuts using spectral relaxation for biomedical segmentation
Almodóvar-Rivera et al. FAST adaptive smoothing and thresholding for improved activation detection in low-signal fMRI

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 18899116

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 18899116

Country of ref document: EP

Kind code of ref document: A1