CN1632830A - Automatic dividing method for cerebral ischemia focus area - Google Patents
Automatic dividing method for cerebral ischemia focus area Download PDFInfo
- Publication number
- CN1632830A CN1632830A CN 200310121174 CN200310121174A CN1632830A CN 1632830 A CN1632830 A CN 1632830A CN 200310121174 CN200310121174 CN 200310121174 CN 200310121174 A CN200310121174 A CN 200310121174A CN 1632830 A CN1632830 A CN 1632830A
- Authority
- CN
- China
- Prior art keywords
- image
- voxel
- local volume
- dti
- disperse
- 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
- 238000000034 method Methods 0.000 title claims abstract description 65
- 206010008118 cerebral infarction Diseases 0.000 title claims description 28
- 201000006474 Brain Ischemia Diseases 0.000 title claims description 27
- 206010008120 Cerebral ischaemia Diseases 0.000 title claims description 27
- 230000011218 segmentation Effects 0.000 claims description 31
- 239000006185 dispersion Substances 0.000 claims description 14
- 230000000694 effects Effects 0.000 claims description 14
- 238000005259 measurement Methods 0.000 claims description 11
- 238000009826 distribution Methods 0.000 claims description 8
- 238000013507 mapping Methods 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 3
- 230000008520 organization Effects 0.000 claims description 2
- 238000003045 statistical classification method Methods 0.000 claims description 2
- 210000004556 brain Anatomy 0.000 abstract description 7
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 abstract description 6
- 201000010099 disease Diseases 0.000 abstract description 5
- 238000004458 analytical method Methods 0.000 abstract description 2
- 239000008280 blood Substances 0.000 abstract 1
- 210000004369 blood Anatomy 0.000 abstract 1
- 210000001519 tissue Anatomy 0.000 description 19
- 238000002597 diffusion-weighted imaging Methods 0.000 description 9
- 230000003044 adaptive effect Effects 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 7
- 210000004885 white matter Anatomy 0.000 description 6
- 238000003384 imaging method Methods 0.000 description 5
- 230000001154 acute effect Effects 0.000 description 4
- 210000003792 cranial nerve Anatomy 0.000 description 4
- 238000003709 image segmentation Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 208000006011 Stroke Diseases 0.000 description 2
- 210000001175 cerebrospinal fluid Anatomy 0.000 description 2
- 238000007596 consolidation process Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 239000000835 fiber Substances 0.000 description 2
- 239000004615 ingredient Substances 0.000 description 2
- 238000002156 mixing Methods 0.000 description 2
- 201000006417 multiple sclerosis Diseases 0.000 description 2
- 210000005036 nerve Anatomy 0.000 description 2
- 238000012113 quantitative test Methods 0.000 description 2
- 206010051290 Central nervous system lesion Diseases 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000007177 brain activity Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 208000026106 cerebrovascular disease Diseases 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 210000004884 grey matter Anatomy 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000005311 nuclear magnetism Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000007170 pathology Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
This invention relates to image process technique, and especially to an automatic division of brain blood shortage area based on multiple-size statistic sorting and local container sorting method, which comprises the following steps: to valuate the DTI image spreading value and direction isomerism; second to computer the size space; third to sort the multiple size; fourth to sort the local container. This invention is of high application value in medical assistant dialogue system, medical image three-dimensional recreation system and clinic disease qualitatively diagnose analysis.
Description
Technical field
The present invention relates to image processing techniques, particularly a kind of based on multiple dimensioned statistical classification and local volume sorting technique automatic division method to the cerebral ischemia focal zone.
Background technology
So-called image segmentation is meant that the zones of different that will have special connotation in the image makes a distinction, and these zones are Uncrossed mutually, and the consistance of specific region is all satisfied in each zone.From the process object angle, cutting apart is the location of definite target of being concerned about in image array.Obviously, have only " interested target object " extracted from the scene of complexity, just might further carry out quantitative test or identification, and then image is understood each sub regions.Image segmentation comprises Threshold Segmentation, rim detection, statistical classification or the like method.Image segmentation can with feature comprise gradation of image, color, texture, partial statistics characteristic or spectrum signature etc., utilize the difference of these features can differentiate between images in the different target object.Since we can only utilize some Partial Feature cut zone in the image information, so the whole bag of tricks must have limitation and specific aim, can only select suitable dividing method at the demand of various practical application area.
At present, the most important practical application area of image partition method is cutting apart of medical image.Medical image comprises CT, MR, and the image that other medical image equipment of reaching is obtained, the research majority of medical image segmentation is an example at the MR image or with the MR image at present.The research of medical image cutting method has two outstanding features, one is the domain knowledge that generally will use in the medical science, general shape as ventricle, the content of encephalic white matter and grey matter and relative position relation or the like, another is the mode that often adopts three-dimensional segmentation, this is because only have 2-D data in the general image, it is the two-dimensional projection that three-dimensional scene obtains by gamma camera or other image device, then directly provided the three-dimensional data with two dimension slicing form tissue in the medical image, this just provides possibility for three-dimensional segmentation.At medical image, concrete dividing method has many kinds, as based on Threshold Segmentation, region growing and division merging, sorter and the cluster in zone and based on method of random field etc.; Also have parallel differentiating operator based on the edge, surface fitting method, boundary curve fitting process etc. in addition.
In the medical imaging mode, the disperse NMR imaging is the important means of clinical diagnosis, especially for the headstroke disease.Headstroke is a kind of very serious disease, may cause people's lifelong deformity, even dead.Diffusion-weighted NMR imaging (diffusion weighted magnetic resonanceimaging is called for short DWI) technology is very important and effectively detects the clinical means of cerebral ischemia disease.Especially DWI can detect focus in the acute stage of falling ill, and this is that conventional NMR imaging technology is incomparable.Accurate position and the size that detects the cerebral ischemia focus helps classification of diseases, estimates disease condition, and instructs treatment.Yet, owing to be subjected to the T2 weighting, and rotation density, the influence of signals such as T1 weighting, and the influence of white matter of brain machine direction, the contrast of DWI image is very low; And only utilize the DWI image that sufficient disperse information can not be provided.Therefore, clinically at present, dispersion tensor NMR imaging (diffusion tensor magnetic resonance imaging is called for short DTI) technology more and more is used for quantitatively estimating, is analyzed the disperse feature of focal zone.
A large amount of scholars utilizes DWI or DTI technology, and research cerebral ischemia patient is along with the transformation of the state of an illness, the variation of focal zone volume; And utilize DTI technology, disperse isotropy, the anisotropy of quantitative test focal zone hydrone.In the former research, adopt manual methods to cut apart the cerebral ischemia focal zone, be very consuming time but manually cut apart, and segmentation result depends on operator's subjective determination more.Because the influence of factors such as noise, local volume effect, intensity non-uniformity and disperse anisotropy, up to the present, the cerebral ischemia zone of the automatic or semi-automatic DWI of cutting apart, DTI image remains the problem of difficulty.The generation of local volume effect is because the limited spatial resolution of nuclear magnetic scanning coil causes; The generation of intensity non-uniformity is because the unevenness (radio frequency inhomogeneities) of radio-frequency pulse causes; Intensity is overlapping be since in DWI or DTI image the similarity of cerebral ischemia focal zone and cranial nerve fibre strength cause.
Document about the automatic or semi-automatic cerebral ischemia focal zone of cutting apart DWI, DTI image is very limited.Martel etc. have proposed a kind of automanual dividing method, and the adaptive threshold that promptly absorbs space constraint is cut apart; And adopt iterated conditional pattern (iterative conditional mode is called for short ICM) method to seek locally optimal solution.But, because the overlapping influence of intensity that the disperse anisotropy is caused, differentiation focal zone and nerve tract that they can not be satisfied with.
For cutting apart of the focal zone of conventional MR image, adopt dividing method usually based on collection of illustrative plates (atlas-based).Utilize non-linear method for registering, dissect the different normal anatomical structures of identification that template can be successful.Yet focal zone can not utilize the normal template of dissecting to generate, and therefore just can not utilize collection of illustrative plates directly to cut apart the acquisition focal zone.Employing statistical classification methods such as Leemput are cut apart the different tissues structure of nuclear-magnetism (MR) image automatically, utilize initial segmentation and the geometric constraints of normal person's collection of illustrative plates as normal structure, and the brain lesion district is made as tissue not in the know (outliers).Cutting apart of the multiple sclerosis that is applied to the MR image of this method success (multiple sclerosis) focal zone.Based on the method for Leemput, Moon etc. have changed the space collection of illustrative plates, and lesion tissue as priori, is roughly set the position of focal zone.But because the low contrast and the anisotropy of disperse image, these methods all are not suitable for the cerebral infarction Region Segmentation of DTI image.
Some solve local volume effect method and are proposed successively.The mixing situation of the different tissues of voxel is represented in utilizations such as Laidlaw based on the grey level histogram in voxel zone, and utilizes the Bayesian probability approach to mate grey level histogram, mixes situation with the most probable different tissues of adjudicating single voxel inside.But the intensity of a voxel is what are organized that class mixes and unclear by actually, and this algorithm has been ignored the influence of the unevenness of intensity.Shattuck etc. and Noe etc. as a kind of new tissue class, and then estimate the mixing situation of pure tissue possible in each voxel with local volume voxel.But this method may cause the excessively level and smooth of classification results.
Rajapakse etc. utilize a statistical model to represent the distribution of mainly organizing class, and the ionization meter model has taken into full account the influence of noise and intensity non-uniformity.The robustness of this algorithm and degree of accuracy have fully been proved for the segmentation result of many clinical conventional brain MR images.But this algorithm still has its binding character: do not consider local volume effect; Segmentation result depends on initial segmentation; This algorithm is not suitable for cutting apart of focal zone, even for conventional brain MR image.
Based on the characteristics and the existing image processing method of DTI image itself, we propose a kind of brand-new adaptive method, cut apart the cerebral ischemia focal zone in the headstroke patient DTI image automatically.This method has taken into full account the influence of noise, disperse anisotropic, local volume effect and intensity non-uniformity.
Summary of the invention
The purpose of this invention is to provide a kind of method of cutting apart DTI image cerebral ischemia focal zone automatically, this method has taken into full account the influence of noise, disperse anisotropic, local volume effect and intensity non-uniformity, proposed the automatic division method based on the classification of multiple dimensioned statistical classification and local volume, assist clinicians is qualitative, quantitative Diagnosis and guidance treatment.
Core concept of the present invention is that we propose a kind of brand-new adaptive method, cuts apart the cerebral ischemia focal zone in the headstroke patient DTI image automatically.This method comprises following step: the image pre-service, calculate the dispersion tensor field, the anisotropic measurement of disperse, self-adapting multi-dimension statistical classification (multi-scale statistical classification, be called for short MSSC), and local volume voxel is classified (partial volume voxel reclassification is called for short PVVR) again.Adaptive M SSC model is considered the DTI image space information, information such as intensity gradient, disperse anisotropic and tissue characteristics; The PVVR model utilizes local parameter information to improve the accuracy that local volume is cut apart.
Based on above-mentioned purpose and thought, cut apart DTI image cerebral ischemia focal zone algorithm automatically based on multiple dimensioned statistical classification and local volume sorting technique and comprise:
(1) image pre-service is carried out filtering to original image;
(2) dispersion tensor field (tensor field) calculates, and obtains the pairing tensor field of each voxel of three dimensions;
(3) the anisotropic measurement of disperse quantizes the anisotropy of each voxel of three dimensions;
(4) calculating of metric space comprises original DTI graphical rule space, and disperse anisotropy mapping graph metric space;
(5) self-adapting multi-dimension statistical classification (multi-scale statistical classification, be called for short MSSC), in the focus zone of seeking based on optimum, overcome the influence that is brought owing to noise, disperse anisotropic and intensity non-uniformity;
(6) local volume voxel is classified (partial volume voxel reclassification is called for short PVVR) again, on the basis of self-adapting multi-dimension statistical classification, further optimizes segmentation result, overcomes the interference that local volume effect brings.
The present invention utilizes self-adapting multi-dimension statistical classification and local volume voxel to classify again, effectively overcomes the influence of noise, disperse anisotropic, local volume effect and intensity non-uniformity, has realized cutting apart automatically DTI image cerebral ischemia focal zone.In fields such as medical image assistant diagnosis system, medical image three-dimensional reconstructing system and clinical pathology qualitative, quantitative diagnostic analysis, there is important use to be worth.
Description of drawings
Fig. 1. carry out the method structured flowchart that DTI image cerebral ischemia focal zone is cut apart automatically based on multiple dimensioned statistical classification and local volume sorting technique;
Fig. 2. a tangent plane of an acute cerebral stroke patient's dispersion tensor volume data; Each image subsection is represented the scalar mapping graph of the ingredient of this aspect dispersion tensor D among the figure;
Fig. 3. to the disperse anisotropy measurement of Fig. 1 tensor tangent plane; Wherein: (a) trace figure; (b) FA figure;
Fig. 4. the metric space diagram;
Border condition after Fig. 5 .MSSC classification between possible different tissues class, wherein: (a) show by two zones of organizing class to form; (b) show by three borders of organizing class to form;
Fig. 6. the segmentation result to the different dividing methods of DTI image compares, wherein: (a) original cerebral ischemia patient's DTI image; (b) the adaptive M AP under the condition of initial segmentation is preferably cut apart, and wherein the arrow points mistake is divided into the cranial nerve bundle of focal zone; (c) MSSC is cut apart; (d) MSSC-PVVR is cut apart;
Fig. 7. the segmentation result to the different dividing methods of DTI image compares, wherein: (a) original cerebral ischemia patient's DTI image; (b) the adaptive M AP under the condition of initial segmentation is preferably cut apart, and wherein the arrow points mistake is divided into the cranial nerve bundle of focal zone; (c) MSSC is cut apart; (d) MSSC-PVVR is cut apart;
Fig. 8. the segmentation result to the different dividing methods of DTI image compares (a) original cerebral ischemia patient's DTI image; (b) the adaptive M AP under the condition of initial segmentation is preferably cut apart, and wherein the arrow points mistake is divided into the cranial nerve bundle of focal zone; (c) MSSC is cut apart; (d) MSSC-PVVR is cut apart.
Embodiment
Describe automatic division method of the present invention in detail below in conjunction with accompanying drawing.As a kind of concrete implementation, structured flowchart is seen Fig. 1, and this dividing method comprises following step: the image pre-service, and calculate the dispersion tensor field, the anisotropic measurement of disperse, the self-adapting multi-dimension statistical classification, and local volume voxel is classified again.
Mainly comprise four steps: (1) estimates the dispersion tensor and the disperse anisotropy of DTI image; (2) calculate metric space; (3) multiple dimensioned statistical classification; (4) local volume voxel is classified again.Below it is made introductions all round.
Step 1: dispersion tensor and the disperse anisotropy of estimating the DTI image
DTI is used for measuring the disperse characteristic of hydrone in biological tissue inside.Disperse is the process of a three-dimensional.But because the design feature of biological tissue, the motion of hydrone is not what equate in all directions of three-dimensional.Usually utilize the kinetic characteristic of the complete description hydrone of dispersion tensor along different coordinate axis, and the interdependence between them.
D represents dispersion tensor.From original DTI image calculation tensor D more complicated, detailed method Bihan is described in detail in list of references.Fig. 2 represents an aspect of an acute cerebral stroke patient's dispersion tensor volume data.Each subimage is represented the scalar mapping graph of the ingredient of this aspect dispersion tensor D among Fig. 2.
Utilize the eigenvalue of tensor D
1, λ
2, λ
3, λ wherein
1>λ
2>λ
3, can obtain different anisotropy measurements.We adopt distribution anisotropy's (fractional anisotropy is called for short FA) herein, quantitatively estimate the disperse anisotropy.
Tr(D)=(λ
1+λ
2+λ
3)/3
Wherein Tr (D) is the mark of D, represents the average dispersity of each different directions.Fig. 3 has showed Fig. 2 tensor tangent plane
Mark (Trace) figure and distribution anisotropy (FA) figure.The anisotropy that can observe normal brain activity white matter zone from Fig. 3 (b) is apparently higher than ectocinerea and cerebrospinal fluid zone.
Step 2: calculate metric space
Metric space can produce based on many different criterions.The key character of the easy blurred picture of linear-scale space technology, for example border of different tissues structure.Non-linear metric space has overcome this shortcoming, makes regional interior level and smooth dynamics greater than interregional level and smooth.Perona and Malik have proposed a partial differential equation of calculating non-linear metric space.Central idea is to have introduced rim detection in filtering, allows different scale inter-step mutual.
c(i,t)=g(|y(i,t)|)
In our dividing method, (i t) represents the DTI image in yardstick level t, the image intensity value of position i to y; (i t) is diffusion parameter to c, depends on the locus and changes, and is the function of the mould that changes of image intensity; Operator is dispersed in the div representative; and Δ are represented gradient operator and Laplce (Laplacian) operator respectively.At this paper, we adopt following spread function g () to produce metric space.
The manual method of being set or using Canny to propose of constant Ks or quilt produces.
Fig. 4 has shown the metric space synoptic diagram.The different different image spatial resolutions of yardstick level t representative.Yardstick level t=0 represents original image, the highest resolution in the metric space.Along with the yardstick level increases, image is more and more fuzzyyer, and it is fewer and feweri that institute comprises image information.The image of high-resolution image and low resolution has close corresponding relation, and this is convenient to our the image extraction one-piece construction information from low resolution, and obtains detailed information from high-resolution image.
Step 3: multiple dimensioned statistical classification (MSSC)
Piece image is to be made of mutual set of voxels adjacent one another are.I represents the coordinate position of voxel in image, y=(y
i, i ∈ I) and the intensity level of representative image, y
iRepresentative image is at the intensity level of voxel location i.The number of organizing class of setting image is K, each organize class with one gather Λ=1, a numeric representation among the 2...k}; To cutting apart in fact exactly of image each voxel being classified as the different classes of organizing.x
i=k representative belongs at the voxel of position i organizes class k, k ∈ Λ.X=(x
i, i ∈ I) and a kind of classification results of representative image y.
The process of cutting apart is exactly to seek suitable x, the correct class of organizing of coming representative image y to belong at each voxel location.We adopt maximum a posteriori estimation (maximum a posteriori is called for short MAP) that original image is cut apart.If x=x
*Represent optimum segmentation
Wherein Ω represents all possible cutting apart, the probability that the acquisition x of p (x|y) representative under known image y condition cut apart.Because prior probability p (y) is independent of and cuts apart x, according to the Bayesian theory,
p(x|y)∝p(x,y)=p(y|x)p(x)
The noise of supposing image meets the white Gaussian noise distribution, if x
i=k, so
μ wherein
K,, i, n
K,, IAnd σ
K,, iRepresent the strength mean value of organizing class k at position i respectively, noise, and the standard variance of noise.θ={ θ
i, i ∈ I} represents the parameter set of measurement model, θ
i={ θ
K, i=(μ
K,, i, σ
K,, i), k ∈ Λ }.The average and the standard variance of class each voxel in image organized in representative, considers the variation of biological tissue and the unevenness of image space intensity.
Suppose R
kRepresent all to belong to the voxel location of organizing class k, conditional probability p (y|x) can write so
We utilize markov random file (Markov Random Field is called for short MRF) to define prior model, and then Probability p (x) is obeyed the Gibbs distribution,
Wherein β is the consolidation constant, and C represents the set of tissue points, and adjacent one another are between the voxel in the set.
We obtain
p(x|y)∝exp{-U(x)}
U (x) is energy function, and is as follows
Preceding two mutual restriction relations of representing between raw image data and the segmentation result of U (x); Last represents the smoothness constraint of prior model to segmentation result.The maximum a posteriori that searching is cut apart estimates that (MAP) estimation problem is equivalent to the minimization problem of energy function U (x).
In step 2, we obtain the sequence image under different resolution.When we increased the yardstick rank, the detailed information of image descended.The image sequence of obfuscation y (t)=(y (i, t), i ∈ I, t ∈ N), N={1,2...n} represents.(t represents the yardstick rank, the yardstick level that the n representative is the highest).Y (t) cut apart accordingly with x (t) represent, x (t)=(x (i, t), i ∈ I).
If obtained to cut apart x (t+1) at yardstick level t+1, estimated parameter θ (t+1) accordingly so just.Use known parameter θ (t+1), we just can obtain new the cutting apart at yardstick level t image y (t).Cut apart at the next one that yardstick level t+1 segmentation result x (t+1) relevant parameters θ (t+1) can be used to obtain at yardstick level t.Finished N the iteration of cutting apart with parameter estimation, when arrival had the original image (t=0) of highest resolution, we had just obtained final optimum segmentation result.Therefore, cut apart and can describe with following two processes:
Known each organize the model parameter of class, estimate most probable cutting apart; Known cutting apart, estimate the optimization model parameter.The selection of measurement model parameter is the maximization principle according to view data.Self-adapting multi-dimension statistical classification (MSSC) process is through iteration repeatedly, the minimum value of coming estimated energy function U (t).U
i(t) representative is at the local energy function of yardstick level t position i, as shown in the formula:
μ wherein
K,, i(t) and σ
K,, i(t) be the μ of DTI image respectively at yardstick level t
K,, iAnd σ
K,, iHere V
c(x
i(t)) be defined as and satisfy x
j(t)=x
i(t), i, j ∈ C, the voxel number.
In the DTI image, the cerebral ischemia focal zone shows high signal, and because anisotropic influence, white matter of brain nerve fibre place also presents high signal.Because intensity level is overlapping, accurately be partitioned into very difficulty of cerebral ischemia focal zone.We solve this problem by add the anisotropic control of disperse in model.
α
i(t)=|y
i(t)-r
i(t)|
r
i(t)=α FA
i(t) r
i(t) the original DTI image of representative is worth FA the distribution anisotropy of the position i of yardstick level t
i(t).Factor a guarantees that FA and y have consistent Strength Space.α
i(t) y is described
i(t) and r
i(t) otherness.We revise energy function U
i(t) formula is as follows:
Wherein γ is the consolidation constant.
Measurement model parameter θ (t)={ θ that we can obtain at yardstick level t
i(t), i ∈ I}, wherein θ
i(t)={ θ
K, i(t)=(μ
K,, i(t), σ
K,, i(t)), k ∈ Λ }.By to log p (y (t) | x (t), θ (t+1)) with respect to μ
K,, i(t) and σ
K,, i(t) ask partial differential, and to make it be zero.Can obtain μ
K,, i(t) and σ
K,, iBeing estimated as follows (t):
R wherein
k(t) all belong to the zone of organizing class k at yardstick level t in representative, i.e. the voxel location collection.
Step 4: local volume voxel classify again (PVVR)
Because the low local volume effect that spatial resolution caused of image may be by mis-classification by the voxel at cut zone edge.Shattuck and Noe are categorized as a kind of new class of organizing with local volume voxel when cutting apart conventional MR image, handle local volume effect.But this method is not suitable for cutting apart of DTI image.The DTI picture contrast is very low, even is difficult to distinguish white matter of brain and ectocinerea; And local volume voxel has and the pure similar intensity level of class of organizing usually.All these make to the conventional MR image of the ratio of division of DTI image difficult more.For addressing this problem, after utilizing self-adapting multi-dimension statistical classification (MSSC) classification, we detect inhomogeneous fringe region, reclassify local volume voxel, with further accurate segmentation result.We utilize the Canny edge detector, detect the fringe region of different tissues class.We utilize multiple dimensioned statistical classification (MSSC) method segmentation result, estimate the organization edge of possible original DTI image.Because different tissues class interior intensity value equates in the segmentation result, rim detection is easy to estimate between the class organizing.The voxel of fringe region is considered to possible local volume voxel, and these voxels are carried out subseries again with local volume voxel (PVVR) method of classifying again.Fig. 5 example after multiple dimensioned statistical classification (MSSC) classification the border condition between possible different tissues class.At Fig. 5 (a) (b) misclassification of local volume voxel all might take place under the situation.
We think that local volume voxel is the linear combined of the different tissues class that has been partitioned into,
θ={ θ
i, i ∈ I}, wherein θ
i=(μ
K,, i, σ
K,, i), k ∈ Λ }, represent known measurement model parameter.θ=θ (0) directly obtains from the segmentation result of multiple dimensioned statistical classification (MSSC).π
I, kRepresentative is at the possible locational weights of local volume voxel i, 0<π
I, k<=1.
Decide different classes shared weight on local volume voxel by finding the solution following formula.According to reality, we think that local volume voxel is by 2 inhomogeneity k
1, k
2Form.
θ
1 i=(μ
1, K, i, σ
1 K, i), k ∈ Λ, representative is at position i adjacent domain R
1 iIn organize mean intensity and the standard deviation of class k.
Wherein | R
1 K, 1| region R in the representative
1 iIn comprised belong to the total number of the voxel of organizing class k.For the unevenness of the intensity that overcomes original DTI image, use θ
1 iSubstitute θ
i
Because known image intensity distributions y, and can obtain measurement model parameter θ, the x that classifies, θ by multiple dimensioned statistical classification (MSSC)
1 i=(μ
1, K, i, σ
1 K, i), k ∈ Λ can obtain easily; By finding the solution π
I, k, reclassify local volume voxel then and be weight bigger organize class.
Operation result
Further verify our algorithm, we have chosen 20 cerebral ischemia patients' DTI image.These DTI images adopt GE 1.5T or 3.0T MR system, use the dispersion tensor imaging to obtain (TR/TE:6000-7000/98ms; Acquisition matrix: 128 * 128; The scan axis position; FOV:24cm; Bed thickness 5mm; Spacing 1.0mm; B value: 1000s/mm
2Disperse direction: 13 directions).In the acute or subacute headstroke stage, in the DTI image, the cerebral ischemia focal zone presents high signal.We are divided into three different class K=3 that organize with image: cerebrospinal fluid, white matter of brain and ectocinerea and cerebral ischemia focal zone.
Self-adapting multi-dimension statistical classification (MSSC)-local volume voxel (PVVR) method of classifying again effectively raises the accuracy of cutting apart of DTI image cerebral ischemia focal zone.We compare several dividing methods: the self-adaptation maximum a posteriori estimates that (MAP) cut apart, and self-adapting multi-dimension statistical classification (MSSC) is cut apart, self-adapting multi-dimension statistical classification (MSSC)-local volume voxel classify again (PVVR) cut apart.As Fig. 6, Fig. 7, shown in Figure 8.(Fig. 2, Fig. 3 and Fig. 6 all derive from same group of DTI data) here, the self-adaptation maximum a posteriori estimates that the initialization that (MAP) cut apart all adopts Threshold Segmentation result preferably; If initial segmentation is poor, the self-adaptation maximum a posteriori estimates that the effect that (MAP) cut apart will be worse than present effect far away.Yet, self-adapting multi-dimension statistical classification (MSSC) is cut apart with self-adapting multi-dimension statistical classification (MSSC)-local volume voxel (PVVR) dividing method of classifying again and is utilized regional split, merge algorithm to finish initial segmentation, the self-adaptation of comparing maximum a posteriori estimation (MAP) is cut apart, more robust, convenience.Can find out obviously that from experimental result self-adaptation maximum a posteriori estimation (MAP) is cut apart can not effectively solve the overlapping problem of intensity level; Self-adapting multi-dimension statistical classification (MSSC) is cut apart with self-adapting multi-dimension statistical classification (MSSC)-local volume voxel and is classified (PVVR) dividing method again owing to absorbed the anisotropic control of disperse in model, has well solved this problem; With the self-adaptation maximum a posteriori estimate (MAP) cut apart and self-adapting multi-dimension statistical classification (MSSC) dividing method relatively, self-adapting multi-dimension statistical classification (MSSC)-local volume voxel classify again (PVVR) cut apart the influence that has effectively reduced local volume effect.
Claims (5)
1. the automatic division method of a cerebral ischemia focal zone is cut apart DTI image cerebral ischemia focal zone automatically based on multiple dimensioned statistical classification and local volume sorting technique, it is characterized in that, may further comprise the steps:
(1) image pre-service is carried out filtering to original image;
(2) calculate the dispersion tensor field, obtains the pairing tensor field of each voxel of three dimensions;
(3) the anisotropic measurement of disperse quantizes the anisotropy of each voxel of three dimensions;
(4) calculating of metric space comprises original DTI graphical rule space, and disperse anisotropy mapping graph metric space;
(5) self-adapting multi-dimension statistical classification in the focus zone of seeking based on optimum, overcomes the influence that is brought owing to noise, disperse anisotropic and intensity non-uniformity;
(6) local volume voxel is classified again, on the basis of self-adapting multi-dimension statistical classification, further optimizes segmentation result, overcomes the interference that local volume effect brings.
2. the method for claim 1 is characterized in that, in (3) step, adopts the distribution anisotropy, quantitatively estimates the disperse anisotropy, and its formula is as follows:
Tr(D)=(λ
1+λ
2+λ
3)/3
Wherein Tr (D) is the mark of D, represents the average dispersity of each different directions.
3. the method for claim 1 is characterized in that, in (4) step, the spread function g () below adopting produces metric space:
Constant Ks produces with the manual method of setting or using Canny to propose.
4. the method for claim 1, it is characterized in that, in (6) step, utilize after the self-adapting multi-dimension statistical classification,, detect the fringe region of different tissues class with edge detector, use the segmentation result of multiple dimensioned statistical classification method again, estimate the organization edge of possible original DTI image, the voxel of fringe region is possible local volume voxel, to these voxels with local volume voxel again sorting technique carry out subseries again.
5. method as claimed in claim 4 is characterized in that, described edge detector is the Canny edge detector.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2003101211745A CN1305008C (en) | 2003-12-22 | 2003-12-22 | Automatic dividing method for cerebral ischemia focus area |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2003101211745A CN1305008C (en) | 2003-12-22 | 2003-12-22 | Automatic dividing method for cerebral ischemia focus area |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1632830A true CN1632830A (en) | 2005-06-29 |
CN1305008C CN1305008C (en) | 2007-03-14 |
Family
ID=34844083
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2003101211745A Expired - Fee Related CN1305008C (en) | 2003-12-22 | 2003-12-22 | Automatic dividing method for cerebral ischemia focus area |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1305008C (en) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010009596A1 (en) * | 2008-07-21 | 2010-01-28 | 深圳先进技术研究院 | Brain characteristic parameter obtaining method and system, thrombolysis indication generating system |
CN100595791C (en) * | 2008-06-18 | 2010-03-24 | 陕西高科医疗信息股份有限公司 | Mixed organization image full-automatic partition method of virtual colonoscope |
CN101067867B (en) * | 2006-05-05 | 2010-07-21 | 美国西门子医疗解决公司 | Diffusion tensor image visualization method and device |
CN101061951B (en) * | 2005-07-15 | 2010-10-13 | 美国西门子医疗解决公司 | Apparatus for classifying tissue using image data |
WO2010142085A1 (en) * | 2009-06-12 | 2010-12-16 | 深圳先进技术研究院 | Method and system for obtaining brain characteristic parameters, thrombolysis decision guideline system and method thereof |
CN101515364B (en) * | 2009-03-26 | 2011-05-04 | 哈尔滨工业大学 | Anisotropic picture flatting method based on most frequency value filtering |
CN102509123A (en) * | 2011-12-01 | 2012-06-20 | 中国科学院自动化研究所 | Brain functional magnetic resonance image classification method based on complex network |
CN101558321B (en) * | 2006-12-11 | 2014-10-29 | 皇家飞利浦电子股份有限公司 | Fibre tracking on the basis of macroscopic information |
CN104504657A (en) * | 2014-12-11 | 2015-04-08 | 中国科学院深圳先进技术研究院 | Method and device for de-noising magnetic resonance diffusion tensor |
CN104749540A (en) * | 2013-12-31 | 2015-07-01 | 西门子(深圳)磁共振有限公司 | Method and device for indicating inter-tissue distinction degree |
CN106780641A (en) * | 2016-11-14 | 2017-05-31 | 西安交通大学 | A kind of Low-dose X-ray CT image reconstruction method |
CN108665465A (en) * | 2018-05-14 | 2018-10-16 | 沈阳东软医疗***有限公司 | Divide method, apparatus, storage medium and the equipment of cerebral ischemic penumbra |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5526446A (en) * | 1991-09-24 | 1996-06-11 | Massachusetts Institute Of Technology | Noise reduction system |
JP2770849B2 (en) * | 1995-09-29 | 1998-07-02 | 日本電気株式会社 | Image processing device |
US5848193A (en) * | 1997-04-07 | 1998-12-08 | The United States Of America As Represented By The Secretary Of The Navy | Wavelet projection transform features applied to real time pattern recognition |
CN1184796C (en) * | 2001-07-26 | 2005-01-12 | 佳能株式会社 | Image processing method and equipment, image processing system and storage medium |
US7428338B2 (en) * | 2002-01-10 | 2008-09-23 | Ricoh Co., Ltd. | Header-based processing of images compressed using multi-scale transforms |
-
2003
- 2003-12-22 CN CNB2003101211745A patent/CN1305008C/en not_active Expired - Fee Related
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101061951B (en) * | 2005-07-15 | 2010-10-13 | 美国西门子医疗解决公司 | Apparatus for classifying tissue using image data |
CN101067867B (en) * | 2006-05-05 | 2010-07-21 | 美国西门子医疗解决公司 | Diffusion tensor image visualization method and device |
CN101558321B (en) * | 2006-12-11 | 2014-10-29 | 皇家飞利浦电子股份有限公司 | Fibre tracking on the basis of macroscopic information |
CN100595791C (en) * | 2008-06-18 | 2010-03-24 | 陕西高科医疗信息股份有限公司 | Mixed organization image full-automatic partition method of virtual colonoscope |
WO2010009596A1 (en) * | 2008-07-21 | 2010-01-28 | 深圳先进技术研究院 | Brain characteristic parameter obtaining method and system, thrombolysis indication generating system |
CN101515364B (en) * | 2009-03-26 | 2011-05-04 | 哈尔滨工业大学 | Anisotropic picture flatting method based on most frequency value filtering |
WO2010142085A1 (en) * | 2009-06-12 | 2010-12-16 | 深圳先进技术研究院 | Method and system for obtaining brain characteristic parameters, thrombolysis decision guideline system and method thereof |
CN102509123B (en) * | 2011-12-01 | 2013-03-20 | 中国科学院自动化研究所 | Brain function magnetic resonance image classification method based on complex network |
CN102509123A (en) * | 2011-12-01 | 2012-06-20 | 中国科学院自动化研究所 | Brain functional magnetic resonance image classification method based on complex network |
CN104749540A (en) * | 2013-12-31 | 2015-07-01 | 西门子(深圳)磁共振有限公司 | Method and device for indicating inter-tissue distinction degree |
CN104749540B (en) * | 2013-12-31 | 2017-11-17 | 西门子(深圳)磁共振有限公司 | The method and apparatus of discrimination between a kind of instruction tissue |
CN104504657A (en) * | 2014-12-11 | 2015-04-08 | 中国科学院深圳先进技术研究院 | Method and device for de-noising magnetic resonance diffusion tensor |
CN104504657B (en) * | 2014-12-11 | 2017-05-03 | 中国科学院深圳先进技术研究院 | Method and device for de-noising magnetic resonance diffusion tensor |
US10620287B2 (en) | 2014-12-11 | 2020-04-14 | Shenzhen Institutes Of Advanced Technology Chinese Academy Of Sciences | Method and apparatus for denoising magnetic resonance diffusion tensor, and computer program product |
CN106780641A (en) * | 2016-11-14 | 2017-05-31 | 西安交通大学 | A kind of Low-dose X-ray CT image reconstruction method |
CN106780641B (en) * | 2016-11-14 | 2020-07-28 | 西安交通大学 | Low-dose X-ray CT image reconstruction method |
CN108665465A (en) * | 2018-05-14 | 2018-10-16 | 沈阳东软医疗***有限公司 | Divide method, apparatus, storage medium and the equipment of cerebral ischemic penumbra |
Also Published As
Publication number | Publication date |
---|---|
CN1305008C (en) | 2007-03-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11250569B2 (en) | Systems and methods for functional imaging follow-up evaluation using deep neural network | |
Van Rikxoort et al. | Automatic lung segmentation from thoracic computed tomography scans using a hybrid approach with error detection | |
US6909797B2 (en) | Density nodule detection in 3-D digital images | |
EP2989609B1 (en) | Processing imaging data to obtain tissue type information | |
US20080292194A1 (en) | Method and System for Automatic Detection and Segmentation of Tumors and Associated Edema (Swelling) in Magnetic Resonance (Mri) Images | |
US20080050020A1 (en) | Density nodule detection in-3-D digital images | |
US20150045651A1 (en) | Method of analyzing multi-sequence mri data for analysing brain abnormalities in a subject | |
WO2007044508A2 (en) | System and method for whole body landmark detection, segmentation and change quantification in digital images | |
CN1305008C (en) | Automatic dividing method for cerebral ischemia focus area | |
EP2987114B1 (en) | Method and system for determining a phenotype of a neoplasm in a human or animal body | |
CN1918601A (en) | Apparatus and method for registering images of a structured object | |
Li et al. | Robust unsupervised segmentation of infarct lesion from diffusion tensor MR images using multiscale statistical classification and partial volume voxel reclassification | |
Ferjaoui et al. | Machine learning for evolutive lymphoma and residual masses recognition in whole body diffusion weighted magnetic resonance images | |
Eskildsen et al. | Detecting Alzheimer’s disease by morphological MRI using hippocampal grading and cortical thickness | |
EP3884462A1 (en) | Characterizing lesions in radiology images | |
Perasso et al. | Spinal canal and spinal marrow segmentation by means of the Hough transform of special classes of curves | |
Ulloa et al. | Improving multiple sclerosis lesion boundaries segmentation by convolutional neural networks with focal learning | |
Mihaylova et al. | A brief survey of spleen segmentation in MRI and CT images | |
Cárdenas-Peña et al. | Adaptive Bayesian label fusion using kernel-based similarity metrics in hippocampus segmentation | |
Zanaty et al. | Improving fuzzy algorithms for automatic magnetic resonance image segmentation. | |
Cherradi et al. | Fully automatic method for 3D T1-weighted brain magnetic resonance images segmentation | |
Palanikumar et al. | Automatic detection of solitary pulmonary nodules using superpixels segmentation based iterative clustering approach | |
Kim et al. | Automatic 3D temporal kinetics segmentation of dynamic emission tomography image using adaptive region growing cluster analysis | |
Fenneteau et al. | A size-adaptative segmentation method for better detection of multiple sclerosis lesions | |
CN117649400B (en) | Image histology analysis method and system under abnormality detection framework |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20070314 Termination date: 20181222 |