A Bilinear Model for Temporally Coherent Respiratory Motion
|
|
- Eugene George
- 5 years ago
- Views:
Transcription
1 A Bilinear Model for Temporally Coherent Respiratory Motion Frank Preiswerk (B) and Philippe C. Cattin Medical Image Analysis Center, University of Basel, Basel, Switzerland Abstract. We propose a bilinear model of respiratory organ motion. The advantages of classical statistical shape modelling are combined with a preconditioned trajectory basis for separately modelling the shape and motion components of the data. The separation of a linear basis into bilinear form leads to a more compact representation of the underlying physical process and the resulting model respects the temporal regularity within the training data, which is an important property for modelling quasi-periodic data. Bilinear modelling is combined with a Bayesian reconstruction algorithm for sparse data under observation noise. By applying the model to liver motion data, we show that our bilinear formulation of respiratory motion is significantly more parsimonious and can even outperform linear PCA-based models. Keywords: Respiratory motion Bilinear model Liver motion 1 Introduction Modelling respiratory organ motion is an active field of research. Having an accurate model of respiratory motion is desirable for many practical applications including motion segmentation, registration, tracking and reconstruction as well as tumour tracking in a clinical scenario. A comprehensive overview of the field is available in [1]. Such models can serve different purposes, two of which are worth mentioning here. First, respiratory motion models are used to estimate the organ position for a given point in time, e.g. based on sparse measurements. Second, they are used to predict the future position based on current and/or past values. Respiratory motion is approximately repetitive, thus it is reasonable to take the repetitive nature into account when designing a model. While respiratory motion is very similar from cycle-to-cycle, there are variations that must be taken into account. On the other hand, we wish the model not to be able to describe arbitrary motion patterns. This trade-off between generalisation and specificity is typically encountered in model selection. Recently, statistical motion models which are learned from data, typically by employing Principal Component Analysis (PCA) techniques, have been proposed in various studies [2 7]. In [4], a kernel-pca is learned to model the relationship between fiducial movement and lung surface motion. In [7] PCA, is computed c Springer International Publishing Switzerland 2014 H. Yoshida et al. (Eds.): ABDI 2014, LNCS 8676, pp , DOI: /
2 222 F. Preiswerk and P.C. Cattin on dense motion fields extracted from 4DMRI of the liver through non-rigid registration. In [2], the coefficients of a PCA model are further restricted by a single governing parameter as prior which models the distribution of the observed data s PCA coefficients. Although all mentioned approaches provide a compact, low-dimensional parameterisation of valid shapes within a breathing cycle, they do not take into account the temporal regularity of respiratory motion. In [8], a bilinear factorisation of the model coefficients of each sample is proposed to compute a statistical model for reconstruction of the heart at discrete points in the cardiac cycle, inspired by the work on separating style and content [9]. Recently, a bilinear spatiotemporal model was proposed for facial animation [10]. This model factorises the basis of a dataset into separate shape and trajectory bases, as compared to the bilinear factorisation of the model coefficients in [8, 9]. In this paper, we propose a general approach to model respiratory organ motion based on a bilinear model. We augment the bilinear model with a Bayesian algorithm for reconstruction from sparse and noisy data, leading to a variety of interesting possible applications. 2 Bilinear Motion Model Let the vector x = (x 1,y 1,...,x n,y n ) T R p be an instance of a 2d shape (w.l.o.g for higher dimensions) and the matrix X =[x 1, x 2,...,x t ], (1) a sequence of t shapes. In a linear model of shape, each instance is represented by a linear combination of base shapes x j = i b i c i,j (2) and the complete sequence can thus be represented as X = B s C s, (3) where B s R p ks consists of k s basis vectors b i R p that span the model space and C s R ks t is a matrix of t coefficient vectors that define linear combinations of shape basis vectors. Conversely, a linear model of motion can be built where the rows (not the columns) of X are modeled as a linear combination of motion trajectories, x i = c i,j b j, (4) j and the data matrix is represented as X = C m B T m, (5) where B m =[b 1, b 2,...,b km ] R t km and C m R p km are k m motion basis and p coefficient vectors, respectively [10]. In practice, the dimensionality of a linear model is often reduced using some type of dimensionality reduction approach,
3 A Bilinear Model for Temporally Coherent Respiratory Motion 223 e.g. principal component analysis (PCA). However, in a linear representation, the temporal regularity of trajectories is ignored. As an illustrative example, consider an ordered sequence of shape vectors [x 1, x 2,...,x t ]. Any permutation [x p(1), x p(2),...,x p(t) ] of the columns in X only results in a permutation of the coefficients in C s but does not affect the shape basis B s. The same also applies forshufflingrowsofx, affecting only rows in C m. In this sense, a linear model is an overparameterisation because it fails to exploit the underlying spatiotemporal structure. For modelling respiratory motion, we would like to respect temporal regularities in the data. A similar problem has been studied in [10], where it is shown that a factorisation X = B s CB T m (6) exists that links together shape and motion bases in a single model. In this bilinear model, the matrix C R ks km defines weights for the outer products of the i-th shape vector b i s and the j-th trajectory vector b j m. Intuitively, it describes how the points of the shape modes b i s vary over time. Figure 1 shows how the coefficients c ij in C are related to the outer product of shape and trajectory basis vectors b i s and b j m. p k s b i s ks k m t k m c ij b j m B s : shape basis C : coefficient matrix B T m : trajectory basis Fig. 1. Illustration of bilinear multiplication. The coefficients c ij define a weighting of each outer product between shape basis vectors b i s and trajectory basis vectors b j m. In general, both bases and the coefficient matrix must be estimated from data using tensor decomposition methods. However, in [10] the problem is considerably simplified by using a conditioned motion basis. It has been shown that the discrete cosine transform (DCT) basis converges to the optimal PCA basis if the data is generated from a stationary first-order Markov process [11]. Moreover, for respiratory motion, we are interested in low-frequency cyclic patterns, which makes the DCT basis a suitable choice for compression. The use of an analytical basis like DCT is also useful for adapting the bilinear model to varying cycle lengths. Since the shape and motion bases factor out separately and as a result, only B m depends on t, the motion basis can be adjusted to appropriate size by changing the length of the basis vectors b j m while leaving the shape basis unchanged. This is not possible with a linear model because shape and motion information are intermingled.
4 224 F. Preiswerk and P.C. Cattin Given a conditioned motion basis B m and a series of respiratory cycle data X 1,...,X N, we can solve for a shape basis that minimises the squared reconstruction error. It is given by the row-space computed through singular value decomposition (SVD) of the matrix Π =[ˆX T 1, ˆX T 2,..., ˆX T N ] T, (7) where ˆX i = X i B m B + m denotes the reconstruction of X from its trajectory projection and + denotes the Moore-Penrose pseudoinverse. This way, the shape basis B s can be computed directly from the data in closed form. It defines a space where variations along the individual dimensions are uncorrelated. In order to simplify the calculations in the next section, we assume that coefficients are scaled to unit variance in each of the two subspaces, i.e. C = Σs 1 B T s XB m Σm 1. 3 Reconstruction from Sparse Data Applications for the bilinear motion model are manifold. A natural application is to use it as a prior for motion segmentation and registration. Here, we focus on another scenario, which is the reconstruction from only a small set of observations, e.g. when measurements from implanted electromagnetic beacons or 2d projections are given and the goal is to reconstruct the entire shape over time. In order to estimate entire breathing cycles from partial observations, we extend the scheme presented in [12] to the bilinear case. Our observation model is as follows, R = L s XL T m, (8) where L s : R p R lp,l p <pacts on the columns (i.e. features) and L m : R t R lt,l t <tacts on the rows (i.e. time frames) of X. Although these matrices can be any linear transformation, in our context L s typically deletes rows (features) over the entire sequence and L m deletes all features in a particular frame. In general, both mappings are not injective, thus there is no unique solution of Eq. 8. The least-squares solution is found by optimising the error function E(C) = Q s CQ + m R F, (9) where Q s = L s Σ s B s, Q m = L m Σ m B T m and F denotes the Frobenius norm. This cost function is minimal for Ĉ = Q + s RQ m. (10) However, the result in Eq. 10 only applies for noiseless measurements. Assuming that measurements are subject to uncorrelated Gaussian noise of variance σn 2, the likelihood of observing R is given by p(r C m )=ν exp( 1 2σN 2 Q s C m R F ), (11)
5 A Bilinear Model for Temporally Coherent Respiratory Motion 225 where we have defined C m = CQ + m and with normalisation constant ν = ((2π) l 2 σn ) 1. Using Bayes rule to compute the posterior probability and maximising it leads to the following minimisation problem, E(C m )= Q s C m R F + η C m F. (12) The regularisation factor η = σn 2 allows to find a trade-off between matching quality and prior probability. It can be shown [12] that the global optimum is obtained via SVD of Q s = U s W s Vs T according to w s,i Ĉ η = V s diag( ws,i 2 + η )UT s R, (13) where w s,i are the diagonal elements of W s. Note that all matrices in Eq. 13 except R are constant and can thus be precomputed. 3.1 Model Size and Compression Any bilinear model can be rewritten in linear form by computing the Kronecker product, B = B m B s R p t, (14) i.e. taking all possible products between the elements of the shape and trajectory bases. Note that this result does not imply that linear and bilinear models are equivalent because the opposite operation, i.e. obtaining a bilinear factorisation from a linear basis, is not always possible. Retaining only the k s basis vectors in B s with highest corresponding singular values leads to optimal dimensionality reduction in a least-squares sense. Since respiratory motion can be encoded primarily by lower-frequency components, the model can be further compressed by only retaining the k m lowest-frequency basis vectors in B m. For example, setting k s = k t = 5 in a model of p = 100 coordinates and 10 time steps leads to a shape basis of size and a motion basis of size 10 5, totalling in 550 values. The corresponding linear model according to Eq. 14 is of size or more than 45 times the number of entries. Formally, the size of the bilinear model grows linearly in the number of model points O(p k s + t k m ) while the size of the corresponding linear model grows according to O(p t k s k m ), which reflects the fact that the bilinear model exploits the structure of the underlying data more efficiently. 4 Experiments We applied the bilinear model to a sequence of sagittal 2d MRI images of the liver of a male subject under free breathing. A total of 120 respiratory cycles was acquired and split into a training set of 40 cycles and a test set of n =80cycles. The images were pre-processed as follows. First, they were non-rigidly registered using NiftyReg [13], an open-source medical image registration suite, by deforming a reference exhalation image to all other images. Due to the discontinuities
6 226 F. Preiswerk and P.C. Cattin MAE [mm] Bilinear model of cycles Linear model of cycles Linear model of states MAE [mm] Bilinear model of cycles Linear model of cycles Linear model of states η (a) Uncompressed models η (b) Compressed models Fig. 2. Mean absolute error (MAE) of reconstruction from full testset over different values of the regularisation factor η. In (a), the complete bases were used. In (b), k s = 5 shape and k m = 5 motion basis vectors were used in the bilinear case and 25 basis vectors were used for the linear models. between the sliding organ and the abdominal wall, the liver was manually masked in the reference image. Additionally, the data was normalised wrt. the number of frames in each cycle, resulting in exactly t = 10 states per cycle. The displacement fields were then equidistandly sampled at n = 85 locations per frame and a bilinear model was built from these vertices as described above. 4.1 Model Compression To compare the proposed bilinear model of respiratory cycles, we computed a linear PCA model of cycles (where each of the n test cycles is concatenated into a column vector of size p t) as well as a linear PCA model of states (where the samples consist of t n column vectors of length p). For the PCA models we use the linear version of the same reconstruction algorithm described above. We reconstructed the full training set without additional noise to see how well the models can describe unseen data. As can be seen in Fig. 2(a), the uncompressed bilinear model has high generalisation capabilities, especially if no regularisation is applied (see η = 0). When the models are compressed, the bilinear model still outperforms the linear model of cycles for all values of η (Fig. 2(a)). It is not surprising that the linear model of states performs well in both experiments because it fits every observed state individually rather than fitting entire cycles. 4.2 Missing Values In order to elaborate on the performance of the models in a setting where data is missing and additionally subject to measurement noise, we performed reconstructions on various percentages of randomly removed data. Figure 3(a) shows the error plot for randomly removing an increasing number of model vertices
7 A Bilinear Model for Temporally Coherent Respiratory Motion Bilinear model of cycles Linear model of cycles Linear model of states 2 Bilinear model of cycles Linear model of cycles MAE [mm] 2 MAE [mm] Removed model points [%] (a) Removed frames [%] (b) Fig. 3. Mean absolute error (MAE) of reconstruction with Gaussian random noise of σ N =1.0 and fixed regularisation factor η =1.0. (a) Model vertices randomly removed over the entire test sequence. (b) Complete time frames randomly removed (linear model of states is not able to reconstruct missing frames). over the entire test sequence. The bilinear model performs better or equal to its linear counterpart, while the linear model of states clearly performs worse in most cases. Lastly, we randomly removed an increasing number of entire frames from the sequence. Figure 3(b) shows the results. Here, the linear model shows a superior error curve that increases exponentially, while the error of the bilinear model increases linearly. This experiment does not apply for the linear model of states because it cannot produce meaningful reconstructions in case of missing frames. 5 Conclusion We presented a novel model for modelling respiratory motion that respects the temporal regularity of the underlying data. Shape and motion information is separated into individual bases and a coefficient matrix defines weights of the basis vectors outer products. This separation leads to significantly more compact models compared to conventional linear models. Furthermore, the model is more flexible because the motion basis can be independently exchanged while leaving the remaining components of the model unchanged in order to generate cycle instances of arbitrary length, although we have restricted our experiments to cycles of the same length for simplicity. We described a Bayesian algorithm for reconstruction from noisy and sparse data in the bilinear case and we could show that the model performs better in most experiments compared to its less compact linear counterparts. All components of the reconstruction algorithm can be precomputed, therefore the algorithm is computationally efficient. The model can potentially be used in registration algorithms as well as for motion compensation in tumour therapy.
8 228 F. Preiswerk and P.C. Cattin References 1. McClelland, J., Hawkes, D.J., Schaeffter, T., King, A.P.: Respiratory motion models: a review. Med. Image Anal. 17, (2013) 2. King, A.P., Blackall, J.M., Penney, G.P., Hawkes, D.J.: Tracking liver motion using 3-D ultrasound and a surface based statistical shape model. In: Proceedings of IEEE Workshop on Mathematical Methods in Biomedical Image Analysis, pp (2001) 3. Zhang, T., Lu, W., Olivera, G.H., Keller, H., et al.: Breathing-synchronized delivery: a potential four-dimensional tomotherapy treatment technique. Int. J. Radiat. Oncol. Biol. Phys. 68, (2007) 4. He, T., Xue, Z., Xie, W., Wong, S.T.C.: Online 4-D CT estimation for patientspecific respiratory motion based on real-time breathing signals. In: Jiang, T., Navab, N., Pluim, J.P.W., Viergever, M.A. (eds.) MICCAI 2010, Part III. LNCS, vol. 6363, pp Springer, Heidelberg (2010) 5. Li, R., Lewis, J., Jia, X., Zhao, T., et al.: PCA-based lung motion model. Technical report (2010) 6. Schneider, M.: Model-based respiratory motion compensation for image-guided cardiac interventions. In: Informatiktage. LNI, vol. S-9, pp GI (2010) 7. Preiswerk, F., Arnold, P., Fasel, B., Cattin, P.C.: A bayesian framework for estimating respiratory liver motion from sparse measurements. In: Yoshida, H., Sakas, G., Linguraru, M.G. (eds.) Abdominal Imaging. LNCS, vol. 7029, pp Springer, Heidelberg (2012) 8. Hoogendoorn, C., Sukno, F.M., Ordás, S., Frangi, A.F.: Bilinear models for spatiotemporal point distribution analysis. Int. J. Comput. Vis. 85, (2009) 9. Tenenbaum, J.B., Freeman, W.T.: Separating style and content with bilinear models. Neural Comput. 12, (2000) 10. Akhter, I., Simon, T., Khan, S., Matthews, I., Sheikh, Y.: Bilinear spatiotemporal basis models. ACM Trans. Graph. 31, 17:1 17:12 (2012) 11. Rao, K.R., Yip, P., Rao, K.R.: Discrete Cosine Transform: Algorithms, Advantages, Applications, vol Academic Press, Boston (1990) 12. Blanz, V., Vetter, T.: Reconstructing the complete 3D shape of faces from partial information. Informationstechnik und Technische Informatik 44, (2002) 13. Modat, M., Ridgway, G.R., Taylor, Z.A., Lehmann, M., et al.: Fast free-form deformation using graphics processing units. Comput. Meth. Prog. Biomed. 98, (2010)
A Bayesian Framework for Estimating Respiratory Liver Motion from Sparse Measurements
A Bayesian Framework for Estimating Respiratory Liver Motion from parse Measurements Frank Preiswerk, Patrik Arnold, Beat Fasel, and Philippe C. Cattin Medical mage Analysis Center, University of Basel,
More informationRobust Tumour Tracking From 2D Imaging Using a Population-Based Statistical Motion Model
Robust Tumour Tracking From 2D Imaging Using a Population-Based Statistical Motion Model Frank Preiswerk, Patrik Arnold, Beat Fasel and Philippe C. Cattin Medical Image Analysis Center University of Basel,
More informationImage Guidance and Beam Level Imaging in Digital Linacs
Image Guidance and Beam Level Imaging in Digital Linacs Ruijiang Li, Ph.D. Department of Radiation Oncology Stanford University School of Medicine 2014 AAPM Therapy Educational Course Disclosure Research
More informationQuantication and analysis of respiratory motion from 4D MRI
Journal of Physics: Conference Series OPEN ACCESS Quantication and analysis of respiratory motion from 4D MRI To cite this article: Ashrani Aizzuddin Abd Rahni et al 2014 J. Phys.: Conf. Ser. 546 012001
More informationFully Automatic Multi-organ Segmentation based on Multi-boost Learning and Statistical Shape Model Search
Fully Automatic Multi-organ Segmentation based on Multi-boost Learning and Statistical Shape Model Search Baochun He, Cheng Huang, Fucang Jia Shenzhen Institutes of Advanced Technology, Chinese Academy
More informationMethodological progress in image registration for ventilation estimation, segmentation propagation and multi-modal fusion
Methodological progress in image registration for ventilation estimation, segmentation propagation and multi-modal fusion Mattias P. Heinrich Julia A. Schnabel, Mark Jenkinson, Sir Michael Brady 2 Clinical
More informationFacial Expression Recognition using Principal Component Analysis with Singular Value Decomposition
ISSN: 2321-7782 (Online) Volume 1, Issue 6, November 2013 International Journal of Advance Research in Computer Science and Management Studies Research Paper Available online at: www.ijarcsms.com Facial
More informationADAPTIVE LOW RANK AND SPARSE DECOMPOSITION OF VIDEO USING COMPRESSIVE SENSING
ADAPTIVE LOW RANK AND SPARSE DECOMPOSITION OF VIDEO USING COMPRESSIVE SENSING Fei Yang 1 Hong Jiang 2 Zuowei Shen 3 Wei Deng 4 Dimitris Metaxas 1 1 Rutgers University 2 Bell Labs 3 National University
More informationA Non-Linear Image Registration Scheme for Real-Time Liver Ultrasound Tracking using Normalized Gradient Fields
A Non-Linear Image Registration Scheme for Real-Time Liver Ultrasound Tracking using Normalized Gradient Fields Lars König, Till Kipshagen and Jan Rühaak Fraunhofer MEVIS Project Group Image Registration,
More informationRegistration of Expressions Data using a 3D Morphable Model
Registration of Expressions Data using a 3D Morphable Model Curzio Basso, Pascal Paysan, Thomas Vetter Computer Science Department, University of Basel {curzio.basso,pascal.paysan,thomas.vetter}@unibas.ch
More informationClassification of Subject Motion for Improved Reconstruction of Dynamic Magnetic Resonance Imaging
1 CS 9 Final Project Classification of Subject Motion for Improved Reconstruction of Dynamic Magnetic Resonance Imaging Feiyu Chen Department of Electrical Engineering ABSTRACT Subject motion is a significant
More informationAlternative Statistical Methods for Bone Atlas Modelling
Alternative Statistical Methods for Bone Atlas Modelling Sharmishtaa Seshamani, Gouthami Chintalapani, Russell Taylor Department of Computer Science, Johns Hopkins University, Baltimore, MD Traditional
More informationUnsupervised learning in Vision
Chapter 7 Unsupervised learning in Vision The fields of Computer Vision and Machine Learning complement each other in a very natural way: the aim of the former is to extract useful information from visual
More informationIntroduction. Biomedical Image Analysis. Contents. Prof. Dr. Philippe Cattin. MIAC, University of Basel. Feb 22nd, of
Introduction Prof. Dr. Philippe Cattin MIAC, University of Basel Contents Abstract 1 Varia About Me About these Slides 2 My Research 2.1 Segmentation Segmentation of Facial Soft Tissues Segmentation of
More informationFMA901F: Machine Learning Lecture 3: Linear Models for Regression. Cristian Sminchisescu
FMA901F: Machine Learning Lecture 3: Linear Models for Regression Cristian Sminchisescu Machine Learning: Frequentist vs. Bayesian In the frequentist setting, we seek a fixed parameter (vector), with value(s)
More informationLearning based face hallucination techniques: A survey
Vol. 3 (2014-15) pp. 37-45. : A survey Premitha Premnath K Department of Computer Science & Engineering Vidya Academy of Science & Technology Thrissur - 680501, Kerala, India (email: premithakpnath@gmail.com)
More informationBlind Sparse Motion MRI with Linear Subpixel Interpolation
Blind Sparse Motion MRI with Linear Subpixel Interpolation Anita Möller, Marco Maaß, Alfred Mertins Institute for Signal Processing, University of Lübeck moeller@isip.uni-luebeck.de Abstract. Vital and
More informationComparison of Variational Bayes and Gibbs Sampling in Reconstruction of Missing Values with Probabilistic Principal Component Analysis
Comparison of Variational Bayes and Gibbs Sampling in Reconstruction of Missing Values with Probabilistic Principal Component Analysis Luis Gabriel De Alba Rivera Aalto University School of Science and
More informationFactorization with Missing and Noisy Data
Factorization with Missing and Noisy Data Carme Julià, Angel Sappa, Felipe Lumbreras, Joan Serrat, and Antonio López Computer Vision Center and Computer Science Department, Universitat Autònoma de Barcelona,
More informationIterative CT Reconstruction Using Curvelet-Based Regularization
Iterative CT Reconstruction Using Curvelet-Based Regularization Haibo Wu 1,2, Andreas Maier 1, Joachim Hornegger 1,2 1 Pattern Recognition Lab (LME), Department of Computer Science, 2 Graduate School in
More informationMotion artifact detection in four-dimensional computed tomography images
Motion artifact detection in four-dimensional computed tomography images G Bouilhol 1,, M Ayadi, R Pinho, S Rit 1, and D Sarrut 1, 1 University of Lyon, CREATIS; CNRS UMR 5; Inserm U144; INSA-Lyon; University
More informationFinite Element Simulation of Moving Targets in Radio Therapy
Finite Element Simulation of Moving Targets in Radio Therapy Pan Li, Gregor Remmert, Jürgen Biederer, Rolf Bendl Medical Physics, German Cancer Research Center, 69120 Heidelberg Email: pan.li@dkfz.de Abstract.
More informationSampling PCA, enhancing recovered missing values in large scale matrices. Luis Gabriel De Alba Rivera 80555S
Sampling PCA, enhancing recovered missing values in large scale matrices. Luis Gabriel De Alba Rivera 80555S May 2, 2009 Introduction Human preferences (the quality tags we put on things) are language
More informationGuide wire tracking in interventional radiology
Guide wire tracking in interventional radiology S.A.M. Baert,W.J. Niessen, E.H.W. Meijering, A.F. Frangi, M.A. Viergever Image Sciences Institute, University Medical Center Utrecht, rm E 01.334, P.O.Box
More informationFACE RECOGNITION USING SUPPORT VECTOR MACHINES
FACE RECOGNITION USING SUPPORT VECTOR MACHINES Ashwin Swaminathan ashwins@umd.edu ENEE633: Statistical and Neural Pattern Recognition Instructor : Prof. Rama Chellappa Project 2, Part (b) 1. INTRODUCTION
More informationSimultaneous Data Volume Reconstruction and Pose Estimation from Slice Samples
Simultaneous Data Volume Reconstruction and Pose Estimation from Slice Samples Manfred Georg *, Richard Souvenir, Andrew Hope, Robert Pless * Washington University in St. Louis * Saint Louis, MO 63130
More informationFast and robust extraction of surrogate respiratory signal from intra-operative liver ultrasound images
DOI 10.1007/s11548-013-0902-y ORIGINAL ARTICLE Fast and robust extraction of surrogate respiratory signal from intra-operative liver ultrasound images Jiaze Wu Cheng Li Su Huang Feng Liu Bien Soo Tan London
More informationReal-time Respiratory Motion Analysis Using Manifold Ray Casting of Volumetrically Fused Multi-View Range Imaging
Real-time Respiratory Motion Analysis Using Manifold Ray Casting of Volumetrically Fused Multi-View Range Imaging Jakob Wasza 1, Sebastian Bauer 1, Joachim Hornegger 1,2 1 Pattern Recognition Lab, Department
More informationAlgebraic Iterative Methods for Computed Tomography
Algebraic Iterative Methods for Computed Tomography Per Christian Hansen DTU Compute Department of Applied Mathematics and Computer Science Technical University of Denmark Per Christian Hansen Algebraic
More informationDetecting Burnscar from Hyperspectral Imagery via Sparse Representation with Low-Rank Interference
Detecting Burnscar from Hyperspectral Imagery via Sparse Representation with Low-Rank Interference Minh Dao 1, Xiang Xiang 1, Bulent Ayhan 2, Chiman Kwan 2, Trac D. Tran 1 Johns Hopkins Univeristy, 3400
More informationA Radiometry Tolerant Method for Direct 3D/2D Registration of Computed Tomography Data to X-ray Images
A Radiometry Tolerant Method for Direct 3D/2D Registration of Computed Tomography Data to X-ray Images Transfer Function Independent Registration Boris Peter Selby 1, Georgios Sakas 2, Stefan Walter 1,
More informationComputer Vision in a Non-Rigid World
2 nd Tutorial on Computer Vision in a Non-Rigid World Dr. Lourdes Agapito Prof. Adrien Bartoli Dr. Alessio Del Bue Non-rigid Structure from Motion Non-rigid structure from motion To extract non-rigid 3D
More informationCPSC 340: Machine Learning and Data Mining. Principal Component Analysis Fall 2016
CPSC 340: Machine Learning and Data Mining Principal Component Analysis Fall 2016 A2/Midterm: Admin Grades/solutions will be posted after class. Assignment 4: Posted, due November 14. Extra office hours:
More informationBayesian Spherical Wavelet Shrinkage: Applications to Shape Analysis
Bayesian Spherical Wavelet Shrinkage: Applications to Shape Analysis Xavier Le Faucheur a, Brani Vidakovic b and Allen Tannenbaum a a School of Electrical and Computer Engineering, b Department of Biomedical
More informationComputational Neuroanatomy
Computational Neuroanatomy John Ashburner john@fil.ion.ucl.ac.uk Smoothing Motion Correction Between Modality Co-registration Spatial Normalisation Segmentation Morphometry Overview fmri time-series kernel
More informationArtefakt-resistente Bewegungsschätzung für die bewegungskompensierte CT
Artefakt-resistente Bewegungsschätzung für die bewegungskompensierte CT Marcus Brehm 1,2, Thorsten Heußer 1, Pascal Paysan 3, Markus Oehlhafen 3, and Marc Kachelrieß 1,2 1 German Cancer Research Center
More informationCLASSIFICATION WITH RADIAL BASIS AND PROBABILISTIC NEURAL NETWORKS
CLASSIFICATION WITH RADIAL BASIS AND PROBABILISTIC NEURAL NETWORKS CHAPTER 4 CLASSIFICATION WITH RADIAL BASIS AND PROBABILISTIC NEURAL NETWORKS 4.1 Introduction Optical character recognition is one of
More informationMultidirectional 2DPCA Based Face Recognition System
Multidirectional 2DPCA Based Face Recognition System Shilpi Soni 1, Raj Kumar Sahu 2 1 M.E. Scholar, Department of E&Tc Engg, CSIT, Durg 2 Associate Professor, Department of E&Tc Engg, CSIT, Durg Email:
More informationImage Segmentation and Registration
Image Segmentation and Registration Dr. Christine Tanner (tanner@vision.ee.ethz.ch) Computer Vision Laboratory, ETH Zürich Dr. Verena Kaynig, Machine Learning Laboratory, ETH Zürich Outline Segmentation
More informationIEEE TRANSACTIONS ON VERY LARGE SCALE INTEGRATION (VLSI) SYSTEMS /$ IEEE
IEEE TRANSACTIONS ON VERY LARGE SCALE INTEGRATION (VLSI) SYSTEMS 1 Exploration of Heterogeneous FPGAs for Mapping Linear Projection Designs Christos-S. Bouganis, Member, IEEE, Iosifina Pournara, and Peter
More informationA fluence convolution method to account for respiratory motion in three-dimensional dose calculations of the liver: A Monte Carlo study
A fluence convolution method to account for respiratory motion in three-dimensional dose calculations of the liver: A Monte Carlo study Indrin J. Chetty, a) Mihaela Rosu, Neelam Tyagi, Lon H. Marsh, Daniel
More informationCHAPTER 9 INPAINTING USING SPARSE REPRESENTATION AND INVERSE DCT
CHAPTER 9 INPAINTING USING SPARSE REPRESENTATION AND INVERSE DCT 9.1 Introduction In the previous chapters the inpainting was considered as an iterative algorithm. PDE based method uses iterations to converge
More informationNon-Rigid Multimodal Medical Image Registration using Optical Flow and Gradient Orientation
M. HEINRICH et al.: MULTIMODAL REGISTRATION USING GRADIENT ORIENTATION 1 Non-Rigid Multimodal Medical Image Registration using Optical Flow and Gradient Orientation Mattias P. Heinrich 1 mattias.heinrich@eng.ox.ac.uk
More informationHierarchical Shape Statistical Model for Segmentation of Lung Fields in Chest Radiographs
Hierarchical Shape Statistical Model for Segmentation of Lung Fields in Chest Radiographs Yonghong Shi 1 and Dinggang Shen 2,*1 1 Digital Medical Research Center, Fudan University, Shanghai, 232, China
More informationEstimating the Information Rate of Noisy Two-Dimensional Constrained Channels
Estimating the Information Rate of Noisy Two-Dimensional Constrained Channels Mehdi Molkaraie and Hans-Andrea Loeliger Dept. of Information Technology and Electrical Engineering ETH Zurich, Switzerland
More informationDimension Reduction CS534
Dimension Reduction CS534 Why dimension reduction? High dimensionality large number of features E.g., documents represented by thousands of words, millions of bigrams Images represented by thousands of
More informationInteractive Deformable Registration Visualization and Analysis of 4D Computed Tomography
Interactive Deformable Registration Visualization and Analysis of 4D Computed Tomography Burak Erem 1, Gregory C. Sharp 2, Ziji Wu 2, and David Kaeli 1 1 Department of Electrical and Computer Engineering,
More informationCS 565 Computer Vision. Nazar Khan PUCIT Lectures 15 and 16: Optic Flow
CS 565 Computer Vision Nazar Khan PUCIT Lectures 15 and 16: Optic Flow Introduction Basic Problem given: image sequence f(x, y, z), where (x, y) specifies the location and z denotes time wanted: displacement
More informationCombination of Markerless Surrogates for Motion Estimation in Radiation Therapy
Combination of Markerless Surrogates for Motion Estimation in Radiation Therapy CARS 2016 T. Geimer, M. Unberath, O. Taubmann, C. Bert, A. Maier June 24, 2016 Pattern Recognition Lab (CS 5) FAU Erlangen-Nu
More informationAn efficient algorithm for sparse PCA
An efficient algorithm for sparse PCA Yunlong He Georgia Institute of Technology School of Mathematics heyunlong@gatech.edu Renato D.C. Monteiro Georgia Institute of Technology School of Industrial & System
More informationABSTRACT 1. INTRODUCTION 2. METHODS
Finding Seeds for Segmentation Using Statistical Fusion Fangxu Xing *a, Andrew J. Asman b, Jerry L. Prince a,c, Bennett A. Landman b,c,d a Department of Electrical and Computer Engineering, Johns Hopkins
More informationRespiratory Motion Compensation for Simultaneous PET/MR Based on Strongly Undersampled Radial MR Data
Respiratory Motion Compensation for Simultaneous PET/MR Based on Strongly Undersampled Radial MR Data Christopher M Rank 1, Thorsten Heußer 1, Andreas Wetscherek 1, and Marc Kachelrieß 1 1 German Cancer
More informationNeural Networks Based Time-Delay Estimation using DCT Coefficients
American Journal of Applied Sciences 6 (4): 73-78, 9 ISSN 1546-939 9 Science Publications Neural Networks Based Time-Delay Estimation using DCT Coefficients Samir J. Shaltaf and Ahmad A. Mohammad Department
More informationLocally Weighted Least Squares Regression for Image Denoising, Reconstruction and Up-sampling
Locally Weighted Least Squares Regression for Image Denoising, Reconstruction and Up-sampling Moritz Baecher May 15, 29 1 Introduction Edge-preserving smoothing and super-resolution are classic and important
More information1498. End-effector vibrations reduction in trajectory tracking for mobile manipulator
1498. End-effector vibrations reduction in trajectory tracking for mobile manipulator G. Pajak University of Zielona Gora, Faculty of Mechanical Engineering, Zielona Góra, Poland E-mail: g.pajak@iizp.uz.zgora.pl
More informationHuman Motion Synthesis by Motion Manifold Learning and Motion Primitive Segmentation
Human Motion Synthesis by Motion Manifold Learning and Motion Primitive Segmentation Chan-Su Lee and Ahmed Elgammal Rutgers University, Piscataway, NJ, USA {chansu, elgammal}@cs.rutgers.edu Abstract. We
More informationICA vs. PCA Active Appearance Models: Application to Cardiac MR Segmentation
ICA vs. PCA Active Appearance Models: Application to Cardiac MR Segmentation M. Üzümcü 1, A.F. Frangi 2, M. Sonka 3, J.H.C. Reiber 1, B.P.F. Lelieveldt 1 1 Div. of Image Processing, Dept. of Radiology
More informationDynamic Texture with Fourier Descriptors
B. Abraham, O. Camps and M. Sznaier: Dynamic texture with Fourier descriptors. In Texture 2005: Proceedings of the 4th International Workshop on Texture Analysis and Synthesis, pp. 53 58, 2005. Dynamic
More informationSUPER RESOLUTION RECONSTRUCTION OF CARDIAC MRI USING COUPLED DICTIONARY LEARNING
SUPER RESOLUTION RECONSTRUCTION OF CARDIAC MRI USING COUPLED DICTIONARY LEARNING Abstract M. Vinod Kumar (M.tech) 1 V. Gurumurthy Associate Professor, M.Tech (Ph.D) 2 Dr.M. Narayana, Professor, Head of
More informationLearning a Manifold as an Atlas Supplementary Material
Learning a Manifold as an Atlas Supplementary Material Nikolaos Pitelis Chris Russell School of EECS, Queen Mary, University of London [nikolaos.pitelis,chrisr,lourdes]@eecs.qmul.ac.uk Lourdes Agapito
More informationFOREWORD TO THE SPECIAL ISSUE ON MOTION DETECTION AND COMPENSATION
Philips J. Res. 51 (1998) 197-201 FOREWORD TO THE SPECIAL ISSUE ON MOTION DETECTION AND COMPENSATION This special issue of Philips Journalof Research includes a number of papers presented at a Philips
More informationIntraoperative Prostate Tracking with Slice-to-Volume Registration in MR
Intraoperative Prostate Tracking with Slice-to-Volume Registration in MR Sean Gill a, Purang Abolmaesumi a,b, Siddharth Vikal a, Parvin Mousavi a and Gabor Fichtinger a,b,* (a) School of Computing, Queen
More informationCOMPUTATIONAL STATISTICS UNSUPERVISED LEARNING
COMPUTATIONAL STATISTICS UNSUPERVISED LEARNING Luca Bortolussi Department of Mathematics and Geosciences University of Trieste Office 238, third floor, H2bis luca@dmi.units.it Trieste, Winter Semester
More informationSalient Region Detection and Segmentation in Images using Dynamic Mode Decomposition
Salient Region Detection and Segmentation in Images using Dynamic Mode Decomposition Sikha O K 1, Sachin Kumar S 2, K P Soman 2 1 Department of Computer Science 2 Centre for Computational Engineering and
More information1 Introduction Motivation and Aims Functional Imaging Computational Neuroanatomy... 12
Contents 1 Introduction 10 1.1 Motivation and Aims....... 10 1.1.1 Functional Imaging.... 10 1.1.2 Computational Neuroanatomy... 12 1.2 Overview of Chapters... 14 2 Rigid Body Registration 18 2.1 Introduction.....
More informationFace Tracking. Synonyms. Definition. Main Body Text. Amit K. Roy-Chowdhury and Yilei Xu. Facial Motion Estimation
Face Tracking Amit K. Roy-Chowdhury and Yilei Xu Department of Electrical Engineering, University of California, Riverside, CA 92521, USA {amitrc,yxu}@ee.ucr.edu Synonyms Facial Motion Estimation Definition
More informationMarkerless Real-Time Target Region Tracking: Application to Frameless Stereotactic Radiosurgery
Markerless Real-Time Target Region Tracking: Application to Frameless Stereotactic Radiosurgery T. Rohlfing,, J. Denzler 3, D.B. Russakoff, Ch. Gräßl 4, and C.R. Maurer, Jr. Neuroscience Program, SRI International,
More informationFully Automatic Model Creation for Object Localization utilizing the Generalized Hough Transform
Fully Automatic Model Creation for Object Localization utilizing the Generalized Hough Transform Heike Ruppertshofen 1,2,3, Cristian Lorenz 2, Peter Beyerlein 4, Zein Salah 3, Georg Rose 3, Hauke Schramm
More informationOutlier Pursuit: Robust PCA and Collaborative Filtering
Outlier Pursuit: Robust PCA and Collaborative Filtering Huan Xu Dept. of Mechanical Engineering & Dept. of Mathematics National University of Singapore Joint w/ Constantine Caramanis, Yudong Chen, Sujay
More informationHYBRID TRANSFORMATION TECHNIQUE FOR IMAGE COMPRESSION
31 st July 01. Vol. 41 No. 005-01 JATIT & LLS. All rights reserved. ISSN: 199-8645 www.jatit.org E-ISSN: 1817-3195 HYBRID TRANSFORMATION TECHNIQUE FOR IMAGE COMPRESSION 1 SRIRAM.B, THIYAGARAJAN.S 1, Student,
More information22 October, 2012 MVA ENS Cachan. Lecture 5: Introduction to generative models Iasonas Kokkinos
Machine Learning for Computer Vision 1 22 October, 2012 MVA ENS Cachan Lecture 5: Introduction to generative models Iasonas Kokkinos Iasonas.kokkinos@ecp.fr Center for Visual Computing Ecole Centrale Paris
More informationNeural Networks. Theory And Practice. Marco Del Vecchio 19/07/2017. Warwick Manufacturing Group University of Warwick
Neural Networks Theory And Practice Marco Del Vecchio marco@delvecchiomarco.com Warwick Manufacturing Group University of Warwick 19/07/2017 Outline I 1 Introduction 2 Linear Regression Models 3 Linear
More informationAlgebraic Iterative Methods for Computed Tomography
Algebraic Iterative Methods for Computed Tomography Per Christian Hansen DTU Compute Department of Applied Mathematics and Computer Science Technical University of Denmark Per Christian Hansen Algebraic
More informationarxiv: v2 [stat.ml] 5 Nov 2018
Kernel Distillation for Fast Gaussian Processes Prediction arxiv:1801.10273v2 [stat.ml] 5 Nov 2018 Congzheng Song Cornell Tech cs2296@cornell.edu Abstract Yiming Sun Cornell University ys784@cornell.edu
More informationPathology Hinting as the Combination of Automatic Segmentation with a Statistical Shape Model
Pathology Hinting as the Combination of Automatic Segmentation with a Statistical Shape Model Pascal A. Dufour 12,HannanAbdillahi 3, Lala Ceklic 3,Ute Wolf-Schnurrbusch 23,JensKowal 12 1 ARTORG Center
More information10-701/15-781, Fall 2006, Final
-7/-78, Fall 6, Final Dec, :pm-8:pm There are 9 questions in this exam ( pages including this cover sheet). If you need more room to work out your answer to a question, use the back of the page and clearly
More informationMULTIVARIATE TEXTURE DISCRIMINATION USING A PRINCIPAL GEODESIC CLASSIFIER
MULTIVARIATE TEXTURE DISCRIMINATION USING A PRINCIPAL GEODESIC CLASSIFIER A.Shabbir 1, 2 and G.Verdoolaege 1, 3 1 Department of Applied Physics, Ghent University, B-9000 Ghent, Belgium 2 Max Planck Institute
More informationNonrigid Registration using Free-Form Deformations
Nonrigid Registration using Free-Form Deformations Hongchang Peng April 20th Paper Presented: Rueckert et al., TMI 1999: Nonrigid registration using freeform deformations: Application to breast MR images
More informationNON-RIGID 3D SHAPE RECOVERY USING STEREO FACTORIZATION
NON-IGID D SHAPE ECOVEY USING STEEO FACTOIZATION Alessio Del Bue ourdes Agapito Department of Computer Science Queen Mary, University of ondon ondon, E1 NS, UK falessio,lourdesg@dcsqmulacuk ABSTACT In
More informationRespiratory Motion Compensation for C-arm CT Liver Imaging
Respiratory Motion Compensation for C-arm CT Liver Imaging Aline Sindel 1, Marco Bögel 1,2, Andreas Maier 1,2, Rebecca Fahrig 3, Joachim Hornegger 1,2, Arnd Dörfler 4 1 Pattern Recognition Lab, FAU Erlangen-Nürnberg
More informationLearning Algorithms for Medical Image Analysis. Matteo Santoro slipguru
Learning Algorithms for Medical Image Analysis Matteo Santoro slipguru santoro@disi.unige.it June 8, 2010 Outline 1. learning-based strategies for quantitative image analysis 2. automatic annotation of
More informationSUPER-RESOLUTION RECONSTRUCTION ALGORITHM FOR BASED ON COUPLED DICTIONARY LEARNING CARDIAC MRI RECONSTRUCTION
SUPER-RESOLUTION RECONSTRUCTION ALGORITHM FOR BASED ON COUPLED DICTIONARY LEARNING CARDIAC MRI RECONSTRUCTION M.PRATHAP KUMAR, Email Id: Prathap.Macherla@Gmail.Com J.VENKATA LAKSHMI, M.Tech, Asst.Prof,
More informationNetwork Traffic Measurements and Analysis
DEIB - Politecnico di Milano Fall, 2017 Introduction Often, we have only a set of features x = x 1, x 2,, x n, but no associated response y. Therefore we are not interested in prediction nor classification,
More informationA System Identification Approach to Estimating a Dynamic Model of Head Motion for MRI Motion Correction
A System Identification Approach to Estimating a Dynamic Model of Head Motion for MRI Motion Correction Burak Erem 1, Onur Afacan 1, Ali Gholipour 1, Sanjay P. Prabhu 2, and Simon K. Warfield 1 1 Computational
More informationBildverarbeitung für die Medizin 2007
Bildverarbeitung für die Medizin 2007 Image Registration with Local Rigidity Constraints Jan Modersitzki Institute of Mathematics, University of Lübeck, Wallstraße 40, D-23560 Lübeck 1 Summary Registration
More informationRobust Realignment of fmri Time Series Data
Robust Realignment of fmri Time Series Data Ben Dodson bjdodson@stanford.edu Olafur Gudmundsson olafurg@stanford.edu December 12, 2008 Abstract FMRI data has become an increasingly popular source for exploring
More informationProbabilistic Tracking and Model-based Segmentation of 3D Tubular Structures
Probabilistic Tracking and Model-based Segmentation of 3D Tubular Structures Stefan Wörz, William J. Godinez, Karl Rohr University of Heidelberg, BIOQUANT, IPMB, and DKFZ Heidelberg, Dept. Bioinformatics
More informationA Robust Image Zero-Watermarking Algorithm Based on DWT and PCA
A Robust Image Zero-Watermarking Algorithm Based on DWT and PCA Xiaoxu Leng, Jun Xiao, and Ying Wang Graduate University of Chinese Academy of Sciences, 100049 Beijing, China lengxiaoxu@163.com, {xiaojun,ywang}@gucas.ac.cn
More informationCompressed Sensing Algorithm for Real-Time Doppler Ultrasound Image Reconstruction
Mathematical Modelling and Applications 2017; 2(6): 75-80 http://www.sciencepublishinggroup.com/j/mma doi: 10.11648/j.mma.20170206.14 ISSN: 2575-1786 (Print); ISSN: 2575-1794 (Online) Compressed Sensing
More informationCompression of Stereo Images using a Huffman-Zip Scheme
Compression of Stereo Images using a Huffman-Zip Scheme John Hamann, Vickey Yeh Department of Electrical Engineering, Stanford University Stanford, CA 94304 jhamann@stanford.edu, vickey@stanford.edu Abstract
More informationAdaptive Multiple-Frame Image Super- Resolution Based on U-Curve
Adaptive Multiple-Frame Image Super- Resolution Based on U-Curve IEEE Transaction on Image Processing, Vol. 19, No. 12, 2010 Qiangqiang Yuan, Liangpei Zhang, Huanfeng Shen, and Pingxiang Li Presented by
More informationarxiv: v1 [cs.cv] 6 Jun 2017
Volume Calculation of CT lung Lesions based on Halton Low-discrepancy Sequences Liansheng Wang a, Shusheng Li a, and Shuo Li b a Department of Computer Science, Xiamen University, Xiamen, China b Dept.
More informationAN EXAMINING FACE RECOGNITION BY LOCAL DIRECTIONAL NUMBER PATTERN (Image Processing)
AN EXAMINING FACE RECOGNITION BY LOCAL DIRECTIONAL NUMBER PATTERN (Image Processing) J.Nithya 1, P.Sathyasutha2 1,2 Assistant Professor,Gnanamani College of Engineering, Namakkal, Tamil Nadu, India ABSTRACT
More informationMachine Learning A W 1sst KU. b) [1 P] Give an example for a probability distributions P (A, B, C) that disproves
Machine Learning A 708.064 11W 1sst KU Exercises Problems marked with * are optional. 1 Conditional Independence I [2 P] a) [1 P] Give an example for a probability distribution P (A, B, C) that disproves
More informationProbabilistic Graphical Models
Overview of Part Two Probabilistic Graphical Models Part Two: Inference and Learning Christopher M. Bishop Exact inference and the junction tree MCMC Variational methods and EM Example General variational
More informationCompressed Sensing Reconstructions for Dynamic Contrast Enhanced MRI
1 Compressed Sensing Reconstructions for Dynamic Contrast Enhanced MRI Kevin T. Looby klooby@stanford.edu ABSTRACT The temporal resolution necessary for dynamic contrast enhanced (DCE) magnetic resonance
More informationSingular Value Decomposition, and Application to Recommender Systems
Singular Value Decomposition, and Application to Recommender Systems CSE 6363 Machine Learning Vassilis Athitsos Computer Science and Engineering Department University of Texas at Arlington 1 Recommendation
More informationPERFORMANCE ANALYSIS OF INTEGER DCT OF DIFFERENT BLOCK SIZES USED IN H.264, AVS CHINA AND WMV9.
EE 5359: MULTIMEDIA PROCESSING PROJECT PERFORMANCE ANALYSIS OF INTEGER DCT OF DIFFERENT BLOCK SIZES USED IN H.264, AVS CHINA AND WMV9. Guided by Dr. K.R. Rao Presented by: Suvinda Mudigere Srikantaiah
More informationPerception and Action using Multilinear Forms
Perception and Action using Multilinear Forms Anders Heyden, Gunnar Sparr, Kalle Åström Dept of Mathematics, Lund University Box 118, S-221 00 Lund, Sweden email: {heyden,gunnar,kalle}@maths.lth.se Abstract
More informationSelective vs. Global Recovery of Rigid and Non-Rigid Motion
Selective vs. Global Recovery of Rigid and Non-Rigid Motion Sami Romdhani University of Basel Departement Informatik Bernoullistrasse 16 CH-456 Basel Switzerland Nikolaos Canterakis Albert-Ludwigs-Universitt
More information