CN110870770A - ICA-CNN classified fMRI space activation map smoothing and broadening method - Google Patents
ICA-CNN classified fMRI space activation map smoothing and broadening method Download PDFInfo
- Publication number
- CN110870770A CN110870770A CN201911144802.9A CN201911144802A CN110870770A CN 110870770 A CN110870770 A CN 110870770A CN 201911144802 A CN201911144802 A CN 201911144802A CN 110870770 A CN110870770 A CN 110870770A
- Authority
- CN
- China
- Prior art keywords
- cnn
- fmri
- space
- ica
- smoothing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
- A61B5/004—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part
- A61B5/0042—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part for the brain
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/40—Detecting, measuring or recording for evaluating the nervous system
- A61B5/4076—Diagnosing or monitoring particular conditions of the nervous system
- A61B5/4088—Diagnosing of monitoring cognitive diseases, e.g. Alzheimer, prion diseases or dementia
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
- A61B5/7267—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2134—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on separation criteria, e.g. independent component analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2135—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2576/00—Medical imaging apparatus involving image processing or analysis
- A61B2576/02—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part
- A61B2576/026—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part for the brain
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Evolutionary Computation (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Pathology (AREA)
- Medical Informatics (AREA)
- Heart & Thoracic Surgery (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Neurology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Psychiatry (AREA)
- Physiology (AREA)
- Radiology & Medical Imaging (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Mathematical Physics (AREA)
- Signal Processing (AREA)
- Child & Adolescent Psychology (AREA)
- High Energy & Nuclear Physics (AREA)
- Developmental Disabilities (AREA)
- Hospice & Palliative Care (AREA)
- Software Systems (AREA)
- Psychology (AREA)
- Neurosurgery (AREA)
- Computational Linguistics (AREA)
Abstract
An ICA-CNN classified fMRI space activation map smoothing and broadening method belongs to the field of biomedical signal processing. Firstly, three-dimensional space smoothing is applied to an fMRI space activation map obtained by ICA separation, a new sample set is generated by amplification, and then the new sample set is sent to CNN for classification, so that the improvement of classification performance of patients and healthy people is realized. The invention is adopted to classify 82 tested resting state fMRI data into patients and healthy people, three Gaussian filters with different FWHM are adopted to separate ICA to obtain a space activation map for space smoothing, three groups of sample sets are generated, and then the three groups of sample sets are sent to CNN for classification. Compared with the existing multi-model order data augmentation method, the method can improve the classification accuracy by 5.76%; by combining the two methods, the classification accuracy can be improved by 21.33%. Therefore, the invention can independently improve the network classification performance, is easy to combine with other augmentation methods, and obviously improves the classification accuracy.
Description
Technical Field
The invention belongs to the field of biomedical signal processing, and particularly relates to a smoothing and amplifying method for an ICA-CNN classified fMRI space activation map.
Background
Convolutional Neural Networks (CNNs) have the advantages of local sensing, weight sharing, and the like, are prominent in tasks such as recognition, detection, classification, and the like, and are also useful in smart medicine. The resting-state fMRI (rs-fMRI) data has the advantages of non-invasiveness, high spatial resolution, easy acquisition on the test of a patient and the like, and is commonly used for analyzing and diagnosing neurological disorders such as schizophrenia. Therefore, CNNs with rs-fMRI data as training data would represent a unique advantage in the task of classifying healthy persons from patients.
In view of the problem that fMRI of patients is difficult to acquire and the data size is not large, linkawa et al propose an ICA-CNN classification framework for complex fMRI data of patients and healthy persons (patent application No. 201910350137.2). In this framework, two-dimensional slices of spatial components of interest are first separated using Independent Component Analysis (ICA), and then fed into two-dimensional CNN to extract features and classify. Compared with a method of directly utilizing observation data or using three-dimensional CNN classification, the framework reduces the requirement on fMRI data volume, and achieves better classification performance under the condition of limited fMRI data volume.
In fact, in the case where the fMRI data amount is constant, the data expansion can serve to further increase the data amount and improve the ICA-CNN classification effect. In the field of image classification, widely used data augmentation methods include rotation, flipping, noise addition, scaling, texture color transformation, and the like. However, these methods are not suitable for rs-fMRI analysis where the spatial structure cannot be varied arbitrarily. The patent application number 201910350137.2 only provides a data augmentation method with multiple model orders, and has important practical value in exploring other methods suitable for fMRI data augmentation under the ICA-CNN classification framework.
Disclosure of Invention
The invention provides an fMRI space activation map augmentation method embedded into an ICA-CNN classification frame, which can effectively expand fMRI sample size and further improve the classification performance of an ICA-CNN network. The specific scheme is that three-dimensional space smoothing is applied to an fMRI space activation map obtained by ICA separation, a new sample set is generated by amplification, and then the new sample set is sent to CNN for classification, so that the improvement of classification performance of patients and healthy people is realized, as shown in figure 1. The method comprises the following concrete steps:
the first step is as follows: inputting the four-dimensional fMRI observation amplitude data of the tested kK is 1, …, K is the total number of trials, T is the number of scan time points, X, Y, Z is the three dimensions of brain space.
The second step is that: mixing X(k)K is 1, …, K, the space dimension is developed into one dimension, namely the size of the space dimension is equal to X multiplied by Y multiplied by Z, then the outside-brain voxels are removed, and only the inside-brain voxels are taken, and the space dimension is obtainedV is the number of endosomes, V<X×Y×Z。
The third step: principal Component Analysis (PCA) (principal Component analysis) pairK is 1, …, K, and reducing vitamin to obtainN is the model order and is less than or equal to T.
The fourth step: ICA isolation and extraction of the component of interest. Using Infmax algorithm pairsAnd (5) carrying out ICA separation when K is 1, … and K; selecting a spatial activation map of the interested component from N ICA separation components based on the principle of maximum correlation coefficient with a spatial component templateThe spatial composition template is selected from the literature (S.M.Smith, P.T.Fox et al, "ceramic of the woven's functional architecture along with actuation and" Proceedings of the National Academy of Sciences of the United States of America, vol.106, No.31Pp.13040-13045,2009).
The fifth step: to pairK is 1, …, K, and carries out brain extrasomatic zero filling to restore to a three-dimensional space activation map, namely
And a sixth step: the space is smooth. fMRI space activation map of test kK1, …, K, convolved with a three-dimensional gaussian filter to achieve spatial smoothing, as shown by:
in the formulaRepresenting a three-dimensional convolution calculation, g (x, y, z) is a three-dimensional zero-mean gaussian function defined as:
where σ is the standard deviation, the width of the filter is determined, and the larger σ is, the larger the smoothing degree is. Full width at Half maximum fwhm (fullbidth at Half maximum) is a common parameter for measuring the size of a gaussian kernel, having a linear relationship with sigma as follows,the FWHM is typically taken to be 2-3 times the voxel size. Smoothing all tested spatial activation maps with J gaussian filters with different FWHM sizes can expand each tested sample set from the original 1 group to J groups. For test k, the J sets of augmented sample sets were recorded asAs shown in fig. 2. The spatial activation map under each smoothing kernel is then further processed as follows:
wherein λ is1、λ2As weight parameters, they satisfy λ1+λ2=1,0<λ1<1,0<λ2<1, and usually takes λ1>λ2。
The seventh step: will be provided withK1, …, K, i 1, …, J, unfolded along the Z-axis as a stack of Z two-dimensional slices of size X Y, denoted asAnd selecting W slices with large activation information according to the effective activation positions of the interested components, wherein W is less than or equal to Z and is marked as L, … and L + W-1.
Eighth step: cutting W pieces of two-dimensional slicesVoxel values of K1, …, K, i 1, …, J, are normalized to [0,1]Is marked asThen attaching a corresponding label to each sliceIt is indicative of a healthy person,representing patients with schizophrenia, establishing an augmented CNN sample setThe sample set of each test sample is expanded from original W to W × J, i.e., the expansion is J times of original. All KThe sample set of each test was expanded to kxwxj.
The ninth step: disorganizing the sequence of the tested sample, and proportionally dividing the sample set into a training set, a verification set and a test set; such as: 60% of the sample set is taken as a training set, 20% is taken as a verification set, and 20% is taken as a test set.
The tenth step: and sending the training set into a CNN network for training. The CNN network structure used in the present invention is consistent with patent application No. 201910350137.2, and as shown in fig. 3, includes 2 convolutional layers, 2 max pooling layers, 1 full connection layer, and an output layer. Wherein, the size of the convolution layer core is 3 multiplied by 3, and the number of the convolution layer core is respectively 8 and 16; the size of the pooling layer core is 2 x 2, and the number of the pooling layer cores is 8 and 16 respectively; the full connection layer has 64 nodes, modified linear unit (ReLU) is used as an activation function, and the output layer uses Softmax to give the category to which each slice belongs. The batch processing size is 64, parameter updating is carried out by adopting an Adam algorithm, a two-class cross entropy function is used as a loss function of the network, and meanwhile, L is used2And (5) regularizing control weight and training a CNN model to have R rounds. At each round, calculating an accuracy ACC using the validation set; and saving the CNN model with the highest verification accuracy in the R round as the CNN model obtained by training. Let CNN verify asAccuracy ACC is defined as follows:
wherein TP represents true positive, i.e.TN indicates true negatives, i.e.FP represents a false positive, i.eFN indicates false negatives, i.e.
The eleventh step: and inputting the test set into the CNN model obtained by the eleventh training step to obtain the test accuracy ACC.
The twelfth step: and 5-fold cross validation is adopted, namely the ninth step to the tenth step are repeated for 5 times to obtain the test accuracy ACC of 5 times, and the average result is calculated and recorded as the final classification accuracy of the CNN.
The invention has the beneficial effects that: the invention focuses on the classification problem in the field of computer-aided diagnosis, and provides a data augmentation method aiming at fMRI data characteristics and embedded in an ICA-CNN classification framework, which is used for training a CNN network and improving classification results. For example, the invention is adopted to classify 82 rs-fMRI data to be tested into schizophrenic patients and healthy people, ICA separation is carried out on fMRI observation data through Infmax to obtain a spatial activation map, the spatial activation map of a Default Mode Network (DMN) is extracted, and the FWHM is respectively 6mm3、8mm3、10mm3The three Gaussian filters carry out three-dimensional spatial smoothing on all tested DMN spatial activation graphs to generate three groups of sample sets, and finally, the three groups of sample sets are sent to the CNN for training and testing. Compared with the multi-model order data amplification method provided by the patent 201910350137.2, the method can improve the classification accuracy by 15.76% under the condition of the same sample number. If the method is combined with a multi-model order augmentation method, three model orders of N-20, 60 and 100 are set for fMRI observation data, DMN space activation graphs are respectively extracted, space smoothing is applied to the DMN space activation graphs obtained by each model order to generate three groups of sample sets, and then all the three groups of sample sets are sent to a CNN network for classification. Compared with the method only using the multiple model order augmentation, the combination method can improve the classification accuracy by 21.33%. Therefore, the invention not only can independently improve the network classification performance, but also can be easily combined with other augmentation methods, thereby obviously improving the classification accuracy.
Drawings
FIG. 1 illustrates the location and usage steps of the present invention in the ICA-CNN classification framework.
FIG. 2 is a schematic representation of the present invention for spatially smooth augmentation of fMRI spatial activation maps.
Fig. 3 is a CNN network architecture used in the present invention.
Detailed Description
An embodiment of the present invention will be described in detail below with reference to the accompanying drawings. There were 82 rs-fMRI magnitude data tested, including 42 schizophrenic patients and 40 healthy people. Each test contained 146 scans, each scan sharing whole-brain data of voxels X × Y × Z53 × 63 × 46 153594, with an intra-brain voxel V62336 and voxel size of 3 × 3 × 3mm3。
The first step is as follows: inputting all 82 tested four-dimensional fMRI observation amplitude datak=1,…,82。
The second step is that: mixing X(k)K is 1, …,82, the spatial dimension is one-dimensional, i.e. the size of the spatial dimension is 53 × 63 × 46 or 153594, then the extra-brain voxels are removed and only the intra-brain voxels are taken, and the method is obtained
The third step: using PCA pairsk is 1, … and 82, dimension reduction is carried out, the model order N is 50, and the obtained product is obtained
The fourth step: ICA isolation and extraction of the component of interest. Using Infmax algorithm pairsAnd (4) carrying out ICA separation when k is 1, … and 82; based on the literature (S.M. Smith, P.T. Fox et al, "Correstondence of the blue's functional architecture during activation and rest," proceedings of the National Academy of Sciences of the United States of America, vol.106, No.31, pp.13040-13045,2009), selecting DMN space activation map from 50 ICA separation components
The fifth step: to pairk is 1, …,82, and carries out brain extrasomatic zero filling to restore to three-dimensional space activation map, namely
And a sixth step: the space is smooth. fMRI space activation map of test kk is 1, …,82, and is convolved with a three-dimensional gaussian filter to achieve spatial smoothing, as shown in equation (1). According to the definition of the formula (2), 3 three-dimensional Gaussian filters are generated, and the FWHM is 6mm3、8mm3、10mm3. Smoothing all tested DMN spatial activation maps by using the 3 Gaussian filters, and expanding each tested sample set from 1 group to 3 groups. For test k, 3 sets of the augmented sample set were recorded ask is 1, …,82, i is 1,2, 3. Then, the space activation map under each smooth kernel is further processed by the formula (3) to obtaink is 1, …,82, i is 1,2,3, and λ is taken1=0.8,λ2=0.2。
The seventh step: will be provided with k 1, …,82, i 1,2,3, unfolded along the z-axis into a stack of 46 two-dimensional slices of size 53 × 63, notedAccording to the effective activation position of the interested component, 25 slices with large activation information quantity are selected, specifically, L is 12, z is 12, … and 36.
Eighth step: 25 two-dimensional slicesVoxel values with k 1, …,82, i 1,2,3 are normalized to [0,1]Is marked asThen attaching a corresponding label to each sliceIt is indicative of a healthy person,representing patients with schizophrenia, establishing an augmented CNN sample setEach sample set tested was expanded from W25 to W × J25 × 3 to 75, i.e., 3 times the original sample set. The number of samples of all 82 samples tested was increased to 82 × W × J, 6150.
The ninth step: the sequence of the tested samples is disturbed, 60% of sample sets are taken as training sets, 20% are taken as verification sets, and 20% are taken as testing sets.
The tenth step: and (5) sending the training set into a CNN network for training, and training the CNN model R for 50 rounds. At each round, calculating the accuracy rate ACC by using the verification set according to the formula (3); and saving the CNN model with the highest verification accuracy in 50 rounds.
The eleventh step: and (4) inputting the test set into the CNN model obtained by the eleventh training step, and calculating according to the formula (3) to obtain the test accuracy ACC.
The twelfth step: and 5-fold cross validation is adopted, namely the tenth step to the twelfth step are repeated for 5 times to obtain the test accuracy ACC for 5 times, and the average result is calculated to obtain the final classification accuracy of 0.8522.
Claims (3)
1. A kind of ICA-CNN classified fMRI space activation map smooth augmentation method, apply three-dimensional space smoothing to fMRI space activation map that ICA separates, augment and produce the new sample set, then send into CNN to classify, realize the improvement to patient and classification performance of the healthy person; the method is characterized by comprising the following steps:
the first step is as follows: inputting the four-dimensional fMRI observation amplitude data of the tested kK is 1, …, K is the total number of trials, T is the number of scan time points, X, Y, Z is the three dimensions of brain space;
the second step is that: mixing X(k)K is 1, …, K, the space dimension is developed into one dimension, namely the size of the space dimension is equal to X multiplied by Y multiplied by Z, then the outside-brain voxels are removed, and only the inside-brain voxels are taken, and the space dimension is obtainedV is the number of endosomes, V<X×Y×Z;
The third step: PCA pairs Using principal component analysisK is 1, …, K, and reducing vitamin to obtainN is the model order, and N is less than or equal to T;
the fourth step: ICA separation and interesting component extraction, and Infmax algorithm pairK is 1, …, K, ICA separation is carried out, and a space activation map of the interested component is selected from N ICA separation components based on the principle that the correlation coefficient with a space component template is maximum
The fifth step: to pairK is 1, …, K, and carries out brain extrasomatic zero filling to restore to a three-dimensional space activation map, namely
And a sixth step: spatial smoothing, fMRI spatial activation map of the examined kK1, …, K, convolved with a three-dimensional gaussian filter to achieve spatial smoothing, as shown by:
in the formulaRepresenting a three-dimensional convolution calculation, g (x, y, z) is a three-dimensional zero-mean gaussian function defined as:
wherein, σ is standard deviation, determines the width of the filter, the larger σ is, the larger smoothing degree is, the full width at half maximum FWHM is a common parameter for measuring the size of Gaussian kernel, and has a linear relation with σ as follows,smoothing all tested space activation maps by using J Gaussian filters with different FWHM sizes, and increasing each tested sample set from 1 group to J groups, and recording J groups of increased sample sets as tested kThe spatial activation map under each smoothing kernel is then further processed as follows:
in the formula of1、λ2Is a weight parameter;
the seventh step: will be provided withK1, …, K, i 1, …, J, unfolded along the Z-axis as a stack of Z two-dimensional slices of size X Y, denoted asSelecting W slices with large activation information according to the effective activation positions of the interested components, wherein W is less than or equal to Z and is recorded as L, … and L + W-1;
eighth step: cutting W pieces of two-dimensional slicesVoxel values of K1, …, K, i 1, …, J, are normalized to [0,1]Is marked asThen attaching a corresponding label to each slice It is indicative of a healthy person,representing patients with schizophrenia, establishing an augmented CNN sample setEach sample testedThe sample set is expanded from the original W to W multiplied by J, namely the expansion is multiplied by J times, and all K tested sample sets are expanded to K multiplied by W multiplied by J;
the ninth step: disorganizing the sequence of the tested sample, and proportionally dividing the sample set into a training set, a verification set and a test set;
the tenth step: sending a training set into a CNN network for training, wherein the CNN network structure comprises 2 convolutional layers, 2 maximum pooling layers, 1 full-connection layer and an output layer, the size of each convolutional layer core is 3 multiplied by 3, the number of each convolutional layer core is 8 and 16, the size of each pooling layer core is 2 multiplied by 2, the number of each pooling layer core is 8 and 16, 64 nodes are arranged in the full-connection layer, a modified linear unit ReLU is used as an activation function, the output layer gives the category to which each slice belongs by using Softmax, the batch processing size is 64, parameter updating is performed by using an Adam algorithm, a two-classification cross entropy function is used as a loss function of the network, and meanwhile, L2Regularizing control weight, training CNN models in R rounds, calculating ACC by using a verification set in each round, storing the CNN model with the highest verification accuracy in the R rounds as the CNN model obtained by training, and making the CNN verification result beAccuracy ACC is defined as follows:
wherein TP represents true positive, i.e.TN indicates true negatives, i.e.FP represents a false positive, i.eFN indicates false negatives, i.e.
The eleventh step: inputting the test set into the CNN model obtained by the eleventh training step to obtain a test accuracy ACC;
the twelfth step: and 5-fold cross validation is adopted, namely the ninth step to the tenth step are repeated for 5 times to obtain the test accuracy ACC of 5 times, and the average result is calculated and recorded as the final classification accuracy of the CNN.
2. The ICA-CNN classified fMRI spatial activation map smooth augmentation method according to claim 1, wherein FWHM is 2-3 times of fMRI voxel size, λ1、λ2Satisfy lambda1+λ2=1,0<λ1<1,0<λ2<1, and λ1>λ2。
3. The ICA-CNN classified fMRI spatial activation map smoothing and augmenting method according to claim 1 or 2, wherein in the ninth step, 60% of the sample set is used as a training set, 20% is used as a verification set, and 20% is used as a test set.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911144802.9A CN110870770B (en) | 2019-11-21 | 2019-11-21 | ICA-CNN classified fMRI space activation map smoothing and broadening method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911144802.9A CN110870770B (en) | 2019-11-21 | 2019-11-21 | ICA-CNN classified fMRI space activation map smoothing and broadening method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110870770A true CN110870770A (en) | 2020-03-10 |
CN110870770B CN110870770B (en) | 2021-05-11 |
Family
ID=69718110
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911144802.9A Active CN110870770B (en) | 2019-11-21 | 2019-11-21 | ICA-CNN classified fMRI space activation map smoothing and broadening method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110870770B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113449820A (en) * | 2021-08-27 | 2021-09-28 | 深圳市商汤科技有限公司 | Image processing method, electronic device, and storage medium |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160292856A1 (en) * | 2015-04-06 | 2016-10-06 | IDx, LLC | Systems and methods for feature detection in retinal images |
US20180082443A1 (en) * | 2016-09-21 | 2018-03-22 | Realize, Inc. | Anomaly detection in volumetric images |
CN109222972A (en) * | 2018-09-11 | 2019-01-18 | 华南理工大学 | A kind of full brain data classification method of fMRI based on deep learning |
CN109376751A (en) * | 2018-08-21 | 2019-02-22 | 北京工业大学 | A kind of human functional network's classification method based on convolutional neural networks |
CN110110776A (en) * | 2019-04-28 | 2019-08-09 | 大连理工大学 | A kind of tranquillization state plural number fMRI data ICA-CNN taxonomy model of patient and Healthy People |
-
2019
- 2019-11-21 CN CN201911144802.9A patent/CN110870770B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160292856A1 (en) * | 2015-04-06 | 2016-10-06 | IDx, LLC | Systems and methods for feature detection in retinal images |
US20180082443A1 (en) * | 2016-09-21 | 2018-03-22 | Realize, Inc. | Anomaly detection in volumetric images |
CN109376751A (en) * | 2018-08-21 | 2019-02-22 | 北京工业大学 | A kind of human functional network's classification method based on convolutional neural networks |
CN109222972A (en) * | 2018-09-11 | 2019-01-18 | 华南理工大学 | A kind of full brain data classification method of fMRI based on deep learning |
CN110110776A (en) * | 2019-04-28 | 2019-08-09 | 大连理工大学 | A kind of tranquillization state plural number fMRI data ICA-CNN taxonomy model of patient and Healthy People |
Non-Patent Citations (1)
Title |
---|
MUHAMMAD NAVEED IQBAL QURESHI ET AL.: "3D-CNN based discrimination of schizophrenia using resting-state fMRI", 《ARTIFICIAL INTELLIGENCE IN MEDICINE》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113449820A (en) * | 2021-08-27 | 2021-09-28 | 深圳市商汤科技有限公司 | Image processing method, electronic device, and storage medium |
CN113449820B (en) * | 2021-08-27 | 2022-01-18 | 深圳市商汤科技有限公司 | Image processing method, electronic device, and storage medium |
Also Published As
Publication number | Publication date |
---|---|
CN110870770B (en) | 2021-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108182441B (en) | Parallel multichannel convolutional neural network, construction method and image feature extraction method | |
CN108491849B (en) | Hyperspectral image classification method based on three-dimensional dense connection convolution neural network | |
CN111242841B (en) | Image background style migration method based on semantic segmentation and deep learning | |
CN107665492B (en) | Colorectal panoramic digital pathological image tissue segmentation method based on depth network | |
Watson et al. | Patterns of response to visual scenes are linked to the low-level properties of the image | |
CN109410219A (en) | A kind of image partition method, device and computer readable storage medium based on pyramid fusion study | |
CN106462771A (en) | 3D image significance detection method | |
CN110110776B (en) | Method for constructing resting state complex fMRI data ICA-CNN classification framework of patient and healthy person | |
CN111090764B (en) | Image classification method and device based on multitask learning and graph convolution neural network | |
CN107292346B (en) | A kind of MR image hippocampus partitioning algorithm based on Local Subspace study | |
CN111009324A (en) | Mild cognitive impairment auxiliary diagnosis system and method based on brain network multi-feature analysis | |
CN104268561B (en) | High spectrum image solution mixing method based on structure priori low-rank representation | |
CN114782298B (en) | Infrared and visible light image fusion method with regional attention | |
CN110870770B (en) | ICA-CNN classified fMRI space activation map smoothing and broadening method | |
Dresp-Langley et al. | Using the quantization error from Self‐Organizing Map (SOM) output for fast detection of critical variations in image time series | |
CN111046918B (en) | ICA-CNN classified fMRI data space pre-smoothing and broadening method | |
CN115131503A (en) | Health monitoring method and system for iris three-dimensional recognition | |
CN107330948B (en) | fMRI data two-dimensional visualization method based on popular learning algorithm | |
CN107818567A (en) | Brain local morphological feature description method based on cortical top point cloud | |
Muralidharan et al. | Diffeomorphic shape trajectories for improved longitudinal segmentation and statistics | |
CN110916661B (en) | ICA-CNN classified fMRI intracerebral data time pre-filtering and amplifying method | |
Chen et al. | A fully convolutional network for MR fingerprinting | |
CN114187475A (en) | Method for explaining CNN classification result of multi-test complex fMRI data based on thermodynamic diagram | |
Li et al. | Collaborative Multi-View Convolutions With Gating For Accurate And Fast Volumetric Medical Image Segmentation | |
CN113095334B (en) | Contour detection method based on retina photopic adaptive mechanism |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |