Embodiment
The present invention relates to be used to detect and mate the method for anatomical structure.To adopt an example of this method to be: by detecting via machine learning or classification and cut apart endocardium of ventricle and epicardial border, and to carry the similar situation of the database of note, detect the zone walls dyskinesia in the heart by identification.Those of ordinary skill in the art it should be understood that the present invention can be used to other and use, and wherein SHAPE DETECTION and coupling are useful, such as but be not limited to identification people's feature, such as facial characteristics or other physical trait.The present invention also can be used to two dimension, the three peacekeeping four-dimension (3D+ time) data analysis, for example such as past that can be in time and the medical analysis of the anatomical structure the heart, lung or the tumour that develop.
In order to describe the present invention, will example be described at the intracardiac wall of the left ventricle that detects human heart.Fig. 1 illustrates the exemplary architecture of ultrasonic cardiography drawing system, and wherein said ultrasonic cardiography drawing system uses according to of the present invention and is used to utilize shape and outward appearance to detect the method for the intracardiac wall of left ventricle.Medical sensor 102 such as ultrasonic transducer is used to patient is carried out inspection.Sensor 102 is used to obtain and the concrete consistent medical measurement of medical inspection.For example, the patient who stands cardiac problems can have the echocardiogram of carrying out for the concrete heart disease of assisted diagnosis.Ultrasonic system provides two, three and four (3D+ times) the dimension image of heart according to various skeleton views.
The information that obtains by sensor 102 is transferred into processor 104, and described processor can be workstation or personal computer.Processor 104 converts sensing data to the image that is transferred into display 108.Display 108 also can transmit other graphical information that relates to image or the table of information.According to the present invention, also provide the data of the initial profile of the intracardiac wall of expression for processor 104.These data can be provided by the user artificially such as doctor or sonographer, or are automatically provided by processor 104.This profile comprises a series of independent points, is illustrated on the display 108 by the motion of processor 104 these points of tracking and with it.
Except the data from medical sensor 102, processor 104 also can receive other data input.For example, processor can receive the data from the database 106 that is associated with processor 104.This data can comprise the subspace model, the potential contour shape of the intracardiac wall of described subspace model representation.These subspace models can be the images of a plurality of patients' of expression left ventricle, perhaps can be the contour shape models based on statistical information that computing machine produces.Utilize known method, such as Bayes nuclear (Bayesian kernel) coupling or based on the method for light stream, processor 104 is followed the tracks of the independent point of contour shapes.Make error accumulation in tracing process up by utilizing multi-template Adaptive matching framework.The form of sentencing covariance matrix at each point is represented the uncertainty of following the tracks of, and wherein uses the subspace shape constraining of nonopiate projection to adopt described covariance matrix subsequently fully.
Fig. 2 illustrates the typical echocardiogram image of heart.Come label that the part of intracardiac wall of the left ventricle of acoustics signal drop-out is arranged with solid line ellipse 208.With dashed lines ellipse 202,204 is represented the estimation of local wall motion.Because acoustic signal is lost, intracardiac wall is the strongest edge in image always not.The feature of echocardiogram image is image fan-shaped of with dashed lines 210,212 expression.The data that do not include usefulness in the zone of fan-shaped outside.
Many detection methods utilize the lifting of Weak Classifier or feature to come object in the detected image.Weak Classifier by making up selected quantity is via the response that promotes, and resulting strong classifier can be realized the high detection rate.Yet known method does not solve the problem of detected object when having other occlusion objects (data outside for example fan-shaped).Weak Classifier is because the errored response of blocking influences the detection of object negatively.
According to one aspect of the present invention, use description in the object detection process, eliminate the method for known influence of blocking now.For example, can handle the echocardiogram image in the mode of not considering the view data (just useless or invalid data) outside fan-shaped.In other words, the data in fan-shaped outside are used as and block processing.
The simple feature that is associated with the image of object is identified as Weak Classifier.The example of this feature is in the rectangular characteristic shown in Fig. 3 a-3d.The value of each rectangular characteristic is poor between the summation of white (the also just being known as) zone of each rectangle and the pixel intensity in grey (also the being known as negative) zone.For the rectangular characteristic shown in Fig. 3 a, negative region is 302, and positive region is 304.For the rectangular characteristic shown in Fig. 3 b, negative region is 308, and positive region is 306.For the rectangular characteristic shown in Fig. 3 c, negative region is 312 and 314, and positive region is 310 and 316.For the rectangular characteristic shown in Fig. 3 d, negative region is 320, and positive region is 318 and 322.
Rectangular characteristic provided complete basis for the fundamental region.For example, if rectangle is 24 * 24 pixels dimensionally, then the quantity of feature is 180000.One of advantage of rectangular characteristic is a computing velocity.By utilizing the intermediate representation that is called integral image (II) as shown in Figure 4, can utilize the operation of little fixed qty to come computation of characteristic values.
Before the calculating of rectangular characteristic, the II of calculating input image (for example echocardiogram image of left ventricle) in advance.(x y), determines brightness value at each pixel among the II.These brightness values are stored in (Fig. 1) in the database 106.In case, just simplify the calculating in all futures greatly for input picture has calculated II.At the position (x in the input picture
0, y
0) each pixel of locating, can be by determining at position (x
0, y
0) on and at position (x
0, y
0) the brightness sum of all pixels on the left side calculate brightness value.In other words, at II (x
0, y
0) locate can following definite II subclass:
Wherein (x y) is in the position (x, the brightness of the pixel of y) locating to I.
Fig. 3 illustrates and how to determine at rectangular characteristic R
fThe calculating of the brightness value of the II of place.The II at calculating location 408 places, it equals the zone within the solid line 410.The another way that is used to limit the II at 408 places, position is rectangle (A+B+C+R
f) the summation of brightness value.In order to obtain R
fSummation, must carry out additional calculations.The II of position 406 provides the summation in the zone that limits by line 412, and it equals the summation of the brightness value of rectangle (A+C).The II that deducts position 406 from the II of position 408 causes rectangle (B+R
f) II '.Then, the II of calculating location 404, it provides the summation in the zone that limits by (A+B).The II that deducts position 404 from II ' causes rectangle (A+R
f) II ".At last, " addition, this provides R with the II and the II of position 402
fSummation.
Yet, at R
fIn pixel comprise that in the situation of blocking, the brightness value of those pixels provides invalid value, this will finally produce the incorrect estimation to rectangular characteristic.Fig. 5 illustrates and comprises the example that blocks 504 integral image 502.Rectangular characteristic 506 is placed in and comprises the position of blocking a part of 504.
According to the present invention, block masks is used to eliminate the effect that is included in the pixel that is blocked in the rectangular characteristic.Figure 6 illustrates the example of the block masks of the II that is used for Fig. 5.When in check environment, obtaining image, block masks can be used, perhaps block masks can be from data, inferred.For example, in surveillance application, known quiescent state background (for example position of door, wall, furniture etc.).Can determine to cause the likelihood of the object in the background of blocking, and use it for the establishment block masks.Another example is a ultrasonoscopy.In ultrasonoscopy, provide fan-shaped position by ultrasound machine, maybe can calculate fan-shaped position, for example the analysis of time variation can produce static inactive area.In case identify fan-shapedly, just can create block masks, in II calculates, to get rid of effectively or to cancel fan-shaped existing.
By will be blocked or otherwise the brightness value of invalid pixel is made as zero, the summation of the brightness value of rectangle will be influenced by incorrect value no longer.Yet, because there are " losing " data now, so summation will be unbalanced.When not existing when losing value, the average brightness value of rectangle summation and rectangle is proportional.Therefore, for the compensating missing value, when existence was blocked, the quantity that has the pixel of effective brightness value by utilization was come approximate average.Can be by at first calculating the quantity that isoboles or block masks obtain valid pixel.
Block masks M comprises Boolean, and wherein the valid pixel value of being endowed 1, and pixel value of being endowed 0 invalid or that be blocked.Can following utilization at current location (x
0, y
0) on and at current location (x
0, y
0) the quantity of valid pixel on the left side calculate the integration mask:
Be similar to the II of equation (1), can calculate the quantity of valid pixel in the rectangle according to the integration mask with the operation of aforesaid equal number.
Weighted difference between the summation of the brightness in the positive and negative image-region will be provided for the equivalent features value of rectangular characteristic 506.If R
+Remarked pixel brightness is with on the occasion of zone of making contributions and R
-The zone that remarked pixel brightness is made contributions with negative value, then eigenwert f is as follows:
N wherein
-, n
+The quantity of representing the valid pixel in negative, positive zone respectively, each zone comprise N pixel.If n
-And n
+All non-zero then passes through N/ (n
-n
+) with final eigenwert standardization.By utilizing block masks to calculate the integral image of rectangular characteristic, obtain more accurate result, it causes better object detection.
Particularly under the situation such as the complex object of face or anatomical structure, owing to the detection for object needs calculated big measure feature or part, instrument is used to reduce needed calculated amount, still produces accurate result simultaneously.Generally a kind of such instrument of Shi Yonging is to promote.Generally, promote a plurality of Weak Classifiers of identification or feature.At each Weak Classifier, can calculate a value, described value compares with predetermined threshold then.If the value of Weak Classifier surpasses threshold value, then keep this sorter.If the value of Weak Classifier is lower than threshold value, then refuse this sorter.By asking weighted sum to surpassing all values threshold value, Weak Classifier, can produce strong classifier, it can be used in the object detection.
The modification that promotes is the cascade that promotes.In this technology, sorter is distinguished priority ranking.Calculate first sorter at window, and if it does not satisfy threshold value, window is moved into the another location so.Only keep those positions that the sorter that is calculated surpasses threshold value.Threshold value typically is set at the level place of appropriateness, to allow roomy margin for error.The quantity of the position by reducing to carry out calculating, this method is effective.Yet these methods abandon the output from before sorter.The training of next stage begins with the homogeneous weighting to new sample group.
According to the present invention, keep the value of each sorter that calculates, and use it for the calculating of following sorter so that strengthen after the classifier calculated in stage.By on the new training group of the current generation as shown in Fig. 7 a and 7b, new threshold value T being set
i *And relevant parity p
i *, the stage is directly used centre strong classifier H before the cascade
I-1Then, before training, with new sorter H based on primitive character
I-1 *Error training sample is weighted.The strong classifier of current generation is H
I-1 *Weighted sum with selected single tagsort device.In testing process, replace throwing away from the sorter output in stage before, it is utilized T
i *With relevant parity p
i *Threshold value is set, and is weighted, and be added in the weighted sum from the output of single feature Weak Classifier of current generation according to its error.
" utilizing the cascade of the lifting of storer " can describing as getting off be training algorithm (BCM):
Initial cascade classifier in the storer: H
0 *=NULL; Error ε=∞ wherein.
The positive sample group of P=, N=negative sample group;
·i=0;
The circulation of stage i by cascade:
1++;
Use P and N, utilize H to use AdaBoost
I-1 *Come training classifier H with the supplementary features of quantification
i, so that satisfy desired wrong positive number and false detection rate;
By with H
iWrong positive number as negative sample and regulate in case of necessity and just organizing reconstruct training group, P → P
*And N → N
*
Again train H
i: revise H
i(for example by selecting optimal threshold T
i *With parity p
i *), with according to new training data P
*And N
*Minimize its error in classification.With resulting H
i *Put into storer with its output.
In the superincumbent algorithm, suppose if use H
I-1 *, then it will be first.H
I-1 *Also can be used in the centre of other single feature Weak Classifier.What all need change is, replaces when the end in stage before, but in the training process of current generation, learns new threshold value and parity.
Can utilize current training group or complete representational checking to organize and carry out H
iAssessment.In last situation, only use H
i, and target is to satisfy the performance objective of current generation (that is per stage 95%); And in one situation of back, should use total cascade classifier, and target is the Comprehensive Performance Objective (that is, 0.95 ') until current generation i.
Another example of the cascade algorithm of the lifting that utilizes storer is shown below.
The user selects f, is every layer of acceptable false positive rate of maximum and d, is the value of every layer of acceptable verification and measurement ratio of minimum.
The total false positive rate F of user's select target
Target
The positive sample group of P=, N=negative sample group
·F
0=1.0;D
0=1.0;i=0;
Initial cascade classifier in the storer: H
0 *=NULL;
Work as F
i>F
TargetThe time
i++
n
i=0;F
i=F
i-1
Work as F
i>f * F
I-1The time
οn
i++
ο uses P and N, utilizes H to use AdaBoost
I-1 *And n
iIndividual supplementary features are come training classifier H
i(when i>1, utilized the result who is stored in the storer to train H
I-1 *If it is stronger than feature, then select H
I-1 *Weak Classifier as stage i.Utilize H
I-1 *Error upgrade weight about training sample so that based on the remaining Weak Classifier of feature selecting.)
The ο assessment is about the current cascade classifier H of checking group
i, to determine F
iAnd D
i
ο reduces i sorter H
iThreshold value (T
i), have d * D at least up to current cascade classifier
I-1Verification and measurement ratio (this also influences F
i)
N=
If F
i>F
Target, assess current cascade detectors H so about non-face image sets
i, and any error-detecting put into the group N in.
Again train H
i *: be H
iSelect the second optimal threshold T
i *With parity p
i *, minimize with the error in classification that makes it according to new training data P and N; With resulting H
i *Put storer into its output.
In the BCM framework, testing process is used for strong classifier in the middle of following each stage of cascade:
Stage 1-H
1: α
1h
1+ α
2h
2+ ... α
N1h
N1
Stage 2-H
2: α
H1p
1 *(H
1-T
1 *)+(α
N1+1h
N1+1+ α
N1+2h
N1+2+ ... α
N1+n2h
N1+n2)
Because in the stage before, assessed H
i, additional calculating only is subtraction and the multiplication at each additional phase, wherein current as shown in Figure 8 specimen will be by described additional phase.For each additional category device of being considered, consider the relevant value of sorter with previous calculating, and it is combined with value at added value.Resulting value compares with threshold value then.Net result provides the value of the more accurate indication of object detection.
In case detected object potentially, just can utilize further treatment technology to obtain additional information about image.According to another aspect of the present invention, can be based on the common use of one group of training image application appearance and shape, with object shapes or the anatomical structure in coupling and the detection test pattern.Outward appearance is used to the object in the test pattern or the location of structure.Adopt matching technique then, find similar situation and to provide shape or CONSTRUCTED SPECIFICATION for the material standed for that is detected to concentrate from positive training data.Can detect, locate and mate with hierarchical approaches, to realize more accurate result.The feature that the matching technique utilization is learnt in testing process is to save computing time and to improve matching performance.
According to the present invention, be used for the general framework that object or anatomical structure detect and shape is recovered and comprise three phases: off-line training step, online detection-phase and matching stage.
In off-line training step, positive training sample is stored as the training data group.For example, under echocardiographic situation, the training group will comprise the image of the left ventricle of human heart.The training group will comprise comprehensive group of sample of the sample of differing formed left ventricle and unusual left ventricle.Preferably, (for example left ventricle occupy in the middle of the image in top condition; Image by standardization with the influence of eliminating size and rotation etc.) under image in the training group is shown.In the training stage process, boosting algorithm is applied to image.Use the data that obtained from these and be stored, and can comprise the proper vector of indication Weak Classifier output with training data.
Handle all positive training samples, with the characteristic that remains unchanged.For example, integral translation, rotation and convergent-divergent are the constant conversion (invariant transform) to the left ventricle of human heart.The aligning of correction data directly influences the design of detecting device, that is, each aiming axis need be expanded in testing process.In other words, for example, if cancellation rotation in training data, then detecting device has to search for a plurality of rotations in testing process.The training data that utilization is aligned, learning algorithm are exported selected feature and just are being used for/are bearing the corresponding decision function of classification.According to the present invention, can all training datas of conversion (for example convergent-divergent and rotation), with the detecting device of training institute conversion.
All features of training data (comprising the training data that is transformed) are stored in (Fig. 1) in the database 106.In addition, be each correction data image calculation proper vector, described correction data image comprises Weak Classifier output and their associated weight.Each Weak Classifier is represented the part relevant with object.The proper vector of correction data subsequently can with compare for proper vector that test pattern calculated, have the correction data of similar features to help identification.The location point corresponding to along the point of object outline of each correction data image also is stored, and is used to the form fit stage.
In online detection-phase, by window in test pattern or the data volume or cube are carried out translation, rotation and/or convergent-divergent, be used for the sweeping scheme of image or data volume, to produce candidate data sticking patch (patch).At each material standed for, during whole scanning process or reach object's position, convergent-divergent and/or rotation afterwards.In some cases, a plurality of detecting devices that are transformed of application are faster than conversion material standed for.Can be used to object or material standed for position in the indicating image based on the detecting device that promotes.
In stage, identification appears to comprise those material standed fors (material standed for of winning victory) of object in form fit.The similarity matching algorithm is applied to these material standed fors,, and their associated shape is applied to material standed for the nearest neighbor of retrieval from corresponding training data group.Proper vector is used to form fit, and described form fit is based on Weak Classifier output h/s and their associated weight α/s.Utilize the distance metric in the following space to mate: { α
1h
1, α
2h
2..., α
Kh
K, wherein K is the quantity of Weak Classifier feature.Also can use the further feature space." pairing (counterpart) " proper vector that the form fit algorithm search is associated with one or more images in the training data.In case obtained coupling, come the data of the relevant matches image in the self-training group can be used to provide details about object shapes and structure.
To utilize two dimension (2D) data to describe example now.Those skilled in the art will appreciate that the present invention also can handle three-dimensional (3D) data and the four-dimension (3D+ time) data.The exemplary framework that Fig. 9 illustrates the data set that utilizes the original image and the endocardial boundary of being followed the tracks of, use the band note of LV frame, detects based on left ventricle (LV) endocardial boundary of study and coupling.
This framework is at first constructed the database of the LV sticking patch orderly group of form, that be aligned (modulus translation, convergent-divergent and rotation) at boundary mark for example or reference mark, wherein its boundary strip note.Utilize this group sample and other negative sticking patch to design and train detection algorithm.For example, can use based on the feature selecting and the sorter construction algorithm that promote.Be all sample calculation proper vectors.
Then, this framework is the detection algorithm of the location employing study of LV.The material standed for that is detected is to comprise local sticking patch LV, that have its correct size and rotation.(by determining this size and rotation with a plurality of possible convergent-divergents and rotation sweep image.)
At last, from the material standed for that is detected, extract proper vector, and compare, to find one or more nearest neighbor with database.For example utilize weighted sum to make up their profile then, to form the profile of the candidate's sticking patch detected, wherein the matching distance of the sticking patch that extremely detected with them of weight is proportional.A possible characteristic type will be the brightness of image (having or do not have double sampling) in the sticking patch.Further feature comprises Weak Classifier output (having or do not have weighting factor).Also can utilize principal component analysis (PCA) (Principal ComponentAnalysis), to select the subclass of this feature.
Figure 10 illustrates the example of the lesion detection in the three-dimensional data volume.This method is similar to method recited above, except 3D neighborhood, 3D feature and 3D scanning are used to detect coupling.
In case detecting aspect shape and the outward appearance and mating object, then just can be along with the past of time the shape of tracing object.Because the rhythmic motion of cardiac muscle, this tracking is important in echocardiogram.In the shape tracing process, measuring uncertainty plays an important role.According to the present invention, with separately the training the part detector applies in image, to develop metastable local appearance, utilize the global shape model to come the constraint portions merging process simultaneously.Use description in automatic shape detection and tracking process, merge best from local detection now, move dynamically and probabilistic Unified frame of subspace shape modeling.The part detecting device that promotes is used to the left ventricle boundary alignment in the ultrasonic cardiography graphic sequence.
Suppose (x, C with N
x), promptly have mean value x and a covariance C
xThe candidate that represents of multidimensional Gaussian distribution before the detection of shape, first step be the situation of the constant conversion of the best following before sample shape x
0Among find and have (x, C by N
x), shape N (m, C
m) and from the predicting shape N (x of previous time step
-, Cx
-) shape before the sample of the common PRML that produces.The equivalence formula is to find x
*, with the summation of the Mahalanobis distance in the shape space of shape space before minimizing and institute's conversion, just,
X wherein
0 '=T (x
0), wherein T is constant conversion.
Under the situation of shape before a plurality of candidates, produce the highest likelihood, also consider to detect the candidate of the likelihood value among the figure before shape win in the judgement time.Equation (5) need be in position and conversion in optimization, and even do not have separating of closing form for simple conversion, such as the similarity conversion that only allows translation, rotation and convergent-divergent yet.Can numerically seek whole the best by iteration, may be too expensive but calculate.
Difficulty comes from the following fact, and promptly the stream shape (shape just) of being crossed over through all possible conversion by shape before arbitrarily generally, particularly the group space dimensionality did not intersect with the shape subspace in relative hour.In the present invention, the shape subspace has from 6 to 12 dimension, and full Euclid (Euclidean) space has 〉=and 34 dimension.Figure 11 illustrates the invariant manifold of the shape alignment that is used for illustrating conceptually this relation, the stream shape that the one dimension Gaussian distribution 1106 of shape vector X and sloping shaft 1104 and expression subspace model was crossed over before wherein bold curve 1102 was described.Generally, stream shape will not given space intersection (just sloping shaft 1104 comprises model barycenter M) with shape.Omit prediction at this, perhaps can regard X as detect and predict amalgamation result.The present invention relates to two step optimization approach of solution as a whole, have separating at the closing form of two steps.Can easily explain this scheme with reference to Figure 11: first step is to utilize C
xIn information go to X from X
*, or in other words, find optimal mapping from X to M.Second step is to be used to from C
MAdditional information from X
*Go to X
MFirst step is known as alignment procedures, and second step is known as the constraint step.
The target of alignment procedures is to consider in the component uncertainty in the process of model transferring of shape and its covariance matrix before.The at first following d that minimizes
2:
Wherein x '=T (x), and C '
x=T (C
x).For contracted notation, suppose prediction N (x
-, C
x) be merged into N (x, C
x).
When T is the similarity conversion, have:
x′=Rx+t, (7)
Wherein t is the translation vector with two free parameters, and R is block diagonal matrix, and wherein each piece is
Utilize pure algebraically (straight algebra), can following rewrite equation (6):
d
2=(R
-1(m-t)-x)
TC
x -1(R
-1(m-t)-x)
=(T
-1(m)-x)
TC
x -1(T
-1(m)-x) (9)
By to four free parameter differentiations among R and the t, can obtain separating of closing form.Figure 12 a and 12b are illustrated in the shape alignment under probabilistic situation of considering and not considering the some position.Figure 12 a is illustrated in the shape alignment under probabilistic situation of not considering the location.Figure 12 b illustrates and utilizes the probabilistic shape alignment of different variance.Oval 1202-1212 describes the covariance about a position, expression piece diagonal angle C
xIn information.Intuition is more to believe to have the more point of high confidence level.
In case shape is aimed at model before, the shape with PRML that is produced by two competition information sources (just, the detection/prediction of aligning is relative with (subspace) model) is determined.Under the situation of total space model, formula directly relates to the information of utilizing Gaussian source or BLUB (best linear unbiased estimator (Best Linear Unbiased Estimator)) and merges.
Suppose two noisy measurements of identical n dimension variable x, each measurement is characterised in that the multidimensional Gaussian distribution, i.e. N (x
1, C
1) and N (x
2, C
2), the maximal possibility estimation of x is the Mahalanobis distance D with revision
2(x, x
2, C
2) the point of minimum summation.Suppose C with being without loss of generality
2Be unusual.Utilize C
2=U Λ U
TSvd, U=[u wherein
1, u
2..., u
n], while u
iBe orthonormal, and Λ=diag{ λ
1, λ
2..., λ
p, 0 ..., 0} is to x
2Mahalanobis distance as follows:
Work as λ
iBe tending towards at 0 o'clock, D
2(x
1, x
2, C
2) the trend infinity, unless U
T 0X=0, wherein U
0=[u
P+1, u
P+2..., u
n].Suppose the initial point of subspace here by luv space with being without loss of generality.Because x
2Occupy in the subspace, so U
0 Tx
2=0.
Because U
0 TSo x=0 is d
2Become now:
Wherein y is 1 * p vector.
To the y differentiation, produce the merging estimator of subspace:
Utilize the equivalent expression in the luv space:
C can be shown
X*And C
Y*Be x
*And y
*Corresponding covariance matrix.
Alternatively, can write equation (12) and equation (13) as getting off:
Here, y
2Be x
2By U
pCoordinate transforming in the subspace of being crossed over, and Λ
p=diag{ λ
1, λ
2..., λ
p.The BLUE that equation (16) can be counted as in the subspace of two Gaussian distribution merges, and one of them Gaussian distribution is N (y
2, Λ
p), and another Gaussian distribution is the N (x in the subspace
1, C
1), N ((U
T pC
1 -1U
p)
-1U
T pC
1 -1x
-1, (U
T pC
1 -1U
p)
-1) common factor.
Above-mentioned subspace merge into (subspace) model constrained, will (it is probabilistic to have different variance) shape measure and principal component analysis (PCA) (PCA) shape be considered as two information sources general formula be provided.Below, add the 3rd source, its expression is according to the performance prediction of following the tracks of.The vital benefit that obtains from tracking except that detecting is additional information and the information merging in time from the system dynamics of domination prediction.Based on top analysis, separating of equation (4) has following form:
This is separated the information from detection, shape and performance prediction is placed in the united frame.When also limiting predicting shape in the subspace, subspace recited above BLUE formula can be applied in the conversion T with nested mode.Prediction N (x
-, C
X-) comprise information from system dynamics.This information is used to global motion trend, such as expansion with shrink and translation slowly and rotation are encoded.Utilize traditional method, the predictive filter in being provided with such as Kalman, can obtain N (x
-, C
X-):
Wherein the system dynamics equation is
x
-=Sx
+,prev+q,
(20)
And Q is the covariance of q, and " prev " expression is from the information of previous time step.
Figure 13 illustrates the synoptic diagram of analytical procedure, wherein propagates the uncertainty of detection in steps by institute.At each frame place,, a plurality of detection material standed fors are assessed by the likelihood of comparison a plurality of detection material standed fors in the context of two shapes and based on the prediction of system dynamics according to former frame.Ellipse such as 1302-1316 illustrates the position uncertainty.In alignment procedures, utilize shape that uncertainty is carried out conversion, in likelihood estimation and tracing process, utilize the prior imformation of model and prediction to merge uncertainty.
Described and be used to utilize outward appearance and shape to detect and mate the embodiment of the method for anatomical structure, but it may be noted that according to top instruction, those skilled in the art can make amendment and change.Therefore it should be understood that and in disclosed specific embodiments of the invention in the scope and spirit of the present invention as defined by the appended claims, to change.Therefore described the present invention, set forth the content of and expectation protection claimed in the appended claims by patent with the desired details of Patent Law and characteristic.