Proton Therapy: 3D Reconstructions from 2D Images
|
|
- Jane Cain
- 5 years ago
- Views:
Transcription
1 Proton Therapy: 3D Reconstructions from 2D Images Ben Herbst Department of Applied Mathematics University of Stellenbosch Colorado School of Mines 16 September 2005
2 2 Collaborators Neil Muller Barry Sherlock Karin Hunter Brian O Kennedy Chris Venter Lourenco Magaia Francois Malan Pieter Rautenbach
3 Overview Proton Therapy and Patient Positioning Stereo Reconstruction 3D Reconstruction from video sequences 3D Reconstruction, single image
4 4 Patient Positioning Nuclear Accelerator
5 Physics Protons from hydrogen atoms, separated from water: Boosted to 20MEV 5 million trips around the accelerator ring (600,000,000rpm) 1 speed of light 2 In live tissue: destroys DNA through ionization Bragg peak implies almost perfect Localization
6 6 Proton Therapy Bragg Peak Robert R Wilson 1946
7 First Patients: facilities world wide 3 in USA 1 facility near Stellenbosch Background Information
8 8 Aligning Proton Beam Many 6000 pound magnets that consume 3000 amperes each Precisely aligned by three-story high 200 ton gantries Aim the beam within a fraction of a millimeter Accelerator pumps a million watts in and out of power plant Deliver 1 5 watt, for 1 3 s Need a cheaper solution
9 Patient Positioning Position the Patient
10 Patient Positioning 10 (a) Mask (b) Reflective markers Patient Mask
11 CT Calibration (c) CT Scanner (d) CT Scanner CT Scanner
12 Patient Positioning 12 (e) Vault (f) Camera Treatment Vault: In development
13 Perspective Cameras Need to know: geometry - camera center, relative orientation corresponding points lines in 3D do not intersect
14 2D Projective Geometry 14 Equation for a line: For line becomes x = ax + by + c = 0 x y 1, l = x T l = 0 x and l are homogeneous vectors: identify all cx, cl all real c 0 x, l P 2 R 2 P 2 : P 2 R 2 : [ x y x y z ] a b c x y 1 [ x/z y/z, ]
15 Intersection of Lines The intersection of l 1 and l 2 : x = l 1 l 2 The line through x 1 and x 2 : l = x 1 x 2 Example l 1 : x 1 = 0 l 2 : x 2 = 0 Intersection x = = Two parallel lines intersect in well-defined point
16 Projective Transformation 16 H : P 2 P 2 preserves straight lines H is a 3 3 nonsingular, homogeneous matrix (g) Transformed image (h) Corrected image
17 Perspective Camera transformation Project X in object (world coordinates), to x on CCD (pixel coordinates) where x = PX P = KR [ I C ] P is 3 4, homogeneous, rank 3 P C = 0 ( C camera center in world coordinates)
18 Camera Calibration 18 (i) Lego Calibration Object (j) Real Calibration Pattern Calibration Object
19 Triangulation For x = PX and P C = 0 X = P + x + λ C For stereo camera pair, this gives X. Find corresponding points from epipolar lines: x T Fx = 0 F is fundamental matrix: calculated from camera matrices F = e P P +
20 Lines do not intersect in 3D 20 Minimize d 2 + d 2 subject to x T Fx = 0 Linearize SVD or Lagrangian multipliers
21 3D Reconstruction from a video sequence
22 Shackleton (Endurance ) 22 Frank Worseley(navigator): On the thirteenth day we were getting nearer to our destination. If we made the tragic mistake of passing it we could never retrace our way on account of the winds and the currents, it therefore became essential that I should get observations. But the morning was foggy, and if you cannot see the horizon it is impossible to measure the altitude of the sun to establish your position. Now, the nearer your eye is to the surface of the sea, the nearer is the horizon. So I adopted the expedient of kneeling on the stones in the bottom of the boat, and by this means succeeded in taking a rough observation. It would have been a bold assumption to say that it was a correct one; but is was the best I could do, and we had to trust it. Two observations are necessary, however, to fix your position, and my troubles were far from over; for at noon, when I wanted to observe our latitude, I found conditions equally difficult. The fog, which before had been on a level with us and therefore did not altogether obscure the sun, had now risen above us and was hovering between the sun and ourselves, so that all I could see was a dim blur. I measured to the centre of this ten times, using the mean of these observations as the sun s altitude.
23 Conditional Densities p(x z) z 5 Given noisy measurement z 0, variance σ 2 z 0 of position x: Conditional density: p(x z 0 ) = N (z 0, σ 2 z 0 ) Best estimate of x: ^x(1 0) = z 0 (estimate at t 1 given measurement z 0 at t 0 )
24 Conditional Densities 24 σ z1 σ z0 z 0 z 1 Given second noisy measurement z 1, variance σ 2 z 1 of x: Conditional density: p(x z 0, z 1 ) = N (x z, σ 2 z)
25 Conditional Densities σ σ z1 σ z0 z 0 µ z 1 Conditional density: p(x z 0, z 1 ) = N (x z, σ 2 z ) [ σ 2 z 1 x z = σ 2 z 0 + σ 2 z 1 1 = σ 2 z 0 σ 2 z 1 σ 2 z ] z 0 + [ σ 2 z 0 σ 2 z 0 + σ 2 z 1 ] z 1 σ 2 z smaller than either σ2 z 0 or σ 2 z 1 less accurate estimates can be used to improve estimate
26 Dynamics 26 dx dt = u + w u: velocity w: noise, variance σ 2 w. Combine best knowledge of position at t 1 with dynamics to get best estimate at t 2 : x(2 1) = x(1 1) + (t 2 t 1 )u. variance σ 2 x(2 1) = σ 2 x(1 1) + (t 2 t 1 )σ 2 w. Disaster, unless measurement z 2 is obtained
27 Dynamics Given: x(2 1) = x(1 1) + (t 2 t 1 )u, with σ 2 x(2 1) = σ 2 x(1 1) + (t 2 t 1 )σ 2 w new measurement z 2, variance σ 2 z 2 Combine and x(2 2) = x(2 1) + K(t 2 )(z 2 x(2 1)), σ 2 x (2 2) = σ2 x (2 1) K(t 2)σ 2 x (2 1). with σ 2 K(t 2 ) = x(2 1) σ 2 x (2 1) + σ2 z 2 If measurement noise σ z2 is large, then K(t 2 ) is small little confidence in noisy measurement If σ z2 then K(t 2 ) 0 ignores measurement If dynamic system noise σ w is large, so is K(t 2 ) little confidence in dynamic system if σ 2 w then K(t 2) 1 and x(2 2) z 2 all information obtained from measurement
28 Kalman Filter 28 Dynamics: x k+1 = A k x k + w k Measurement: z k+1 = H k+1 x k+1 + ν k+1 Noise w k and ν k have normal distributions: N (0, Q k ) and N (0, R k ) Dynamic update Given x(k k) and covariance P(k k) at time t k prior to the measurement at t k+1 : (a) Predicted state at t k+1 : Covariance: (b) Predicted measurement: Covariance: x(k + 1 k) = A k x(k k) P(k + 1 k) = A k P(k k)a T k + Q k ẑ(k + 1 k) = H k x(k + 1 k) S(k + 1 k) = H k+1 P(k + 1 k)h T k+1 + R k+1
29 Kalman Filter Dynamics: x k+1 = A k x k + w k Measurement: z k+1 = H k+1 x k+1 + ν k+1 Noise w k and ν k have normal distributions: N (0, Q k ) and N (0, R k ) Measurement update Given x(k + 1 k) and P(k + 1 k), and measurement z k+1 at t k+1 : ẑ(k + 1 k) and S(k + 1 k) (a) Kalman gain: W(k + 1) = P(k + 1 k)h T k+1s(k + 1 k) 1 (b) State update: x(k + 1 k + 1) = x(k + 1 k) + W(k + 1)(z k+1 ẑ(k + 1 k)) Covariance: P(k + 1 k + 1) = P(k + 1 k) W(k + 1)S(k + 1 k)w(k + 1) T.
30 Application 30 State variables: Structure: X j, j = 1,..., N (time independent) Motion: Rotation R and translation t Dynamic update: x k+1 = x k + w k Measurement: (Nonlinear) Error covariances: update, measurement Need nonlinear extension of Kalman, EKF, UKF, etc
1 Kalman Filter Introduction
1 Kalman Filter Introduction You should first read Chapter 1 of Stochastic models, estimation, and control: Volume 1 by Peter S. Maybec (available here). 1.1 Explanation of Equations (1-3) and (1-4) Equation
More informationCS 532: 3D Computer Vision 6 th Set of Notes
1 CS 532: 3D Computer Vision 6 th Set of Notes Instructor: Philippos Mordohai Webpage: www.cs.stevens.edu/~mordohai E-mail: Philippos.Mordohai@stevens.edu Office: Lieb 215 Lecture Outline Intro to Covariance
More informationTime Series Analysis
Time Series Analysis hm@imm.dtu.dk Informatics and Mathematical Modelling Technical University of Denmark DK-2800 Kgs. Lyngby 1 Outline of the lecture State space models, 1st part: Model: Sec. 10.1 The
More informationImage Alignment and Mosaicing Feature Tracking and the Kalman Filter
Image Alignment and Mosaicing Feature Tracking and the Kalman Filter Image Alignment Applications Local alignment: Tracking Stereo Global alignment: Camera jitter elimination Image enhancement Panoramic
More informationKalman Filter. Predict: Update: x k k 1 = F k x k 1 k 1 + B k u k P k k 1 = F k P k 1 k 1 F T k + Q
Kalman Filter Kalman Filter Predict: x k k 1 = F k x k 1 k 1 + B k u k P k k 1 = F k P k 1 k 1 F T k + Q Update: K = P k k 1 Hk T (H k P k k 1 Hk T + R) 1 x k k = x k k 1 + K(z k H k x k k 1 ) P k k =(I
More informationFinal Exam Due on Sunday 05/06
Final Exam Due on Sunday 05/06 The exam should be completed individually without collaboration. However, you are permitted to consult with the textbooks, notes, slides and even internet resources. If you
More informationMultiple View Geometry in Computer Vision
Multiple View Geometry in Computer Vision Prasanna Sahoo Department of Mathematics University of Louisville 1 Scene Planes & Homographies Lecture 19 March 24, 2005 2 In our last lecture, we examined various
More informationL06. LINEAR KALMAN FILTERS. NA568 Mobile Robotics: Methods & Algorithms
L06. LINEAR KALMAN FILTERS NA568 Mobile Robotics: Methods & Algorithms 2 PS2 is out! Landmark-based Localization: EKF, UKF, PF Today s Lecture Minimum Mean Square Error (MMSE) Linear Kalman Filter Gaussian
More informationEE Camera & Image Formation
Electric Electronic Engineering Bogazici University February 21, 2018 Introduction Introduction Camera models Goal: To understand the image acquisition process. Function of the camera Similar to that of
More informationMiscellaneous. Regarding reading materials. Again, ask questions (if you have) and ask them earlier
Miscellaneous Regarding reading materials Reading materials will be provided as needed If no assigned reading, it means I think the material from class is sufficient Should be enough for you to do your
More informationThe Kalman Filter ImPr Talk
The Kalman Filter ImPr Talk Ged Ridgway Centre for Medical Image Computing November, 2006 Outline What is the Kalman Filter? State Space Models Kalman Filter Overview Bayesian Updating of Estimates Kalman
More informationState Estimation of Linear and Nonlinear Dynamic Systems
State Estimation of Linear and Nonlinear Dynamic Systems Part I: Linear Systems with Gaussian Noise James B. Rawlings and Fernando V. Lima Department of Chemical and Biological Engineering University of
More informationthe robot in its current estimated position and orientation (also include a point at the reference point of the robot)
CSCI 4190 Introduction to Robotic Algorithms, Spring 006 Assignment : out February 13, due February 3 and March Localization and the extended Kalman filter In this assignment, you will write a program
More informationEKF, UKF. Pieter Abbeel UC Berkeley EECS. Many slides adapted from Thrun, Burgard and Fox, Probabilistic Robotics
EKF, UKF Pieter Abbeel UC Berkeley EECS Many slides adapted from Thrun, Burgard and Fox, Probabilistic Robotics Kalman Filter Kalman Filter = special case of a Bayes filter with dynamics model and sensory
More informationEKF, UKF. Pieter Abbeel UC Berkeley EECS. Many slides adapted from Thrun, Burgard and Fox, Probabilistic Robotics
EKF, UKF Pieter Abbeel UC Berkeley EECS Many slides adapted from Thrun, Burgard and Fox, Probabilistic Robotics Kalman Filter Kalman Filter = special case of a Bayes filter with dynamics model and sensory
More informationLecture 5. Epipolar Geometry. Professor Silvio Savarese Computational Vision and Geometry Lab. 21-Jan-15. Lecture 5 - Silvio Savarese
Lecture 5 Epipolar Geometry Professor Silvio Savarese Computational Vision and Geometry Lab Silvio Savarese Lecture 5-21-Jan-15 Lecture 5 Epipolar Geometry Why is stereo useful? Epipolar constraints Essential
More informationModel based optimization and estimation of the field map during the breakdown phase in the ITER tokamak
Model based optimization and estimation of the field map during the breakdown phase in the ITER tokamak Roberto Ambrosino 1 Gianmaria De Tommasi 2 Massimiliano Mattei 3 Alfredo Pironti 2 1 CREATE, Università
More informationUsing the Kalman Filter to Estimate the State of a Maneuvering Aircraft
1 Using the Kalman Filter to Estimate the State of a Maneuvering Aircraft K. Meier and A. Desai Abstract Using sensors that only measure the bearing angle and range of an aircraft, a Kalman filter is implemented
More informationTowards efficient and accurate particle transport simulation in medical applications
Towards efficient and accurate particle transport simulation in medical applications L. Grzanka1,2, M. Kłodowska1, N. Mojżeszek1, N. Bassler3 1 Cyclotron Centre Bronowice, Institute of Nuclear Physics
More informationMulti-Frame Factorization Techniques
Multi-Frame Factorization Techniques Suppose { x j,n } J,N j=1,n=1 is a set of corresponding image coordinates, where the index n = 1,...,N refers to the n th scene point and j = 1,..., J refers to the
More informationBayes Filter Reminder. Kalman Filter Localization. Properties of Gaussians. Gaussians. Prediction. Correction. σ 2. Univariate. 1 2πσ e.
Kalman Filter Localization Bayes Filter Reminder Prediction Correction Gaussians p(x) ~ N(µ,σ 2 ) : Properties of Gaussians Univariate p(x) = 1 1 2πσ e 2 (x µ) 2 σ 2 µ Univariate -σ σ Multivariate µ Multivariate
More informationData assimilation with and without a model
Data assimilation with and without a model Tim Sauer George Mason University Parameter estimation and UQ U. Pittsburgh Mar. 5, 2017 Partially supported by NSF Most of this work is due to: Tyrus Berry,
More informationKalman Filtering in a Mass-Spring System
Department of Mathematics Center for Quantitative Sciences in Biomedicine North Carolina State University July 30, 2016 Practical Overview Goal: Implement Kalman filter for linear system 1. Solve a linear
More informationUNC Charlotte 2004 Algebra with solutions
with solutions March 8, 2004 1. Let z denote the real number solution to of the digits of z? (A) 13 (B) 14 (C) 15 (D) 16 (E) 17 3 + x 1 = 5. What is the sum Solution: E. Square both sides twice to get
More informationData assimilation with and without a model
Data assimilation with and without a model Tyrus Berry George Mason University NJIT Feb. 28, 2017 Postdoc supported by NSF This work is in collaboration with: Tim Sauer, GMU Franz Hamilton, Postdoc, NCSU
More informationQuaternion based Extended Kalman Filter
Quaternion based Extended Kalman Filter, Sergio Montenegro About this lecture General introduction to rotations and quaternions. Introduction to Kalman Filter for Attitude Estimation How to implement and
More informationRadiation protection issues in proton therapy
Protons IMRT Tony Lomax, Centre for Proton Radiotherapy, Paul Scherrer Institute, Switzerland Overview of presentation 1. Proton therapy: An overview 2. Radiation protection issues: Staff 3. Radiation
More informationVideo and Motion Analysis Computer Vision Carnegie Mellon University (Kris Kitani)
Video and Motion Analysis 16-385 Computer Vision Carnegie Mellon University (Kris Kitani) Optical flow used for feature tracking on a drone Interpolated optical flow used for super slow-mo optical flow
More informationState Observers and the Kalman filter
Modelling and Control of Dynamic Systems State Observers and the Kalman filter Prof. Oreste S. Bursi University of Trento Page 1 Feedback System State variable feedback system: Control feedback law:u =
More information2D Image Processing (Extended) Kalman and particle filter
2D Image Processing (Extended) Kalman and particle filter Prof. Didier Stricker Dr. Gabriele Bleser Kaiserlautern University http://ags.cs.uni-kl.de/ DFKI Deutsches Forschungszentrum für Künstliche Intelligenz
More informationLecture 2: From Linear Regression to Kalman Filter and Beyond
Lecture 2: From Linear Regression to Kalman Filter and Beyond Department of Biomedical Engineering and Computational Science Aalto University January 26, 2012 Contents 1 Batch and Recursive Estimation
More informationDielectric Wall Accelerator (DWA) and Distal Edge Tracking Proton Delivery System Rock Mackie Professor Dept of Medical Physics UW Madison Co-Founder
Dielectric Wall Accelerator (DWA) and Distal Edge Tracking Proton Delivery System Rock Mackie Professor Dept of Medical Physics UW Madison Co-Founder and Chairman of the Board or TomoTherapy Inc I have
More informationCamera calibration Triangulation
Triangulation Perspective projection in homogenous coordinates ~x img I 0 apple R t 0 T 1 ~x w ~x img R t ~x w Matrix transformations in 2D ~x img K R t ~x w K = 2 3 1 0 t u 40 1 t v 5 0 0 1 Translation
More informationOnline monitoring of MPC disturbance models using closed-loop data
Online monitoring of MPC disturbance models using closed-loop data Brian J. Odelson and James B. Rawlings Department of Chemical Engineering University of Wisconsin-Madison Online Optimization Based Identification
More informationComputer Vision Motion
Computer Vision Motion Professor Hager http://www.cs.jhu.edu/~hager 12/1/12 CS 461, Copyright G.D. Hager Outline From Stereo to Motion The motion field and optical flow (2D motion) Factorization methods
More informationState Estimation of Linear and Nonlinear Dynamic Systems
State Estimation of Linear and Nonlinear Dynamic Systems Part III: Nonlinear Systems: Extended Kalman Filter (EKF) and Unscented Kalman Filter (UKF) James B. Rawlings and Fernando V. Lima Department of
More informationarxiv: v1 [math.oc] 18 Sep 2007
Report no. 07/18 arxiv:0709.2791v1 math.oc 18 Sep 2007 Kalman Filtering with Equality and Inequality State Constraints Nachi Gupta Raphael Hauser Oxford University Computing Laboratory, Numerical Analysis
More informationState Estimation using Moving Horizon Estimation and Particle Filtering
State Estimation using Moving Horizon Estimation and Particle Filtering James B. Rawlings Department of Chemical and Biological Engineering UW Math Probability Seminar Spring 2009 Rawlings MHE & PF 1 /
More informationGaussian Process Approximations of Stochastic Differential Equations
Gaussian Process Approximations of Stochastic Differential Equations Cédric Archambeau Dan Cawford Manfred Opper John Shawe-Taylor May, 2006 1 Introduction Some of the most complex models routinely run
More informationEnsemble Kalman Filter
Ensemble Kalman Filter Geir Evensen and Laurent Bertino Hydro Research Centre, Bergen, Norway, Nansen Environmental and Remote Sensing Center, Bergen, Norway The Ensemble Kalman Filter (EnKF) Represents
More informationCamera Models and Affine Multiple Views Geometry
Camera Models and Affine Multiple Views Geometry Subhashis Banerjee Dept. Computer Science and Engineering IIT Delhi email: suban@cse.iitd.ac.in May 29, 2001 1 1 Camera Models A Camera transforms a 3D
More informationA Comparitive Study Of Kalman Filter, Extended Kalman Filter And Unscented Kalman Filter For Harmonic Analysis Of The Non-Stationary Signals
International Journal of Scientific & Engineering Research, Volume 3, Issue 7, July-2012 1 A Comparitive Study Of Kalman Filter, Extended Kalman Filter And Unscented Kalman Filter For Harmonic Analysis
More informationFundamentals of attitude Estimation
Fundamentals of attitude Estimation Prepared by A.Kaviyarasu Assistant Professor Department of Aerospace Engineering Madras Institute Of Technology Chromepet, Chennai Basically an IMU can used for two
More informationThe T2K experiment Results and Perspectives. PPC2017 Corpus Christi Mai 2017 Michel Gonin On behalf of the T2K collaboration
The T2K experiment Results and Perspectives PPC2017 Corpus Christi Mai 2017 Michel Gonin On behalf of the T2K collaboration 1 Overview Neutrino oscillations The T2K off-axis experiment Oscillation results
More informationLecture 4: Extended Kalman filter and Statistically Linearized Filter
Lecture 4: Extended Kalman filter and Statistically Linearized Filter Department of Biomedical Engineering and Computational Science Aalto University February 17, 2011 Contents 1 Overview of EKF 2 Linear
More information4 Derivations of the Discrete-Time Kalman Filter
Technion Israel Institute of Technology, Department of Electrical Engineering Estimation and Identification in Dynamical Systems (048825) Lecture Notes, Fall 2009, Prof N Shimkin 4 Derivations of the Discrete-Time
More informationIntroduction to Unscented Kalman Filter
Introduction to Unscented Kalman Filter 1 Introdution In many scientific fields, we use certain models to describe the dynamics of system, such as mobile robot, vision tracking and so on. The word dynamics
More informationII) Experimental Design
SLAC Experimental Advisory Committee --- September 12 th, 1997 II) Experimental Design Theory and simulations Great promise of significant scientific and technological achievements! How to realize this
More informationNew Fast Kalman filter method
New Fast Kalman filter method Hojat Ghorbanidehno, Hee Sun Lee 1. Introduction Data assimilation methods combine dynamical models of a system with typically noisy observations to obtain estimates of the
More informationElectromagnetic characterization of big aperture magnet used in particle beam cancer therapy
Electromagnetic characterization of big aperture magnet used in particle beam cancer therapy Jhonnatan Osorio Moreno M.Pullia, C.Priano Presented at Comsol conference 2012 Milan Milan 10 th October 2012
More informationPDE Model Reduction Using SVD
http://people.sc.fsu.edu/ jburkardt/presentations/fsu 2006.pdf... John 1 Max Gunzburger 2 Hyung-Chun Lee 3 1 School of Computational Science Florida State University 2 Department of Mathematics and School
More informationENGR352 Problem Set 02
engr352/engr352p02 September 13, 2018) ENGR352 Problem Set 02 Transfer function of an estimator 1. Using Eq. (1.1.4-27) from the text, find the correct value of r ss (the result given in the text is incorrect).
More informationFunctions and Equations
Canadian Mathematics Competition An activity of the Centre for Education in Mathematics and Computing, University of Waterloo, Waterloo, Ontario Euclid eworkshop # Functions and Equations c 006 CANADIAN
More informationSystem of Linear Equations. Slide for MA1203 Business Mathematics II Week 1 & 2
System of Linear Equations Slide for MA1203 Business Mathematics II Week 1 & 2 Function A manufacturer would like to know how his company s profit is related to its production level. How does one quantity
More informationTowards Proton Computed Tomography
SCIPP Towards Proton Computed Tomography L. R. Johnson, B. Keeney, G. Ross, H. F.-W. Sadrozinski, A. Seiden, D.C. Williams, L. Zhang Santa Cruz Institute for Particle Physics, UC Santa Cruz, CA 95064 V.
More informationDESIGNING A KALMAN FILTER WHEN NO NOISE COVARIANCE INFORMATION IS AVAILABLE. Robert Bos,1 Xavier Bombois Paul M. J. Van den Hof
DESIGNING A KALMAN FILTER WHEN NO NOISE COVARIANCE INFORMATION IS AVAILABLE Robert Bos,1 Xavier Bombois Paul M. J. Van den Hof Delft Center for Systems and Control, Delft University of Technology, Mekelweg
More informationLocal Positioning with Parallelepiped Moving Grid
Local Positioning with Parallelepiped Moving Grid, WPNC06 16.3.2006, niilo.sirola@tut.fi p. 1/?? TA M P E R E U N I V E R S I T Y O F T E C H N O L O G Y M a t h e m a t i c s Local Positioning with Parallelepiped
More informationVision 3D articielle Session 2: Essential and fundamental matrices, their computation, RANSAC algorithm
Vision 3D articielle Session 2: Essential and fundamental matrices, their computation, RANSAC algorithm Pascal Monasse monasse@imagine.enpc.fr IMAGINE, École des Ponts ParisTech Contents Some useful rules
More informationTechnical University of Denmark
Technical University of Denmark Page 1 of 10 pages Written test, 12 December 2012 Course name: Introduction to medical imaging Course no. 31540 Aids allowed: None. Pocket calculator not allowed "Weighting":
More informationInduced Planar Homologies in Epipolar Geometry
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 12, Number 4 (2016), pp. 3759 3773 Research India Publications http://www.ripublication.com/gjpam.htm Induced Planar Homologies in
More informationDATA MINING LECTURE 8. Dimensionality Reduction PCA -- SVD
DATA MINING LECTURE 8 Dimensionality Reduction PCA -- SVD The curse of dimensionality Real data usually have thousands, or millions of dimensions E.g., web documents, where the dimensionality is the vocabulary
More informationLinear Dynamical Systems (Kalman filter)
Linear Dynamical Systems (Kalman filter) (a) Overview of HMMs (b) From HMMs to Linear Dynamical Systems (LDS) 1 Markov Chains with Discrete Random Variables x 1 x 2 x 3 x T Let s assume we have discrete
More informationA Study of the Kalman Filter applied to Visual Tracking
A Study of the Kalman Filter applied to Visual Tracking Nathan Funk University of Alberta Project for CMPUT 652 December 7, 2003 Abstract This project analyzes the applicability of the Kalman filter as
More informationWhat is Image Deblurring?
What is Image Deblurring? When we use a camera, we want the recorded image to be a faithful representation of the scene that we see but every image is more or less blurry, depending on the circumstances.
More informationSIGNAL AND IMAGE RESTORATION: SOLVING
1 / 55 SIGNAL AND IMAGE RESTORATION: SOLVING ILL-POSED INVERSE PROBLEMS - ESTIMATING PARAMETERS Rosemary Renaut http://math.asu.edu/ rosie CORNELL MAY 10, 2013 2 / 55 Outline Background Parameter Estimation
More informationData Fusion Kalman Filtering Self Localization
Data Fusion Kalman Filtering Self Localization Armando Jorge Sousa http://www.fe.up.pt/asousa asousa@fe.up.pt Faculty of Engineering, University of Porto, Portugal Department of Electrical and Computer
More information1. A moving kaon decays into two pions, one of which is left at rest. (m K
Physics Qualifying Examination Part I September 12, 2015 7-Minute Questions 1. A moving kaon decays into two pions, one of which is left at rest. (m K 500 MeV, m π = 140 MeV). a. What is the total energy
More informationLecture 5: Control Over Lossy Networks
Lecture 5: Control Over Lossy Networks Yilin Mo July 2, 2015 1 Classical LQG Control The system: x k+1 = Ax k + Bu k + w k, y k = Cx k + v k x 0 N (0, Σ), w k N (0, Q), v k N (0, R). Information available
More informationStatus Report on the Survey and Alignment efforts at DESY
Status Report on the Survey and Alignment efforts at DESY Johannes J. Prenting Deutsches Elektronen-Synchrotron DESY, MEA2 IWAA 2008 KEK Japan Major New Projects @ DESY - : European X-ray Free Electron
More informationSequential Monte Carlo Methods for Bayesian Computation
Sequential Monte Carlo Methods for Bayesian Computation A. Doucet Kyoto Sept. 2012 A. Doucet (MLSS Sept. 2012) Sept. 2012 1 / 136 Motivating Example 1: Generic Bayesian Model Let X be a vector parameter
More informationA Theoretical Overview on Kalman Filtering
A Theoretical Overview on Kalman Filtering Constantinos Mavroeidis Vanier College Presented to professors: IVANOV T. IVAN STAHN CHRISTIAN Email: cmavroeidis@gmail.com June 6, 208 Abstract Kalman filtering
More informationExtended Kalman Filter based State Estimation of Wind Turbine
Extended Kalman Filter based State Estimation of Wind Turbine Kavitha N*, Vijayachitra S# *PG Scholar,.E. Control and Instrumentation Engineering #Professor, Department of Electronics and Instrumentation
More informationGraduate Accelerator Physics. G. A. Krafft Jefferson Lab Old Dominion University Lecture 1
Graduate Accelerator Physics G. A. Krafft Jefferson Lab Old Dominion University Lecture 1 Course Outline Course Content Introduction to Accelerators and Short Historical Overview Basic Units and Definitions
More informationExam in Automatic Control II Reglerteknik II 5hp (1RT495)
Exam in Automatic Control II Reglerteknik II 5hp (1RT495) Date: August 4, 018 Venue: Bergsbrunnagatan 15 sal Responsible teacher: Hans Rosth. Aiding material: Calculator, mathematical handbooks, textbooks
More informationState Estimation with a Kalman Filter
State Estimation with a Kalman Filter When I drive into a tunnel, my GPS continues to show me moving forward, even though it isn t getting any new position sensing data How does it work? A Kalman filter
More informationTracking for VR and AR
Tracking for VR and AR Hakan Bilen November 17, 2017 Computer Graphics University of Edinburgh Slide credits: Gordon Wetzstein and Steven M. La Valle 1 Overview VR and AR Inertial Sensors Gyroscopes Accelerometers
More informationKalman Filter Computer Vision (Kris Kitani) Carnegie Mellon University
Kalman Filter 16-385 Computer Vision (Kris Kitani) Carnegie Mellon University Examples up to now have been discrete (binary) random variables Kalman filtering can be seen as a special case of a temporal
More informationLecture 2: From Linear Regression to Kalman Filter and Beyond
Lecture 2: From Linear Regression to Kalman Filter and Beyond January 18, 2017 Contents 1 Batch and Recursive Estimation 2 Towards Bayesian Filtering 3 Kalman Filter and Bayesian Filtering and Smoothing
More informationThe Project. National Schools Observatory
Sunspots The Project This project is devised to give students a good understanding of the structure and magnetic field of the Sun and how this effects solar activity. Students will work with sunspot data
More informationEEE 187: Take Home Test #2
EEE 187: Take Home Test #2 Date: 11/30/2017 Due : 12/06/2017 at 5pm 1 Please read. Two versions of the exam are proposed. You need to solve one only. Version A: Four problems, Python is required for some
More informationMaximal Entropy for Reconstruction of Back Projection Images
Maximal Entropy for Reconstruction of Back Projection Images Tryphon Georgiou Department of Electrical and Computer Engineering University of Minnesota Minneapolis, MN 5545 Peter J Olver Department of
More informationImprovements of a Thermal Method for the Determination of Solar Absorptance and Thermal Emittance of Spacecraft Coatings
Improvements of a Thermal Method for the Determination of Solar Absorptance and Thermal Emittance of Spacecraft Coatings B. Biering, D. Hagelschuer German Aerospace Research Establishment Institute for
More informationFrom Bayes to Extended Kalman Filter
From Bayes to Extended Kalman Filter Michal Reinštein Czech Technical University in Prague Faculty of Electrical Engineering, Department of Cybernetics Center for Machine Perception http://cmp.felk.cvut.cz/
More informationGroßer Beleg. Real-Time Structure from Motion Using Kalman Filtering
Großer Beleg Real-Time Structure from Motion Using Kalman Filtering By Jeannette Bohg Born on 15th March 1981 in Cottbus Submitted on 21st March 2005 Overseeing Professor: Steffen Hölldobler Supervisor:
More informationInitial estimates for R 0 and Γ0 R Γ. Γ k. State estimate R Error covariance. k+1 k+1
HERA-B -32 VDET -1 CATS track fitting algorithm based on the discrete Kalman filter D. Emeliyanov 12 I. Kisel 34 1 Institute of Control Sciences, Profsoyuznaya ul., 65, Moscow, 11786 Russia 2 Deutsches
More informationAdaptive ensemble Kalman filtering of nonlinear systems
Adaptive ensemble Kalman filtering of nonlinear systems Tyrus Berry George Mason University June 12, 213 : Problem Setup We consider a system of the form: x k+1 = f (x k ) + ω k+1 ω N (, Q) y k+1 = h(x
More informationPhotometric Stereo: Three recent contributions. Dipartimento di Matematica, La Sapienza
Photometric Stereo: Three recent contributions Dipartimento di Matematica, La Sapienza Jean-Denis DUROU IRIT, Toulouse Jean-Denis DUROU (IRIT, Toulouse) 17 December 2013 1 / 32 Outline 1 Shape-from-X techniques
More informationA Study of Kruppa s Equation for Camera Self-calibration
Proceedings of the International Conference of Machine Vision and Machine Learning Prague, Czech Republic, August 14-15, 2014 Paper No. 57 A Study of Kruppa s Equation for Camera Self-calibration Luh Prapitasari,
More informationValidity of the Analysis of Radiochromic Film Using Matlab Code
Validity of the Analysis of Radiochromic Film Using Matlab Code Contact Samuel.millington@stfc.ac.uk S.J.Millington J.S.Green D.C.Carroll Introduction Radiochromic film is typically used in experiments
More informationA Study of Covariances within Basic and Extended Kalman Filters
A Study of Covariances within Basic and Extended Kalman Filters David Wheeler Kyle Ingersoll December 2, 2013 Abstract This paper explores the role of covariance in the context of Kalman filters. The underlying
More informationQuis custodiet ipsos custodes?
Quis custodiet ipsos custodes? James B. Rawlings, Megan Zagrobelny, Luo Ji Dept. of Chemical and Biological Engineering, Univ. of Wisconsin-Madison, WI, USA IFAC Conference on Nonlinear Model Predictive
More informationAircraft Turbofan Engine Health Estimation Using Constrained Kalman Filtering. Funded by the NASA Aviation Safety Program June 16, 2003
Aircraft Turbofan Engine Health Estimation Using Constrained Kalman Filtering Dan Simon Electrical Engineering Dept. Cleveland State University Cleveland, Ohio Donald L. Simon US Army Research Laboratory
More informationLecture 2 and 3: Controllability of DT-LTI systems
1 Lecture 2 and 3: Controllability of DT-LTI systems Spring 2013 - EE 194, Advanced Control (Prof Khan) January 23 (Wed) and 28 (Mon), 2013 I LTI SYSTEMS Recall that continuous-time LTI systems can be
More informationPARTICLE BEAMS, TOOLS FOR MODERN SCIENCE AND MEDICINE Hans-H. Braun, CERN
5 th Particle Physics Workshop National Centre for Physics Quaid-i-Azam University Campus, Islamabad PARTICLE BEAMS, TOOLS FOR MODERN SCIENCE AND MEDICINE Hans-H. Braun, CERN 2 nd Lecture Examples of Modern
More informationNonlinear State Estimation! Particle, Sigma-Points Filters!
Nonlinear State Estimation! Particle, Sigma-Points Filters! Robert Stengel! Optimal Control and Estimation, MAE 546! Princeton University, 2017!! Particle filter!! Sigma-Points Unscented Kalman ) filter!!
More information[variable] = units (or dimension) of variable.
Dimensional Analysis Zoe Wyatt wyatt.zoe@gmail.com with help from Emanuel Malek Understanding units usually makes physics much easier to understand. It also gives a good method of checking if an answer
More informationFitting the Bragg peak for accurate proton range determination
Fitting the Bragg peak for accurate proton range determination Koen Lambrechts July 10, 2015 Abstract This paper focusses on the uncertainties in proton range determination in the framework of optimizing
More informationA Practical Method for Decomposition of the Essential Matrix
Applied Mathematical Sciences, Vol. 8, 2014, no. 176, 8755-8770 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.410877 A Practical Method for Decomposition of the Essential Matrix Georgi
More informationA Crash Course on Kalman Filtering
A Crash Course on Kalman Filtering Dan Simon Cleveland State University Fall 2014 1 / 64 Outline Linear Systems Probability State Means and Covariances Least Squares Estimation The Kalman Filter Unknown
More information6.4 Kalman Filter Equations
6.4 Kalman Filter Equations 6.4.1 Recap: Auxiliary variables Recall the definition of the auxiliary random variables x p k) and x m k): Init: x m 0) := x0) S1: x p k) := Ak 1)x m k 1) +uk 1) +vk 1) S2:
More information