Disclosure of Invention
In order to solve the defects of the prior art, a first object of the present invention is to provide a functional magnetic resonance imaging brain partitioning method based on a quantum potential energy model, which can quickly and accurately partition fMRI brain images containing large data volume.
The invention discloses a functional magnetic resonance imaging brain partitioning method based on a quantum potential energy model, which comprises the following steps:
preprocessing the fMRI brain image;
and partitioning the preprocessed fMRI brain image by using a quantum potential energy model.
Further, the specific process of partitioning the preprocessed fMRI brain image by using the quantum potential energy model is as follows:
extracting a gray matter part in the preprocessed fMRI brain image to obtain a gray matter image;
converting the gray matter image into a two-dimensional matrix form M; the number of the rows is equal to the total prime number of the grey matter images, the number of the columns is an integral multiple of 4, and the multiple is equal to the number of fMRI brain images to be partitioned; each row represents a voxel, and each 4 columns represent the x-axis position coordinate, the y-axis position coordinate and the z-axis position coordinate of the voxel and the gray value of the voxel respectively;
the two-dimensional matrix M or the value of each column of the matrix formed after the principal component analysis of the two-dimensional matrix M is used as the input value of the quantum potential energy model to construct a quantum wave function, and further the voxel position coordinates and the gray value of the fMRI brain image are converted into the particle distribution of a quantum system;
inputting the constructed quantum wave function into a Schrodinger equation to obtain an expression of particle space potential energy of a quantum system, and further changing the process of searching for the center of a brain partition in a fMRI brain image into searching for a minimum value point of the potential energy in the quantum system;
calculating all minimum value points of the particle space potential energy of the quantum system as the center of the brain subarea;
setting a potential energy range, associating the voxels within the potential energy range from the center of the brain partition to the corresponding partition, and modifying the gray value of the voxels to be the label value corresponding to the center of the partition;
outputting a result matrix; the number of columns of the result matrix is integral multiple of 4, the multiple is equal to the number of fMRI brain images to be partitioned, and each 4 columns are respectively the x-axis position coordinate, the y-axis position coordinate and the z-axis position coordinate of a voxel and a corresponding label value; the number of rows is equal to the number of voxels, and each row corresponds to one voxel;
and projecting the corresponding voxels to a three-dimensional curved surface standard MNI brain model according to the x-axis, y-axis and z-axis position coordinates of the voxels in the result matrix, and displaying different colors by using different label values to obtain a brain partition result.
Further, the method further comprises: calling an OpenGL interface, and accurately presenting a brain partition result on a three-dimensional MNI standard brain model through three-dimensional rendering; and simultaneously displaying three views of an XZ plane, a YZ plane and an XY plane so as to obtain three views of a sagittal position, a coronal position and a transverse position of a brain partition result.
Further, the method further comprises: and outputting the partition result, wherein the partition result is output in an NIfTI format and is applied to fMRI functional connection and brain functional network analysis.
Further, the preprocessed fMRI brain image is subjected to dot multiplication with a gray matter template in a standard MNI space, and a gray matter part in the preprocessed fMRI brain image is extracted to obtain a gray matter image.
Further, the quantum wave function is:
wherein xi=[β1,β2,…,βn],β1,β2,…βnN principal component values respectively corresponding to the ith voxel, i.e. n elements of the ith row of the matrix M, x being an n-dimensional column vector space constructed from the input values, xiThe matrix is formed after principal component analysis is carried out on the two-dimensional matrix M; σ is the width parameter of the wave function, and the size of σ is related to the number of partitions.
Further, the method further comprises:
controlling the number of minimum value points by setting a threshold of the extreme value points, recording the number of the minimum value points after the threshold is set as l, and assigning the minimum value points to continuous label values 1,2 and … l from small to large according to x-axis, y-axis and z-axis position coordinates; wherein l is a positive integer greater than or equal to 1.
Further, the process of preprocessing the fMRI brain image includes:
correcting a temporal layer of the fMRI brain images;
correcting a spatial layer of the fMRI brain image;
registering individual fMRI brain images to an MNI spatial template;
performing spatial smoothing on the fMRI brain image registered to the MNI spatial template;
and filtering the fMRI brain image after the spatial smoothing processing.
The second purpose of the invention is to provide a functional magnetic resonance imaging brain partition method based on a quantum potential energy model, which is used for outputting partition results in an NIfTI format, and is applied to fMRI functional connection and brain functional network analysis, so that the accuracy of fMRI analysis is improved.
The invention also provides a functional magnetic resonance imaging brain partition system based on the quantum potential energy model.
The invention relates to a functional magnetic resonance imaging brain partition system based on a quantum potential energy model, which comprises a memory and a processor; the processor configured to perform the steps of:
preprocessing the fMRI brain image;
and partitioning the preprocessed fMRI brain image by using a quantum potential energy model.
Further, in the processor, the specific process of partitioning the preprocessed fMRI brain image by using the quantum potential energy model is as follows:
extracting a gray matter part in the preprocessed fMRI brain image to obtain a gray matter image;
converting the gray matter image into a two-dimensional matrix form M; the number of the rows is equal to the total prime number of the grey matter images, the number of the columns is an integral multiple of 4, and the multiple is equal to the number of fMRI brain images to be partitioned; each row represents a voxel, and each 4 columns represent the x-axis position coordinate, the y-axis position coordinate and the z-axis position coordinate of the voxel and the gray value of the voxel respectively;
the two-dimensional matrix M or the value of each column of the matrix formed after the principal component analysis of the two-dimensional matrix M is used as the input value of the quantum potential energy model to construct a quantum wave function, and further the voxel position coordinates and the gray value of the fMRI brain image are converted into the particle distribution of a quantum system;
inputting the constructed quantum wave function into a Schrodinger equation to obtain an expression of particle space potential energy of a quantum system, and further changing the process of searching for the center of a brain partition in a fMRI brain image into searching for a minimum value point of the potential energy in the quantum system;
calculating all minimum value points of the particle space potential energy of the quantum system as the center of the brain subarea;
setting a potential energy range, associating the voxels within the potential energy range from the center of the brain partition to the corresponding partition, and modifying the gray value of the voxels to be the label value corresponding to the center of the partition;
outputting a result matrix; the number of columns of the result matrix is integral multiple of 4, the multiple is equal to the number of fMRI brain images to be partitioned, and each 4 columns are respectively the x-axis position coordinate, the y-axis position coordinate and the z-axis position coordinate of a voxel and a corresponding label value; the number of rows is equal to the number of voxels, and each row corresponds to one voxel;
and projecting the corresponding voxels to a three-dimensional curved surface standard MNI brain model according to the x-axis, y-axis and z-axis position coordinates of the voxels in the result matrix, and displaying different colors by using different label values to obtain a brain partition result.
Further, the processor is configured to perform:
and outputting the partition result, wherein the partition result is output in an NIfTI format and is applied to the subsequent analysis of fMRI functional connection and brain functional network.
Further, in the processor, the preprocessing the fMRI brain image includes:
correcting a temporal layer of the fMRI brain images;
correcting a spatial layer of the fMRI brain image;
registering individual fMRI brain images to an MNI spatial template;
performing spatial smoothing on the fMRI brain image registered to the MNI spatial template;
and filtering the fMRI brain image after the spatial smoothing processing.
Compared with the prior art, the invention has the beneficial effects that:
according to the functional magnetic resonance imaging brain partitioning method and system based on the quantum potential energy model, advantages of a quantum algorithm in analysis and processing of large data and high-dimensional data are utilized, the brain partitioning is converted into the quantum potential energy model, and clustering partitioning is performed on single tested data or multiple tested fMRI data through quantum calculation.
Detailed Description
It should be noted that the following detailed description is exemplary and is intended to provide further explanation of the disclosure. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of example embodiments according to the present application. As used herein, the singular forms "a", "an" and "the" are intended to include the plural forms as well, and it should be understood that when the terms "comprises" and/or "comprising" are used in this specification, they specify the presence of stated features, steps, operations, devices, components, and/or combinations thereof, unless the context clearly indicates otherwise.
FIG. 1 is a flow chart of a functional magnetic resonance imaging brain partitioning method based on a quantum potential energy model.
As shown in fig. 1, a functional magnetic resonance imaging brain partition method based on a quantum potential energy model of the present invention includes:
step 1: the fMRI brain images are preprocessed.
The fMRI brain image input into the fMRI brain partition system based on the quantum potential energy model needs to be preprocessed. The fMRI brain image is a magnetic resonance brain image acquired by a fast imaging sequence such as a gradient echo planar imaging sequence and a spin echo sequence by using a blood oxygen level dependent technique.
Preprocessing employs Statistical Parametric Mapping (SPM) software commonly used for functional magnetic resonance image processing. In a particular embodiment, the process of pre-processing the fMRI brain image includes:
step 1.1: correcting a temporal layer of the fMRI brain images;
before correcting the time layer of the fMRI brain image, the method further comprises the following steps: and (3) converting the data format, namely converting the original DICOM format data into NIfTI format data commonly used by a computer.
Specifically, a Slice Timing module of the SPM is used to perform temporal layer correction on the data.
Step 1.2: the spatial layer of the fMRI brain image is corrected.
Specifically, the Realign module of the SPM is adopted to correct the spatial layer of the data.
Step 1.3: individual fMRI brain images are registered to the MNI spatial template.
Specifically, by using the normalation module of the SPM, the individual brain fMRI data can be registered to the Montreal Nerve Institute (MNI) standard brain space by three methods, namely, an EPI template, a structural image joint registration, or a differential homoembryo image fusion algorithm.
Step 1.4: and performing spatial smoothing on the fMRI brain image registered to the MNI spatial template.
Specifically, a smoothening module of the SPM is adopted, and a gaussian kernel of full width at half maximum can be selected to perform convolution processing on the image, and the full width at half maximum is generally 6-8 mm.
Step 1.5: and filtering the fMRI brain image after the spatial smoothing processing.
Specifically, the image data may be filtered using a band-selectable low-pass filter.
Step 2: and partitioning the preprocessed fMRI brain image by using a quantum potential energy model.
In a specific embodiment, as shown in fig. 2, a specific process of partitioning the preprocessed fMRI brain image by using the quantum potential energy model is as follows:
step 2.1: extracting a gray matter part in the preprocessed fMRI brain image to obtain a gray matter image;
specifically, the preprocessed fMRI brain image is dot-multiplied with a gray matter template in a standard MNI space, and a gray matter portion in the preprocessed fMRI brain image is extracted to obtain a gray matter image, as shown in fig. 3.
For example: and performing point multiplication on each preprocessed fMRI image to be tested and a binary gray matter template in a standard MNI space, extracting a gray matter part in the image, performing brain partition processing by adopting the gray matter image, and performing first-step gray matter extraction on the residual 67531 voxels in the image.
Step 2.2: converting the gray matter image into a two-dimensional matrix form M; the number of the rows is equal to the total prime number of the grey matter images, the number of the columns is an integral multiple of 4, and the multiple is equal to the number of fMRI brain images to be partitioned; each row represents a voxel, and each 4 columns represent the x-axis position coordinate, the y-axis position coordinate and the z-axis position coordinate of the voxel and the gray value of the voxel respectively.
Specifically, the three-dimensional gray matter image is sequentially converted into a two-dimensional matrix form of 67531 rows and 4m columns according to the x-axis position coordinates, y-axis position coordinates and the gray value of each voxel, wherein m represents the number of fMRI images, and when only one tested fMRI image exists, m is 1. Each row corresponds to a voxel, each 4 columns correspond to the x-axis position coordinate, the y-axis position coordinate and the z-axis position coordinate of the voxel and the gray value of the voxel, and the 67531 row-4M-column two-dimensional matrix is marked as M.
Step 2.3: and (3) taking the value of each column of the two-dimensional matrix M or the matrix formed after the principal component analysis of the two-dimensional matrix M as an input value of the quantum potential energy model, constructing a quantum wave function, and further converting the voxel position coordinates and the gray value of the fMRI brain image into particle distribution of a quantum system.
When the input data is a plurality of pieces of fMRI data to be tested (m >1), principal component analysis may be employed, the main steps of which are:
(1) normalization: calculating the mean value of each column of the matrix M, subtracting the mean value corresponding to the column from each element of the column, and normalizing, wherein the normalized matrix is marked as M';
(2) calculating a covariance matrix C of the matrix M', calculating an eigenvalue lambda of the covariance matrix CiAnd orthogonalizing the eigenvector ai;
(3) The characteristic value lambda is measurediArranged from large to small such that λ1>λ2>…λiAnd the orthogonalized eigenvectors corresponding to the eigenvalues are also arranged according to the order of the eigenvalues, so that the ith principal component corresponding to the original matrix M is Fi=ai'M'. N main components of the two-dimensional matrix M can be extracted as required for further analysis, the matrix output in the third step is F, and F has 67531 rows and n columns; when m is 1, the principal component analysis may not be performed, and this step may be skipped.
Specifically, the quantum wave function is:
wherein xi=[β1,β2,…,βn],β1,β2,…βnOf n principal components corresponding to ith voxelThe value, i.e. n elements of the ith row of the matrix M, x is an n-dimensional column vector space constructed from the input values, and if not specified, x is defaultediThe matrix is formed after principal component analysis is carried out on the two-dimensional matrix M; σ is the width parameter of the wave function, and the size of σ is related to the number of partitions.
Step 2.4: and inputting the constructed quantum wave function into a Schrodinger equation to obtain an expression of the particle space potential energy of the quantum system, and further changing the process of searching the center of the brain partition in the fMRI brain image into searching a minimum value point of the potential energy in the quantum system.
The form of Schrodinger equation is shown in formula (1)
Wherein H is a Hamiltonian, sigma is a width parameter of the wave function,
for Laplace operator, V (x) is space potential energy, psi is wave function describing quantum system, E is energy eigenvalue, E ═ d/2 is minimum possible eigenvalue, d is input data x
iOf (c) is calculated. Constructing the quantum wave function in the fourth step
Inputting the formula (1), and solving the formula (1) to obtain the expression of V (x)
According to the principle of quantum mechanics, the lowest point of potential energy, i.e. the point where the quantum system is relatively stable, is often the point where the wave function takes the maximum value, and can be used as the center of the brain partition. Therefore, the process of searching the brain partition center in the fMRI image is changed into the process of searching the minimum value point of the potential energy in the quantum system.
Step 2.5: calculating all minimum value points of the particle space potential energy of the quantum system as the center of the brain subarea;
specifically, the number of minimum points is controlled by setting a threshold of the minimum points, the number of minimum points after the threshold is set is recorded as l, and the minimum points are sequentially assigned to continuous label values 1,2 and … l from small to large according to x-axis, y-axis and z-axis position coordinates; wherein l is a positive integer greater than or equal to 1.
As shown in fig. 4(a) -4 (d), wherein fig. 4(a) is a partial diagram of a quantum potential energy model, and the quantum potential energy model in the partial diagram has a total of 1,2, 3, 4, and 5 potential energy minimum points, which can be used as the brain partition center; fig. 4(b) is a potential energy threshold, by setting the potential energy threshold, the minimum value points 2 and 4 do not reach the threshold, and only three minimum value points 1, 3, and 5 are left as the brain partition center after threshold screening. Fig. 4(c) shows the xy plane when z is 50, i.e., the 50 th layer transverse position view, the outer black boxes indicate fMRI scan fields, black is blank voxels, and white indicates effective fMRI voxels; fig. 4(d) is an equipotential surface of a cross-sectional view of layer 50, where the box is the potential minimum point, i.e., the center of the partition, passing a set threshold. The lines with the same gray level in the graph are planes formed by points with the same potential energy, namely equipotential surfaces, and the points in a certain potential energy range in the center of the partition can be associated to the center of the partition by setting the potential energy range, namely the equipotential surfaces, so that the partition of each voxel of the whole brain is completed.
Step 2.6: setting a potential energy range, associating the voxels within the potential energy range from the center of the brain partition to the corresponding partition, and modifying the gray value of the voxels to be the label value corresponding to the center of the partition;
step 2.7: outputting a result matrix; the number of columns of the result matrix is integral multiple of 4, the multiple is equal to the number of fMRI brain images to be partitioned, and each 4 columns are respectively the x-axis position coordinate, the y-axis position coordinate and the z-axis position coordinate of a voxel and a corresponding label value; the number of rows is equal to the number of voxels, and each row corresponds to one voxel;
step 2.8: and projecting the corresponding voxels to a three-dimensional curved surface standard MNI brain model according to the x-axis, y-axis and z-axis position coordinates of the voxels in the result matrix, and displaying different colors by using different label values to obtain a brain partition result.
In a specific implementation, the method further comprises:
step 2.9: and storing the result matrix as an NIfTI format file for output, and applying the NIfTI format file to subsequent analysis of fMRI functional connection and brain functional network.
In another embodiment, the method further comprises: calling an OpenGL interface, and accurately presenting a brain partition result on a three-dimensional MNI standard brain model through three-dimensional rendering; simultaneously displaying three views of an XZ plane, a YZ plane and an XY plane, and further obtaining three views of a sagittal position, a coronal position and a transverse position of a brain partition result, as shown in fig. 5(a) -5 (d), fig. 5(a) is a three-dimensional brain partition result graph; FIG. 5(b) is a cross-sectional view angle result chart; FIG. 5(c) is a sagittal view angle results plot; FIG. 5(d) is a coronary view angle result chart.
The functional magnetic resonance imaging brain partitioning method based on the quantum potential energy model utilizes the advantages of a quantum algorithm in analyzing and processing large data and high-dimensional data to convert brain partitioning into the quantum potential energy model, and clustering and partitioning are carried out on single tested fMRI data or a plurality of tested fMRI data through quantum calculation.
The invention also provides a functional magnetic resonance imaging brain partition system based on the quantum potential energy model.
The invention relates to a functional magnetic resonance imaging brain partition system based on a quantum potential energy model, which comprises a memory and a processor; the processor configured to perform the steps of:
(1) preprocessing the fMRI brain image;
specifically, the process of preprocessing the fMRI brain image includes:
the fMRI brain image input into the fMRI brain partition system based on the quantum potential energy model needs to be preprocessed. The fMRI brain image is a magnetic resonance brain image acquired by a fast imaging sequence such as a gradient echo planar imaging sequence and a spin echo sequence by using a blood oxygen level dependent technique.
Preprocessing employs Statistical Parametric Mapping (SPM) software commonly used for functional magnetic resonance image processing. In a particular embodiment, the process of pre-processing the fMRI brain image includes:
step 1.1: correcting a temporal layer of the fMRI brain images;
before correcting the time layer of the fMRI brain image, the method further comprises the following steps: and (3) converting the data format, namely converting the original DICOM format data into NIfTI format data commonly used by a computer.
Specifically, a Slice Timing module of the SPM is used to perform temporal layer correction on the data.
Step 1.2: the spatial layer of the fMRI brain image is corrected.
Specifically, the Realign module of the SPM is adopted to correct the spatial layer of the data.
Step 1.3: individual fMRI brain images are registered to the MNI spatial template.
Specifically, by using the normalation module of the SPM, the individual brain fMRI data can be registered to the Montreal Nerve Institute (MNI) standard brain space by three methods, namely, an EPI template, a structural image joint registration, or a differential homoembryo image fusion algorithm.
Step 1.4: and performing spatial smoothing on the fMRI brain image registered to the MNI spatial template.
Specifically, a smoothening module of the SPM is adopted, and a Gaussian kernel with the full width at half maximum can be selected to perform spatial processing on the image, wherein the full width at half maximum is 6-8 mm generally.
Step 1.5: and filtering the fMRI brain image after the spatial smoothing processing.
Specifically, the image data may be filtered using a band-selectable low-pass filter.
(2) And partitioning the preprocessed fMRI brain image by using a quantum potential energy model.
Specifically, as shown in fig. 2, the specific process of partitioning the preprocessed fMRI brain image by using the quantum potential energy model is as follows:
in a specific embodiment, as shown in fig. 2, a specific process of partitioning the preprocessed fMRI brain image by using the quantum potential energy model is as follows:
step 2.1: extracting a gray matter part in the preprocessed fMRI brain image to obtain a gray matter image;
specifically, the preprocessed fMRI brain image is dot-multiplied with a gray matter template in a standard MNI space, and a gray matter portion in the preprocessed fMRI brain image is extracted to obtain a gray matter image, as shown in fig. 3.
For example: and performing point multiplication on each preprocessed fMRI image to be tested and a binary gray matter template in a standard MNI space, extracting a gray matter part in the image, performing brain partition processing by adopting the gray matter image, and performing first-step gray matter extraction on the residual 67531 voxels in the image.
Step 2.2: converting the gray matter image into a two-dimensional matrix form M; the number of the rows is equal to the total prime number of the grey matter images, the number of the columns is an integral multiple of 4, and the multiple is equal to the number of fMRI brain images to be partitioned; each row represents a voxel, and each 4 columns represent the x-axis position coordinate, the y-axis position coordinate and the z-axis position coordinate of the voxel and the gray value of the voxel respectively.
Specifically, the three-dimensional gray matter image is sequentially converted into a two-dimensional matrix form of 67531 rows and 4m columns according to the x-axis position coordinates, y-axis position coordinates and the gray value of each voxel, wherein m represents the number of fMRI images, and when only one tested fMRI image exists, m is 1. Each row corresponds to a voxel, each 4 columns correspond to the x-axis position coordinate, the y-axis position coordinate and the z-axis position coordinate of the voxel and the gray value of the voxel, and the 67531 row-4M-column two-dimensional matrix is marked as M.
Step 2.3: and (3) taking the value of each column of the two-dimensional matrix M or the matrix formed after the principal component analysis of the two-dimensional matrix M as an input value of the quantum potential energy model, constructing a quantum wave function, and further converting the voxel position coordinates and the gray value of the fMRI brain image into particle distribution of a quantum system.
When the input data is a plurality of pieces of fMRI data to be tested (m >1), principal component analysis may be employed, the main steps of which are:
(1) normalization: calculating the mean value of each column of the matrix M, subtracting the mean value corresponding to the column from each element of the column, and normalizing, wherein the normalized matrix is marked as M';
(2) calculating a covariance matrix C of the matrix M', calculating an eigenvalue lambda of the covariance matrix CiAnd orthogonalizing the eigenvector ai;
(3) The characteristic value lambda is measurediArranged from large to small such that λ1>λ2>…λiAnd the orthogonalized eigenvectors corresponding to the eigenvalues are also arranged according to the order of the eigenvalues, so that the ith principal component corresponding to the original matrix M is Fi=ai'M'. N main components of the two-dimensional matrix M can be extracted as required to be analyzed in the next step, the matrix output in the third step is recorded as F, and the F has 67531 rows and n columns; when m is 1, the principal component analysis may not be performed, and this step may be skipped.
Specifically, the quantum wave function is:
wherein xi=[β1,β2,…,βn],β1,β2,…βnN principal component values respectively corresponding to the ith voxel, i.e. n elements in the ith row of the matrix M, x is an n-dimensional column vector space constructed according to input values, and if no special description is given, x is defaultediA matrix F formed after principal component analysis is carried out on the two-dimensional matrix M; σ is the width parameter of the wave function, and the size of σ is related to the number of partitions.
Step 2.4: and inputting the constructed quantum wave function into a Schrodinger equation to obtain an expression of the particle space potential energy of the quantum system, and further changing the process of searching the center of the brain partition in the fMRI brain image into searching a minimum value point of the potential energy in the quantum system.
The form of Schrodinger equation is shown in formula (1)
Wherein H is a Hamiltonian, sigma is a width parameter of the wave function,
is Laplace operator, V (x) is space potential energy, psi is wave function describing quantum system, E is energy eigenvalueD/2 is the smallest possible eigenvalue, d is the input data x
iOf (c) is calculated. Constructing the quantum wave function in the fourth step
Inputting the formula (1), and solving the formula (1) to obtain the expression of V (x)
According to the principle of quantum mechanics, the lowest point of potential energy, i.e. the point where the quantum system is relatively stable, is often the point where the wave function takes the maximum value, and can be used as the center of the brain partition. Therefore, the process of searching the brain partition center in the fMRI image is changed into the process of searching the minimum value point of the potential energy in the quantum system.
Step 2.5: calculating all minimum value points of the particle space potential energy of the quantum system as the center of the brain subarea;
specifically, the number of minimum points is controlled by setting a threshold of the minimum points, the number of minimum points after the threshold is set is recorded as l, and the minimum points are sequentially assigned to continuous label values 1,2 and … l from small to large according to x-axis, y-axis and z-axis position coordinates; wherein l is a positive integer greater than or equal to 1.
As shown in fig. 4(a) -4 (d), wherein fig. 4(a) is a partial diagram of a quantum potential energy model, and the quantum potential energy model in the partial diagram has a total of 1,2, 3, 4, and 5 potential energy minimum points, which can be used as the brain partition center; fig. 4(b) is a potential energy threshold, by setting the potential energy threshold, the minimum value points 2 and 4 do not reach the threshold, and only three minimum value points 1, 3, and 5 are left as the brain partition center after threshold screening. Fig. 4(c) shows the xy plane when z is 50, i.e., the 50 th layer transverse position view, the outer black boxes indicate fMRI scan fields, black is blank voxels, and white indicates effective fMRI voxels; fig. 4(d) is an equipotential surface of a cross-sectional view of layer 50, where the box is the potential minimum point, i.e., the center of the partition, passing a set threshold. The lines with the same gray level in the graph are planes formed by points with the same potential energy, namely equipotential surfaces, and the points in a certain potential energy range in the center of the partition can be associated to the center of the partition by setting the potential energy range, namely the equipotential surfaces, so that the partition of each voxel of the whole brain is completed.
Step 2.6: setting a potential energy range, associating the voxels within the potential energy range from the center of the brain partition to the corresponding partition, and modifying the gray value of the voxels to be the label value corresponding to the center of the partition;
step 2.7: outputting a result matrix; the number of columns of the result matrix is integral multiple of 4, the multiple is equal to the number of fMRI brain images to be partitioned, and each 4 columns are respectively the x-axis position coordinate, the y-axis position coordinate and the z-axis position coordinate of a voxel and a corresponding label value; the number of rows is equal to the number of voxels, and each row corresponds to one voxel;
step 2.8: and projecting the corresponding voxels to a three-dimensional curved surface standard MNI brain model according to the x-axis, y-axis and z-axis position coordinates of the voxels in the result matrix, and displaying different colors by using different label values to obtain a brain partition result.
The processor is further configured to:
and storing the result matrix as an NIfTI format file for output, and applying the NIfTI format file to subsequent analysis of fMRI functional connection and brain functional network.
In a further embodiment of the method according to the invention,
the processor is further configured to:
calling an OpenGL interface, and accurately presenting a brain partition result on a three-dimensional MNI standard brain model through three-dimensional rendering; simultaneously displaying three views of an XZ plane, a YZ plane and an XY plane, and further obtaining three views of a sagittal position, a coronal position and a transverse position of a brain partition result, as shown in fig. 5(a) -5 (d), fig. 5(a) is a three-dimensional brain partition result graph; FIG. 5(b) is a cross-sectional view angle result chart; FIG. 5(c) is a sagittal view angle results plot; FIG. 5(d) is a coronary view angle result chart.
The functional magnetic resonance imaging brain partition system based on the quantum potential energy model utilizes the advantages of a quantum algorithm in analyzing and processing large data and high-dimensional data to convert the brain partition into the quantum potential energy model, and conducts clustering partition on single tested data or a plurality of tested fMRI data through quantum calculation.
Although the embodiments of the present invention have been described with reference to the accompanying drawings, it is not intended to limit the scope of the present invention, and it should be understood by those skilled in the art that various modifications and variations can be made without inventive efforts by those skilled in the art based on the technical solution of the present invention.