Abstract:
The invention relates to a method for evaluating treatment-relevant spatial anatomical information among different data sets of the heart, the method comprising the steps of:—determining a reference anatomical 3 dimensional data set of the heart,—providing a first anatomical 3 dimensional data set of the heart, the first anatomical 3 dimensional data set comprising first treatment-relevant spatial anatomical information,—providing a second anatomical 3 dimensional data set of the heart, the second anatomical 3 dimensional data set comprising second treatment-relevant spatial anatomical information,—registering the reference data set to the first and the second data sets,—transferring the treatment relevant spatial anatomical information of the first and the second data set to the reference data set in order to generate a first transferred treatment-relevant spatial anatomical information on the reference data set and a second transferred treatment-relevant spatial anatomical information on the reference data set—evaluating the first and the second transferred treatment-relevant spatial anatomical information.

Description:
[0001]    The invention relates to a method for evaluating treatment-relevant spatial anatomical information among different data sets of the heart and to a system therefor 
         [0002]    Statistical shape models are widely used today in different fields of medical image processing. A common application is the use as prior information for segmentation of 3-D medical image data [1]. The left atrium is a challenging structure, as it shows a large amount of variation in surface topology and shape across different patients. In addition to the anatomical differences among subjects, there is the aspect of cardiac motion which leads to deformations of the left atrium in different heart phases. In electrophysiology ablation procedures, a model of the underlying anatomical structure could be used for planning of the intervention as, e.g., suggested by Keustermans et al. using patient-specific 3-D data sets for planning of atrial fibrillation treatment [2]. 
         [0003]    Since electrophysiology (EP) procedures involve fluoroscopic guidance, C-arm X-ray devices with image overlay functionality, fusing soft-tissue information with live fluoro images, offer an alternative navigation approach to mapping systems [3, 4]. For EP procedures, these systems can be used to merge soft-tissue heart information derived from pre-procedurally acquired 3-D data sets as well as planning information with live fluoroscopy during the intervention [5]. A first approach for planning of atrial fibrillation has been presented by Hastenteufel et al. [6]. More recently, Brost et al. [7] proposed a system for interactive planning of EP procedures, especially atrial fibrillation ablation procedures involving cryo-balloons. An example of pre-planned ablation lines for PVI is shown in  FIG. 1 . Ablation lines are one example for treatment-relevant spatial anatomical information. The use of pre-planned ablation lines is not restricted to fluoro overlay navigation. This kind of annotation can also be helpful when performing PVI using a mapping system, e.g. CARTO 3 (Biosense Webster, Diamong Bar, Calif., USA) or Ensite Velocity (St. Jude Medical, St. Paul, Minn., USA). 
         [0004]    Pulmonary vein isolation is the treatment of choice for paroxysmal atrial fibrillation with a reported success rate of over 70% [8]. PVI is also applied to persistent AFib, however with lower success rate of just over 50%. The structure of the left atrium is very complex. There are critical areas within the left atrium, e.g. the left atrial appendage or the area close to the esophagus, that should be avoided or treated with special care during ablation procedures. There are different strategies to perform PVI. Two common strategies are segmental ablation and ipsilateral pulmonary vein (PV) isolation. Arentz et al. investigated the influence of the isolation area around the PVI on the procedure outcome. Based on their findings, the ipsilateral ablation strategy is advantageous [9]. 
         [0005]    Personalized planning of ablation lines for electrophysiology procedures offers the potential to improve procedure success while lowering the risk to the patient. As of today, there has, however, not been an investigation on how ipsilateral ablation lines differ across similar cases. 
         [0006]    Thus, the need exists to further improve the evaluation of treatment-relevant spatial anatomical information of ablation lines. This need is met by the features of the independent claims. Further embodiments are described in the dependent claims. 
         [0007]    According to a first aspect, a method for evaluating treatment-relevant spatial anatomical information among different data sets of the heart is provided. In one step, a reference anatomical three-dimensional data set of the heart is determined. Furthermore, a first anatomical three-dimensional data set of the heart is provided, which comprises first treatment-relevant spatial anatomical information. A second anatomical three-dimensional data set of the heart with second treatment-relevant spatial anatomical information is provided. The reference data set is registered to the first and the second data set. In a further step, the treatment-relevant spatial anatomical information of the first and of the second data sets are transferred to the reference data set in order to generate a first transferred treatment-relevant spatial anatomical information on the reference data set and a second transferred treatment-relevant spatial anatomical information on the reference data set. Additionally, the first and the second transferred treatment-relevant spatial anatomical information is evaluated. 
         [0008]    As the treatment-relevant spatial anatomical information ablation planning lines of the left atrium may be used. 
         [0009]    The step of evaluating the first and the second transferred treatment-relevant spatial anatomical information can comprise the step of determining at least one of a mean value or standard deviation of the first and the second transferred treatment-relevant spatial anatomical information on the reference data set. Furthermore, the step of evaluating the transferred treatment-relevant spatial anatomical information may contain the step of indicating both pieces of information on the reference data set, e.g. the indication of both transferred ablation lines on the reference data set. 
         [0010]    According to a further aspect, a plurality of further anatomical three-dimensional data sets of the heart may be provided with corresponding treatment-relevant spatial anatomical information. The reference anatomical three-dimensional data set can be registered to each of the further data sets and the treatment-relevant spatial anatomical information of each of the further data sets is transferred to the reference data set in order to determine transferred treatment-relevant spatial anatomical information on the reference data set for each of the further data sets. 
         [0011]    The part of the heart that is of special interest in the present invention can be the left atrium of the heart. 
         [0012]    Furthermore, it is possible to determine an average anatomical three-dimensional data set of the atrium based on at least the first and second data sets and an average treatment-relevant spatial anatomical information can be determined on the reference data set taking into account at least the first and the second transferred treatment-relevant spatial anatomical information. 
         [0013]    Preferably the reference anatomical three-dimensional data set, the first and second data set and the further anatomical three-dimensional data sets of the atrium are each provided as a mesh structure, and a non-rigid mesh registration is used to register one data set to the other data sets. One possible registration method is the coherent point drift method. However, any other registration method may be used. 
         [0014]    The invention may also be used for planning an interventional treatment. A new three-dimensional data set of the atrium is provided and registered to the reference data set in order to determine treatment-relevant spatial anatomical information for the new data set based on the transferred treatment-relevant spatial anatomical information of the reference data set. 
         [0015]    Furthermore, it is possible that different reference data sets of the atrium are generated in order to take into account different heart anatomies, and for each of the different reference data sets transferred treatment-relevant spatial anatomical information is evaluated. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0016]    The invention will be described in further detail with reference to the accompanying drawings. In the figures, 
           [0017]      FIG. 1  shows a schematic view of a left atrium with planned ipsilateral ablation lines, 
           [0018]      FIG. 2  shows a schematic view of a system that can be used to plan ablation treatment areas, 
           [0019]      FIG. 3  shows a schematic view of the left atrium with several pre-planned ipsilateral pulmonary vein ablation lines, 
           [0020]      FIG. 4  shows the schematic left atrium of  FIG. 3  with pre-planned ablation lines transferred from different data sets and an estimated mean ablation planning line, 
           [0021]      FIG. 5  shows a mean and standard deviation of the residual mesh-to-mesh error, a landmark error after registration and initial landmark error per data set and the mean and standard deviation of the initial and residual landmark error, 
           [0022]      FIG. 6  shows a mean and standard deviation for residual mesh-to-mesh error after non-rigid registration of the reference model to the respective template mesh, 
           [0023]      FIG. 7  shows the mean and standard deviation for the distance of pre-planned ablation lines per line segment for left-sided planning lines and right-sided planning lines, 
           [0024]      FIG. 8  shows the orientation and labeling of planning line segments, 
           [0025]      FIG. 9  shows a schematic flow chart including exemplary steps to evaluate treatment areas of several data sets on a reference data set, 
           [0026]      FIG. 10  shows a flow chart comprising steps that are carried out to compute a reference data set, and 
           [0027]      FIG. 11  shows a flow chart comprising steps of applying reference treatment-relevant spatial anatomical information to a new data set. 
       
    
    
     DETAILED DESCRIPTION 
       [0028]    We present a novel method to investigate this problem by comparing pre-planned ablation lines defined by an experienced clinician. Although actual ablation lines may differ from their targets, a comparison among desirable ablation lines does provide insight into how much ablation lines may vary along certain parts of the left atrium. 
         [0029]    We propose a framework for transfer of ablation planning information among different left atria heart models. Initially, a statistical shape model of the left atrium is generated to serve as a common basis. This approach is possible to transfer (spatially organized) annotations from input data to the previously computed common model. By projecting the different annotations onto to the common model, we can compare the different annotations and extract common features. By way of example, this makes it possible to compare different planning structures used in the context for atrial fibrillation ablation. Further on, it can be used as basis for an automatic planning approach. When applied to planning an atrial fibrillation ablation, different anatomical configurations of the left atrium have to be considered. This could be addressed by using a common basis for each configuration, e.g. for pulmonary veins, a common pulmonary vein or supernumerary vein. 
         [0030]    The left atria are segmented from 3-D data sets (CT, C-arm CT, MRI) and modeled as triangle meshes. The correspondence between different left-atrial models is calculated via non-rigid registration of the model. Transfer of planning lines from the mesh-model A to mesh-model B is carried out as follows: First, model A is registered to model B. Afterwards, the planning data from model A is projected to model B. 
         [0031]    The proposed solution can contain the following contributions: First of all, it uses a non-rigid point cloud based registration algorithm, for example coherent point drift, to generate a statistical shape model of the left atrium. 
         [0032]    Secondly, spatially organized information, e.g. pulmonary vein isolation ablation planning lines are transferred and merged across multiple specific anatomies. 
         [0033]    The statistical shape model of the left atrium can be used as a reference model to perform RF ablation planning, e.g. for pulmonary vein isolation (PVI). The second aspect enables us to transfer these annotations onto new unseen mesh models. This functionality can be used for automatic procedure planning. 
         [0034]    This invention proposes a method for left atrium shape modeling using non-rigid point cloud registration. In our approach, we generate the shape model from 3-D magnetic resonance imaging (MRI) volume data sets. We exclusively used data sets of left atria with four pulmonary veins, which reflects the most common anatomic configuration [10]. First, the relevant structure was segmented and represented as triangle mesh. Then we used the Coherent Point Drift (CPD) algorithm [11] to pairwise align meshes via non-rigid point cloud registration. Basically, CPD registration is performed based on a Gaussian Mixture Model (GMM) framework and a regularization of the displacement field. Benefits of the CPD algorithm are the generation of smooth deformation fields while being robust against noise and outliers [11]. 
         [0035]    Left atrium mesh models of ten subjects were extracted from contrast enhanced 3-D MRI volume data sets. The MRI data sets were acquired with a resolution of 256×256×68 voxels. The in-plane pixel spacing was 1.23×1.23 mm and the slice thickness 1.5 mm. The left atrium was segmented from MRI voxel data sets using a semi-automatic segmentation software (syngo InSpace EP, Siemens AG, Forchheim, Germany). The segmentation process is initialized by manually selecting a point inside the left atrium. Based on this seedpoint, the complete left atrium is segmented automatically. The segmentation results are represented as triangle meshes. 
         [0036]    For registration, let us consider the mesh as a point cloud M consisting of N points x i ∈R 3    
         [0000]        M≡m=[x   1   T   , . . . ,x   N   T ] T   ∈R   3N   (1)
 
         [0037]    In a first step, we selected one left atrium mesh model as a reference mesh. The reference mesh was chosen based on visual inspection to clearly express the LA anatomy. The reference mesh m Ref  is then registered to a set of sample meshes {m t } T   t=1 , with T=9, using the CPD algorithm. All meshes have the same anatomical orientation, and are zero centered before applying the registration. 
         [0038]    We used the coherent point drift (CPD) algorithm to register the reference mesh to the set of sample meshes. CPD follows a probabilistic approach by considering the alignment of the two point sets as a probability density estimation problem. The basic idea is to fit the GMM centroids, represented by the points of the reference mesh m Ref , to the sample mesh m t , by maximizing the likelihood. This optimization is performed with the expectation maximization algorithm. During the optimization process, the GMM centroids are forced to move coherently as a group, to ensure preservation of the topological structure of the point set. The displacement function v for the reference mesh is defined as 
         [0000]        {circumflex over (m)}   Ref   =m   Ref   +v ( m   Ref )  (2)
 
         [0000]    with m Ref  as the initial centroid positions, {circumflex over (m)} Ref  and v, respectively are obtained by minimizing the following energy function [12]: 
         [0000]    
       
         
           
             
               
                 
                   
                     E 
                      
                     
                       ( 
                       
                         
                           m 
                           ^ 
                         
                         Ref 
                       
                       ) 
                     
                   
                   = 
                   
                     
                       - 
                       
                         
                           ∑ 
                           
                             n 
                             = 
                             1 
                           
                           N 
                         
                          
                         
                           log 
                            
                           
                             
                               ∑ 
                               
                                 m 
                                 = 
                                 1 
                               
                               M 
                             
                              
                             
                                
                               
                                 
                                   - 
                                   
                                     1 
                                     2 
                                   
                                 
                                  
                                 
                                   
                                      
                                     
                                       
                                         
                                           x 
                                           m 
                                         
                                         - 
                                         
                                           y 
                                           m 
                                         
                                       
                                       σ 
                                     
                                      
                                   
                                   2 
                                 
                               
                             
                           
                         
                       
                     
                     + 
                     
                       
                         λ 
                         2 
                       
                        
                       
                         Φ 
                          
                         
                           ( 
                           v 
                           ) 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   3 
                   ) 
                 
               
             
           
         
       
     
         [0000]    where Φ(v) is a regularization to ensure the displacement field to be smooth. x n  denotes a point of the mesh m t , y m  a point of the transformed mesh {circumflex over (m)} Ref , respectively. N and M refer to the number of points within the respective mesh. The parameter λ determines the trade-off between data fitting and smoothness of the deformation field. We empirically determined a suitable value for this parameter (λ=2.0). 
         [0039]    The reference mesh m Ref  is registered to every sample mesh m t . The transformed mesh {circumflex over (m)} Ref  is labeled y t  for ease of use. The training set is defined as V={m Ref , v 1 , . . . , v T }. We used a Principle Component Analysis (PCA) approach [13] to compute the modes of variation. Applying PCA to the covariance matrix of the centered version of V yields a set of eigenvectors e 1  describing the principle modes of variation in the training data set. The eigenvectors are ordered in descending order based on the value of their corresponding eigenvalue. The P largest eigenvectors are stored in the matrix Φ=[e 1 , . . . , e p ]∈R 3N×P . A linear combination of the P principal modes of variation, with b∈R P as weighting factors, spans a subset of linearized mesh models composed of the given modes of variation: 
         [0000]        m′= ν +Φb   (4)
 
         [0040]    The mean shape  ν  is defined as 
         [0000]    
       
         
           
             
               
                 
                   
                     v 
                     → 
                   
                   = 
                   
                     
                       1 
                       
                         T 
                         + 
                         1 
                       
                     
                      
                     
                       ( 
                       
                         
                           m 
                           Ref 
                         
                         + 
                         
                           
                             ∑ 
                             
                               i 
                               = 
                               1 
                             
                             T 
                           
                            
                           
                             v 
                             i 
                           
                         
                       
                       ) 
                     
                   
                 
               
               
                 
                   ( 
                   5 
                   ) 
                 
               
             
           
         
       
     
         [0041]    For quantitative evaluation of the proposed framework, we used ten clinical data sets with manually annotated pulmonary veins (PV) ostia. These landmarks are labeled RSPV (Right Superior Pulmonary Veins), RIPV (Right Inferior Pulmonary Veins), LSPV (Left Superior Pulmonary Veins), and LIPV (Left Inferior Pulmonary Veins). The quality of the registration is measured based on residual landmark distances and mesh-to-mesh distance. The residual landmark error is defined as the Euclidean distance of the center of corresponding PV ostia and measured after non-rigid CPD registration. 
         [0042]    Quantitative registration results are shown in  FIG. 5  and discussed in more detail below.  FIG. 5( a )  shows the average mesh-to-mesh distance, as well as the initial and residual landmark offset per data set. The initial and residual offsets per landmark are shown in  FIG. 5( b ) . The average mesh-to-mesh distance is between 2.5 and 5.1 mm, the average landmark offset is between 2.9 and 13.9 mm. The highest mesh-to-mesh distance occurs at the end of the pulmonary veins. In this case especially on the right inferior PV. 
         [0043]    We described a method for left atrium shape modeling using non-rigid point cloud registration. The overall performance of the mesh registration shows a mean mesh-to-mesh error of 3.4 mm over all data sets. The coherent point drift algorithm was capable of dealing with high variations in anatomy. The highest residual mesh-to-mesh distance results from different extents of the pulmonary veins. The average landmark offset was 8.5 mm. Landmarks on the right side of the left atrium, namely RSPV and RIPV, show a lower residual error compared to left sided landmarks LSPV and LIPV. This might be due to the additional pouch on the left side of the left atrium, the left atrial appendage, which is anterior to the PV ostia. 
         [0044]    The mesh models of the left atrium also contained a large part of the connected pulmonary veins. Removing or trimming these extensions might improve the accuracy, since these structures show a high variation in shape and size. For the modeling of the atrium, short pulmonary vein ostia would be sufficient. This work is a first step towards our goal of automatic planning of ablation regions for atrial fibrillation procedures. Planning structures could be transferred to augmented fluoroscopy systems used to guide the procedure and overlaid to the X-ray images. 
         [0045]    It has been found that the highest mesh-to-mesh distance occurs at the end of the pulmonary veins. 
         [0046]    In the following, the comparison of pre-planned ablation lines for the treatment of atrial fibrillation using a common reference model is discussed in more detail. 
         [0047]    We used left atrial (LA) surface models of seven different patients to build up a common reference shape. The mesh models are represented as a triangulated mesh structure. The models cover the left atrium as well as a certain extent of the attached pulmonary veins. The length of the pulmonary veins varies from data set to data set For consistency, a pre-processing step has been applied to all LA mesh models. The attached pulmonary veins have been removed from the mesh about 2 cm distal to the LA antrum. 
         [0048]    The common reference shape is generated via pairwise non-rigid registration of the LA models. For registration, the mesh is seen as point cloud M consisting of N points x i ∈R 3    
         [0000]        M≡m=[x   1   T   , . . . , x   N   T ] T   ∈R   3N .  (6)
 
         [0049]    All meshes have the same anatomical orientation, and are zero centered before applying the registration. 
         [0050]    We used the Coherent Point Drift (CPD) algorithm [12] to perform the non-rigid point cloud registration. The main benefit of the CPD algorithm is the robustness against nose and outliers while generating smooth deformation fields. A pivot mesh m pivot  is registered to template meshes m t     i   . The variables t i  are used to refer to each of the T=7 template meshes. 
         [0051]    Below, we describe how to derive a common reference model. Then we establish mean ablation lines. Afterwards, we evaluate how individual pre-planned ablation lines vary around their mean as we move around the left and right ipsilateral PVs. 
         [0052]    In the first step, we selected one left atrium mesh model as pivot element, 
         [0000]      m pivot =m t     1     (7)
 
         [0053]    In the next step, we registered the pivot mesh to the remaining template meshes {m t     2   , . . . , m t     T   }. We used the coherent point drift algorithm to compute the non-rigid transformation between the pivot mesh and the other left atrial mesh models. CPD follows a probabilistic approach by considering the alignment of the two point sets as a probability density estimation problem. The basic idea is to fit the Gaussian mixture model (GMM) centroids, represented by the points of the pivot mesh m pivot , to the template mesh m t     i   , by maximizing the likelihood. This optimization is performed with the expectation maximization algorithm. During the optimization the GMM centroids are forced to move coherently as a group, to ensure preservation of the topological structure of the point set. 
         [0054]    For each template mesh m t     i   , the estimated deformation field μ t     i   ∈R 3N  is calculated by registration of the pivot mesh m pivot  to m t     i   . The transformed pivot mesh can then be described as 
         [0000]        {tilde over (m)}   t     i     =m   pivot +μ t     i     (8)
 
         [0055]    The common reference model or reference data set is defined as the mean shape given by 
         [0000]    
       
         
           
             
               
                 
                   
                     m 
                     ref 
                   
                   = 
                   
                     
                       1 
                       
                         T 
                         - 
                         1 
                       
                     
                      
                     
                       ( 
                       
                         
                           ∑ 
                           
                             i 
                             = 
                             2 
                           
                           T 
                         
                          
                         
                           
                             m 
                             ~ 
                           
                           
                             t 
                             i 
                           
                         
                       
                       ) 
                     
                   
                 
               
               
                 
                   ( 
                   9 
                   ) 
                 
               
             
           
         
       
     
         [0056]    Putting (8) and (9) together, it is easy to see that the reference mesh is comprised of the selected pivot mesh, and a mean deformation field. In other words, 
         [0000]    
       
         
           
             
               
                 
                   
                     m 
                     ref 
                   
                   = 
                   
                     
                       m 
                       pivot 
                     
                     + 
                     
                       
                         1 
                         
                           T 
                           - 
                           1 
                         
                       
                        
                       
                         
                           ( 
                           
                             
                               ∑ 
                               
                                 i 
                                 = 
                                 2 
                               
                               T 
                             
                              
                             
                               u 
                               
                                 t 
                                 i 
                               
                             
                           
                           ) 
                         
                         . 
                       
                     
                   
                 
               
               
                 
                   ( 
                   10 
                   ) 
                 
               
             
           
         
       
     
         [0057]    The selection of a proper pivot mesh is important, because it determines the basic shape of the resulting mean mesh. This is why we carefully selected the pivot mesh, making sure that all relevant landmarks, namely the pulmonary veins and left atrium appendage, were clearly expressed. 
         [0058]    Planning lines l are represented as a set of points l={x 1 , . . . x p } with x∈R 3 . Each template mesh m t     i    has two planning lines l t     i     ,R  and l t     i     ,L  attached. They represent desirable ablation lines for right and left sided ipsilateral pulmonary veins, respectively. The planning lines are a subset of the corresponding template mesh, i.e. {l t     i     ,R ,l t     i     ,L }⊂m t     i   . 
         [0059]    To transfer the planning lines from a template mesh onto the reference model, m ref  is registered to m t     i    using CPD. After registration, the two mesh models are optimally aligned based on the optimization criterion stated in (5). The planning lines {l t     i     ,R ,l t     i     ,L } are now projected onto the transformed reference model {tilde over (m)} t     i    with a nearest neighbor approach. By applying the inverse deformation field, the lines can be mapped onto the reference model m ref . There, they are labeled {circumflex over (l)} t     i     ,R  and {circumflex over (l)} t     i     ,L  for right and left sided planning, respectively. 
         [0060]    The mean planning lines are derived from the set of re-mapped ablation lines defined as 
         [0000]      L={{circumflex over (l)} t     i     ,R ,{circumflex over (l)} t     i     ,L } T   i= .  (11)
 
         [0061]    Initially, each planning line consists of an arbitrary number of points. For consistency, each line l was interpolated with a cubic spline, and equidistantly sampled with a fixed number of sample points P L . To investigate the spread and distribution among the pre-planned ablation lines, a common orientation and labeling was enforced. Each planning line represents a closed loop encircling the LA. The point closest to the top is defined as the starting point, and the remaining points are traversed in anterior direction. 
         [0062]    After correct alignment of L, mean reference lines l ref,R  and l ref,L  are generated by averaging corresponding points along the interpolated lines. 
         [0063]    In  FIG. 1 , a schematic view of a left atrium  10  is shown including the pulmonary veins  11 ,  12 ,  13 ,  14 . An example of ipsilateral ablation planning lines  15  and  16  is shown. 
         [0064]    In  FIG. 3 , several pre-planned ablation lines transferred from different meshes are shown. In  FIG. 3  three different ablation lines  15   a - 15   c  and  16   a - 16   c  are shown. 
         [0065]    In  FIG. 4  in addition to the ablation lines already shown in  FIG. 3 , estimated mean ablation lines  17  and  18  are shown, which were calculated as discussed above. 
         [0066]    We evaluated our approach on T=7 LA meshes with attached planning lines. The pre-planned ablation lines were placed by an experienced clinician. The registration accuracy of the non-rigid registration of the reference model to the template meshes in terms of residual average mesh-to-mesh error is shown in  FIG. 6 . The mean residual mesh-to-mesh error calculated over all meshes was 2.0 mm. Data set  1  was selected as pivot mesh for the reference mesh generation. As can be seen in (8), the reference mesh is comprised of the pivot mesh with an additive deformation term. Hence, non-rigid registration of the reference mesh onto data set  1  is possible with a very low residual error, as illustrated in the first column of  FIG. 6 . 
         [0067]    The deviation of the re-mapped planning lines from their respective mean is evaluated as follows. For each point of the reference planning lines l ref,R  and l ref,L , the distance to the re-mapped planning lines {circumflex over (l)} t     i     ,R , {circumflex over (l)} t     i     ,L  with i=1 . . . T is calculated. 
         [0068]    For further, more anatomically oriented evaluations, the planning lines were divided into eight equally spaced segments, as depicted in  FIG. 8 , in which a clock schema is used to divide the substantial circular planning lines into different segments. Based on this convention, we present the mean deviation of the individual planning lines per line segment in  FIG. 7 . Line Segment  1  and  2  are located on the roof of the LA.  5  and  6  are below the inferior PV, and  7  and  8  are on the posterior side of the LA. Note that there is an additional pouch on the left side of the LA, the left atrial appendage (LAA). The LAA is located anterior the left-sided PVs, separated from them by an arrow ridge. The LAA runs on the left side along segments  2  to  5 . 
         [0069]    The average distance of the re-mapped ablation lines to the reference ablation line was 2.9±1.9 mm and 1.8±1.5 mm for right and left sided planning lines, respectively. The maximum distance of a single planning point to the reference planning line was 8.5 mm and 7.6 mm for right and left side, respectively. 
         [0070]    According to our analysis, the average deviation over all individual planning lines that were part of our data set was 2.4±1.8 mm. This results in a region for pre-planned ablation lines that is approximately twice this width. The smallest deviation was encountered in line segments  2  and  3 . We believe that this is due to the presence of the ridge between left atrial appendage and left PVs, that leaves limited space for setting up pre-planned ablation lines. 
         [0071]    Based on feedback from physicians in this field, our findings appear plausible. In a nutshell, by comparing pre-planned ablation lines placed by an experienced clinician on LA models of actual patients, we found that one seems to have some 5 mm of “wiggle room” despite the presence of prominent anatomical structures. After analyzing the inter-patient variance of manually placed planning lines and learning a general planning pattern, the design of an algorithm for fully automatic pre-planning of ablation lines will be the next step. 
         [0072]    In  FIG. 9 , we summarize the method for evaluating treatment-relevant spatial anatomical information on a reference data set. An example for the treatment-relevant spatial anatomical information are the ablation planning lines in the left atrium. 
         [0073]    The method starts in step S 90 . In a step S 91 , a reference three-dimensional data set of the atrium is provided. In step S 92 , further three-dimensional data sets with treatment-relevant spatial anatomical information are provided, e.g. a first data set, a second data set and eventually further data sets. In step S 93 , the reference data set is registered to each of the data sets provided in step S 92 . 
         [0074]    In step S 94 , the treatment-relevant spatial anatomical information of each of the data sets provided in step S 92  is transferred on the reference data set so that transferred treatment-relevant spatial anatomical information is attached to the reference data set. In step S 95 , the transferred treatment-relevant spatial anatomical information attached to the reference data set is compared and evaluated. This means that by way of example, a mean and a standard deviation of the planning lines transferred to the reference data set is determined. In the next step S 96 , a reference treatment-relevant spatial anatomical information is determined, e.g. reference treatment planning lines. The method ends in step S 97 . 
         [0075]    In  FIG. 10 , a flow chart is shown comprising the steps which can be carried out to compute a reference three-dimensional data set. The method starts in step S 100 . In another step S 101 , the different three-dimensional anatomical data sets are provided. In step S 102 , a pivot data set is selected. This pivot data set may be selected based on visual inspection of the different data sets and serves as basis for the generated reference data set. In step S 103 , the pivot data set is registered to every of the other three-dimensional data sets. In step S 104 , it is possible to determine a statistical shape model, e.g. a mean shape as explained in more detail above. In step S 105 , a reference model or reference data set of the left atrium can be determined as mentioned above in equation 10. As discussed above, the reference model can be comprises of the selected pivot data set and a mean deformation field. 
         [0076]    In  FIG. 11 , a flow chart is shown comprising the steps which can be used to plan an interventional treatment in an atrium. The method starts in step  110 . In the method shown in  FIG. 11 , the transferred treatment-relevant spatial anatomical information is applied to a new data set. In step S 111 , the reference three-dimensional data set is provided. Furthermore, reference treatment-relevant spatial anatomical information on the reference data set is provided, e.g. as reference treatment planning lines (step S 112 ). In step S 113 , the new three-dimensional image data set is provided for which the planning lines should be determined. In step S 114 , the left atrium in the new data set is segmented and in step S 115 , the reference data set is registered to the new data set. In step S 116 , the reference treatment-relevant spatial anatomical information is transferred to the new data set to determine planning lines for the new data set. The method ends in step S 117 . 
         [0077]    In  FIG. 2 , a system is schematically shown with which ablation lines can be evaluated and with which a planning of ablation lines is possible. The three-dimensional images can be generated with an imaging system  21 , which can be an MR or CT scanner. The different three-dimensional data sets generated by the imaging system are stored in data base  22 . The system may not necessarily contain the imaging system  21 . It is also possible that the system only contains the data base  22  with the prestored three-dimensional data sets of the heart or of the atrium. With a processing module  23 , the different steps discussed above and the calculations discussed in detail above can be carried out. A human machine interface, HMI,  24  is provided with which a user can interact with the system, e.g. can draw ablation lines, select data sets, etc. The information provided by the system may be displayed on a display  25 . 
         [0078]    The different components shown in  FIG. 2  need not be provided in the form shown in  FIG. 2 . The functional modules shown in  FIG. 2  may be incorporated by another combination of modules. Furthermore, the different modules may be incorporated by hardware or software. 
       REFERENCES 
       [0000]    
       
         [1] Heimann T, Meinzer H P, “Statistical shape models for 3D medial image segmentation”,  A review. Med Image Anal.  2009; 13(4): 543-563 
         [2] Keustermans J, De Buck S, Heidbuechel H, Suetens P, “Automated planning of ablation targets in atrial fibrillation treatment”,  Proc SPIE.  2011; 7962:796207 
         [3] Stijn De Buck, Frederik Maes, Joris Ector, Jan Bogaert, Steven Dymarkowski, Hein Heidbuechel, and Paul Suetens, “An augmented reality system for patient-specified guidance of cardiac catheter ablation procedures”.  IEEE Transactions on Medical Imaging, vol.  24, no. 11, pp. 1512 - 1524, November 2005 
         [4] Martin Koch, Arne Langenkamp, Atilla Kiraly, Alexander Brost, Norbert Strobel, and Joachim Hornegger, “Navigation system with contact force assessment to guide pulmonary vein isolation procedures,” in 23 rd    Conference of the Society for Medical Innovation and Technology  ( SMIT ), Tel Aviv, Sep. 13-16, 2011 
         [5] Felix Bourier, Dejan Vukajlovic, Alexander Brost, Joachim Hornegger, Norbert Strobel and Klaus Kurzidim, “Pulmonary vein isolation supported by MRI-derived 3D-augmented biplane fluoroscopy: A feasibility study and a quantitative analysis of the accuracy of the technique,”  Journal of Cardiovascular Electrophysiology,  2012 [electronic publication ahead of print] 
         [6] M. Hastenteufel, I. Wolf, C. Christoph, S. Yang, T. Boettger, M. Vetter, and H P Meinzer, “A novel method for planning and visualization of ablation lines for atrial fibrillation treatment,” in  Computers in Cardiology . IEEE, 2004, pp. 13-16 
         [7] Alexander Brost, Felix Bourier, Andreas Kleinoeder, Jens Raab, Martin Koch, Marc Stamminger, Joachim Hornegger, Norbert Strobel, and Klaus Kurzidim, “AFiT—Atrial Fibrillation Ablation Planning Tool,” in  Vision, Modeling and Visualization  ( VMV ), 2011, pp. 223-230 
         [8] Roderick Tung, Eric Buch, and Kalyanam Shivkumar, “Catheter ablation of atrial fibrillation,”  Circulation , vol. 126, no. 2, pp. 223-229, 2012 
         [9] Thomas Arentz, Reinhold Weber, Gerd Brkle, Claudia Herrera, Thomas Blum, Jochem Stockinger, Jan Minners, Franz Josef Neumann, and Dietrich Kalusche, “Small or large isolation areas around the pulmonary veins for the treatment of atrial fibrillation?”,  Circulation , vol. 115, no. 24, pp. 3057-3063, 2007 
         [10] Kautzner J, Micochova H, Peichl P, “Anatomy of the Left Atrium and Pulmonary Veins—Lessons Learned from Novel Imaging Techniques”.  Eur Cardiol.  2006; 2(1): 89-90 
         [11] Myronenko A, Song X, “Point set registration: Coherent Point Drift”.  IEEE Trans Pattern Anal Mach Intell.  2010; 32(12):2262-2275 
         [12] Myronenko A, Song X, Carreira-Perpinan M, “Non-rigid point set registration: Coherent Point Drift”.  Adv Neural Inf Process Syst.  2007; 19:1009-1016 
         [13] Jolliffe I, “Principal component analysis”,  Wiley Online Library;  2005