CN103006220A - Method for detecting brain function communicated area based on signal sparse approximation - Google Patents

Method for detecting brain function communicated area based on signal sparse approximation Download PDF

Info

Publication number
CN103006220A
CN103006220A CN2012105826911A CN201210582691A CN103006220A CN 103006220 A CN103006220 A CN 103006220A CN 2012105826911 A CN2012105826911 A CN 2012105826911A CN 201210582691 A CN201210582691 A CN 201210582691A CN 103006220 A CN103006220 A CN 103006220A
Authority
CN
China
Prior art keywords
signal
functional
sparse
data
brain function
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
Application number
CN2012105826911A
Other languages
Chinese (zh)
Other versions
CN103006220B (en
Inventor
王倪传
曾卫明
王晗
姚胜南
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Maritime University
Original Assignee
Shanghai Maritime University
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 Shanghai Maritime University filed Critical Shanghai Maritime University
Priority to CN201210582691.1A priority Critical patent/CN103006220B/en
Publication of CN103006220A publication Critical patent/CN103006220A/en
Application granted granted Critical
Publication of CN103006220B publication Critical patent/CN103006220B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention discloses a method for detecting a brain function communicated area based on signal sparse approximation. The method comprises the following steps: 1, performing wavelet packet decomposition on each time point data in a functional magnetic resonance signal to obtain a wavelet tree on original time point data; 2, carrying out sparsity measurement on nodes in the wavelet tree, and selecting the node with strongest sparsity in the wavelet tree to form an effective sparse approximation data set on original functional magnetic resonance data; and 3, carrying out hybrid matrix optimization on the formed sparse approximation data set by adopting independent content analysis, and reconstructing a source signal in a functional area in combination with an original functional magnetic resonance mixed signal to finish accurate positioning and detection of the functional area. According to the method, sparse approximation of original mixed data is obtained by utilizing the more general supposed sparsity of the functional magnetic resonance signal, and the source signal is separated through the independent content analysis and the signal reconstruction, so that the purpose of accurately positioning the brain function communicated area is achieved.

Description

A kind of based on the sparse approximate brain function connected region detection method of signal
Technical field
The present invention relates to a kind of human brain function connected region detection method based on Functional magnetic resonance imaging, be specifically related to a kind of based on accurate location and the detection method of the sparse approximate blind source signal separation method of functional MRI mixed signal to the human brain function connected region.
Background technology
Functional magnetic resonance imaging is a kind of Novel magnetic resonance imaging technique that begins to rise the nineties in 20th century.This technology combines function, dissection and image three aspects: information, study the systematic study that combines with function for traditional mr techniques from single morphosis strong technical support is provided, simultaneously for providing favourable technical guarantee to diagnosis of research, brain section disease and the mental illness of human brain function connected region, neuro-cognitive etc.Utilizing Functional magnetic resonance imaging to carry out in the research of human brain function connected region, the method for processing effectively, accurately the higher-dimension MR data has played vital effect to the research of human brain function connected region.Yet current Data Management Analysis method although can finish to a certain extent the detection of functional area, is all deposited more deficiency and defective, and the accuracy of its Orientation of functions is also treated further raising.For example, method of fuzzy cluster analysis is limited by the restriction that number is estimated in iteration speed, Fuzzy Exponential and functional areas; Independent component analysis is limited by the separate hypothesis of stronger functional areas source signal fully, has limited the detection in Functional connectivity zone.
Therefore, existing human brain function connected region detection technique is also waited to further develop and improve.
Summary of the invention
The invention provides a kind of based on the sparse approximate brain function connected region detection method of signal, it is sparse approximate by the functional MRI mixed signal is carried out, obtain sparse approximate data collection, adopt again Independent Component Analysis to carry out that brain function signal is blind to be cut apart, realize more accurately detection and localization brain function connected region.
For achieving the above object, the invention provides and a kind ofly be characterized in based on the sparse approximate brain function connected region detection method of signal, the method includes the steps of:
Step 1, each the time point data in the functional magnetic resonance signal are carried out WAVELET PACKET DECOMPOSITION, obtain the wavelet tree about former time point data;
Step 2, the node in the wavelet tree is carried out sparse property tolerance, and select the strongest wavelet tree node of sparse property, so that form effective sparse property approximate data collection about former functional MRI data;
Step 3, the sparse property approximate data centralized procurement that forms is carried out the hybrid matrix optimizing with independent component analysis, and carry out the reconstruct of functional areas source signal in conjunction with former functional MRI mixed signal, to finish accurate location and the detection of functional area.
In the above-mentioned steps 1, the time point data are carried out respectively three layers of one dimension WAVELET PACKET DECOMPOSITION.
The present invention is a kind of to be compared based on the sparse approximate brain function connected region detection method of signal and the Data Management Analysis method of prior art, its advantage is, the present invention is by the sparse property of the more general hypothesis of at first utilizing functional magnetic resonance signal, obtain the sparse approximate of former blended data, again by independent component analysis and signal reconstruction to realize separating of source signal, thereby reach the purpose of accurate location brain function connected region, be conducive to the research to aspects such as brain science, neuroscience and brain section diseases.
Description of drawings
Fig. 1 is a kind of method flow diagram based on the sparse approximate brain function connected region detection method of signal of the present invention;
Fig. 2 is a kind of flow chart based on the sparse approximate brain function connected region detection method of signal of the present invention.
The specific embodiment
Below in conjunction with accompanying drawing, further specify specific embodiments of the invention.
Such as Fig. 1 and in conjunction with shown in Figure 2, the present invention discloses a kind of based on the sparse approximate brain function connected region detection method of signal, and the method includes the steps of:
Step 1, WAVELET PACKET DECOMPOSITION: to each time point data of functional MRI data, carry out respectively 3 layers of one dimension WAVELET PACKET DECOMPOSITION, the wherein selection of wavelet basis has been chosen the Daubechies(that possesses simultaneously anti symmetry, orthogonality and biorthogonality in the little wave system and slightly has been written as: db) the db2 wavelet basis in the family among the present invention; After WAVELET PACKET DECOMPOSITION, can obtain the corresponding wavelet tree of each time point data (as shown in Figure 2), each node has different sparse attributes in the wavelet tree., and formed a kind of effectively sparse approximate about original signal based on the wavelet tree node of gained
Step 2, the sparse approximate data collection operation of formation: this operation has mainly comprised the work of two aspects.
Step 2.1, one are the sparse property tolerance device of the wavelet tree node of a distance metric norm based on normed linear space of structure (lp norm), sparse property to each wavelet tree node is measured, obtain one about the sparse property amount vector Q of each node, it satisfies, and wherein its expression is about the wavelet packet tree of former each time point of functional MRI data the lThe common wavelet coefficient vector that forms of individual wavelet tree node, and n represents the length of vector.
Step 2.2, it is two according to sparse property size, selects the wavelet tree node of sparse property maximum to form a kind of effectively sparse approximate collection about former mixed signal, thereby has improved reconstruct and the stationkeeping ability of brain function signal.
Step 3, independent component analysis and reconstruct: this flow process mainly is that the sparse approximate collection about former mixed signal to the formation in the step 2 carries out independent component analysis, so that estimate about each source signal mixed coefficint A in the former functional MRI data and the mixed sparse W of solution, unite again former mixed function MR data signal, carry out the source signal reconstruct of (or being called functional network), finally obtain needed Functional connectivity network, i.e. independent element image and corresponding time course.
Although content of the present invention has been done detailed introduction by above preferred embodiment, will be appreciated that above-mentioned description should not be considered to limitation of the present invention.After those skilled in the art have read foregoing, for multiple modification of the present invention with to substitute all will be apparent.Therefore, protection scope of the present invention should be limited to the appended claims.

Claims (2)

1. one kind based on the sparse approximate brain function connected region detection method of signal, it is characterized in that, the method includes the steps of:
Step 1, each the time point data in the functional magnetic resonance signal are carried out WAVELET PACKET DECOMPOSITION, obtain the wavelet tree about former time point data;
Step 2, the node in the wavelet tree is carried out sparse property tolerance, and select the strongest wavelet tree node of sparse property, so that form effective sparse property approximate data collection about former functional MRI data;
Step 3, the sparse property approximate data centralized procurement that forms is carried out the hybrid matrix optimizing with independent component analysis, and carry out the reconstruct of functional areas source signal in conjunction with former functional MRI mixed signal, to finish accurate location and the detection of functional area.
2. as claimed in claim 1ly it is characterized in that based on the sparse approximate brain function connected region detection method of signal, in the described step 1, the time point data are carried out respectively three layers of one dimension WAVELET PACKET DECOMPOSITION.
CN201210582691.1A 2012-12-28 2012-12-28 Method for detecting brain function communicated area based on signal sparse approximation Expired - Fee Related CN103006220B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210582691.1A CN103006220B (en) 2012-12-28 2012-12-28 Method for detecting brain function communicated area based on signal sparse approximation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210582691.1A CN103006220B (en) 2012-12-28 2012-12-28 Method for detecting brain function communicated area based on signal sparse approximation

Publications (2)

Publication Number Publication Date
CN103006220A true CN103006220A (en) 2013-04-03
CN103006220B CN103006220B (en) 2015-02-11

Family

ID=47955716

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210582691.1A Expired - Fee Related CN103006220B (en) 2012-12-28 2012-12-28 Method for detecting brain function communicated area based on signal sparse approximation

Country Status (1)

Country Link
CN (1) CN103006220B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103211597A (en) * 2013-04-27 2013-07-24 上海海事大学 Resting brain function connected region detecting method based on affine clustering
CN103961103A (en) * 2014-05-07 2014-08-06 大连理工大学 Method for performing phase correction on ICA estimation components of plural fMRI data
CN104850863A (en) * 2015-05-29 2015-08-19 上海海事大学 Human brain functional activity state classification method
CN108197661A (en) * 2018-01-19 2018-06-22 上海海事大学 Cognitive activities state classification system and method based on voxel level brain function information

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090264785A1 (en) * 2008-04-18 2009-10-22 Brainscope Company, Inc. Method and Apparatus For Assessing Brain Function Using Diffusion Geometric Analysis
WO2010045724A1 (en) * 2008-10-20 2010-04-29 Neurochip Corporation Method and rhythm extractor for detecting and isolating rhythmic signal features from an input signal using the wavelet packet transform
US20100191139A1 (en) * 2009-01-28 2010-07-29 Brainscope, Inc. Method and Device for Probabilistic Objective Assessment of Brain Function
US20110130675A1 (en) * 2009-12-01 2011-06-02 Neurowave Systems Inc. Multi-channel brain or cortical activity monitoring and method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090264785A1 (en) * 2008-04-18 2009-10-22 Brainscope Company, Inc. Method and Apparatus For Assessing Brain Function Using Diffusion Geometric Analysis
WO2010045724A1 (en) * 2008-10-20 2010-04-29 Neurochip Corporation Method and rhythm extractor for detecting and isolating rhythmic signal features from an input signal using the wavelet packet transform
US20100191139A1 (en) * 2009-01-28 2010-07-29 Brainscope, Inc. Method and Device for Probabilistic Objective Assessment of Brain Function
US20110130675A1 (en) * 2009-12-01 2011-06-02 Neurowave Systems Inc. Multi-channel brain or cortical activity monitoring and method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
曾卫明等: "基于小波分析的功能磁共振图像的去噪预处理", 《应用科学学报》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103211597A (en) * 2013-04-27 2013-07-24 上海海事大学 Resting brain function connected region detecting method based on affine clustering
CN103211597B (en) * 2013-04-27 2014-12-17 上海海事大学 Resting brain function connected region detecting method based on affine clustering
CN103961103A (en) * 2014-05-07 2014-08-06 大连理工大学 Method for performing phase correction on ICA estimation components of plural fMRI data
CN103961103B (en) * 2014-05-07 2015-12-30 大连理工大学 A kind of ICA to plural fMRI data is estimated the method being divided into line phase and correcting
CN104850863A (en) * 2015-05-29 2015-08-19 上海海事大学 Human brain functional activity state classification method
CN104850863B (en) * 2015-05-29 2017-11-17 上海海事大学 A kind of human brain function activity state classification method
CN108197661A (en) * 2018-01-19 2018-06-22 上海海事大学 Cognitive activities state classification system and method based on voxel level brain function information
CN108197661B (en) * 2018-01-19 2022-03-01 上海海事大学 Cognitive activity state classification system and method based on voxel horizontal brain function information

Also Published As

Publication number Publication date
CN103006220B (en) 2015-02-11

Similar Documents

Publication Publication Date Title
Wang et al. Single slice based detection for Alzheimer’s disease via wavelet entropy and multilayer perceptron trained by biogeography-based optimization
Akiki et al. Determining the hierarchical architecture of the human brain using subject-level clustering of functional networks
Lu et al. A pathological brain detection system based on radial basis function neural network
Catani et al. Comment on “The geometric structure of the brain fiber pathways”
Chen et al. OctopusNet: a deep learning segmentation network for multi-modal medical images
CN103006220B (en) Method for detecting brain function communicated area based on signal sparse approximation
CN103093455A (en) Diffusion tensor imaging white matter fiber clustering method
CN116503680B (en) Brain image structured analysis and brain disease classification system based on brain atlas
Ye et al. The intrinsic geometry of the human brain connectome
Pardakhti et al. Brain age estimation based on 3D MRI images using 3D convolutional neural network
Goceri Comparison of the impacts of dermoscopy image augmentation methods on skin cancer classification and a new augmentation method with wavelet packets
Yeswanth et al. Super resolution-based leaf disease detection in potato plant using broad deep residual network (BDRN)
Kushwaha et al. 3D medical image fusion using dual tree complex wavelet transform
Liang et al. Development of a framework for large scale three-dimensional pathology and biomarker imaging and spatial analytics
CN115984257A (en) Multi-modal medical image fusion method based on multi-scale transform
Niu et al. Fusion of multiple spatial networks derived from complex-valued fMRI data via CNN classification
Li et al. Decoding motor imagery with a simplified distributed dipoles model at source level
Xi et al. Multimodal medical volumetric image fusion based on multi-feature in 3-D shearlet transform
Chen et al. A new approach for functional connectivity via alignment of blood oxygen level-dependent signals
CN103211597B (en) Resting brain function connected region detecting method based on affine clustering
Dai et al. Network-based classification using cortical thickness of AD patients
CN103838942A (en) Cerebral function imaging diagnostic method based on knowledge base
Zawadzka et al. Graph representation integrating signals for emotion recognition and analysis
Zaabi et al. Deep and statistical-based methods for alzheimer's disease detection: A survey
Bae et al. Structural and functional diversity of a dense sample of retinal ganglion cells

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150211

Termination date: 20181228

CF01 Termination of patent right due to non-payment of annual fee