Mathematical analysis of ultrafast ultrasound imaging
|
|
- Jemimah Lester
- 5 years ago
- Views:
Transcription
1 Mathematical analysis of ultrafast ultrasound imaging Giovanni S Alberti Department of Mathematics, University of Genoa Joint work with H. Ammari (ETH), F. Romero (ETH) and T. Wintz (Sony). IPMS 18, May 1-5 Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May / 33
2 From ultrasonography to ultrafast ultrasonography Filippo, 17 Errico et. al, Nature 57, 99 5, 15 Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May 1-5 / 33
3 Conventional ultrasonography Conventional ultrasound imaging: focused ultrasonic waves high spatial resolution long acquisition time very low contrast: soft biological tissues are almost acoustically homogeneous, due to the high water concentration fine details (such as blood vessels) are completely invisible Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May / 33
4 Ultrafast ultrasound imaging Use of plane waves instead of focused waves High frame rate: up to, frames per second Lots of data to post-process Demené et al., IEEE Trans Med Imaging, 15. Errico et. al, Nature, 15 Single frame of Ultrafast ultrasound brain of a thinned skull rat Power doppler image obtained via a SVD filter applied to 5 frames. Superresolution: 75, frames with blinking microbubbles Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May 1-5 / 33
5 Summary 1 Mathematical modeling of Ultrafast Ultrasonography and blood flow imaging Super-resolution in Ultrafast Ultrasound Localization Microscopy Giovanni S. Alberti, Habib Ammari, Francisco Romero, and Timothée Wintz, Mathematical analysis of ultrafast ultrasound imaging, SIAM J. Appl. Math. 77 (17), no. 1, 1 5. Giovanni S. Alberti, Habib Ammari, Francisco Romero, and Timothée Wintz, Dynamic spike super-resolution and applications to ultrafast ultrasound imaging, arxiv preprint arxiv: (18). Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May / 33
6 Blood flow imaging The main issue is the removal of the clutter signal (the scattering coming from the tissue) Ultrafast ultrasonography allows us to overcome this issue, thanks to the very high frame rate. Idea: blood moves, tissue does not (in general). Temporal filters (Bercoff et al., 11): high-pass filtering the data to remove clutter signals. Drawback: not applicable when the clutter and blood velocities are close. Idea: tissue movement is spatially coherent, while blood flow is not. Spatiotemporal method based on the SVD of the data (Demene et al., 15): exploits the different spatial coherence of the clutter and blood scatterers. Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May / 33
7 The static direct problem y x Receptor array z Imaging plane The imaging system Time (s) 1-6 The pulse f(t) = e πiν t χ (ν t) Incident field in the direction k θ = (sin θ, cos θ): u i (x, y, z, t) = A z (y)f ( t c 1 k θ (x, z) ) c : background speed of sound. c(x): speed of sound. Perturbation: n(x) = 1 c (x) 1 c In the Born approximation, the scattered field takes the form: u s (u, t) = ˆR (π) 1 u x f ( t x k θ + u x c ) n (x ) dx, u = (u, ) Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May / 33
8 The static inverse problem: beamforming Scattered field u x x u k θ x x u s (u, t), u = (u, ), t > Travel time from the receptor array Γ to a point x and back to a receptor in u : τ θ x(u) = c 1 (k θ x + x u ) k θ Beamforming: averaging the signals z s θ (x, z) := ˆ x+f z x F z Inserting the expression for u s obtained before we obtain u s ( u, τ θ x (u) ) du ˆ ˆ x+f z s θ (x) = n (x ) (π) 1 x R x F z x u f ( τx θ (u) τx θ (u)) du dx }{{} =g θ (x,x ) Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May / 33
9 The static inverse problem: the point spread function The static image s θ may be rewritten as ˆ s θ (x) = g θ (x, x ) n (x ) dx, x R where g θ is the point spread function of the system: ˆ x+f z g θ (x, x (π) 1 ) = x F z x u f ( τx θ (u) τx θ (u)) du The PSF may be approximated with a convolution g θ (x, x ) g θ (x x ), s θ = g θ n, where (f = ν c 1 and χ = πiχ + χ ) g θ (x) iν F χ (f z) e πifz e πifθx sinc(πf F x) Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May / 33
10 The point spread function In the particular case θ = : g (x) iν F χ (f z) e πifz sinc(πf F x) The real part of the PSF g The real part of the PSF g (The size of the square is mm mm, and the horizontal and vertical axes are the x and z axes.) Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May / 33
11 Angle compounding In order to improve the decay in the x direction, (Montaldo et al., 9) introduced angle compounding: s ac Θ (x) = 1 ˆ Θ s θ (x) dθ, gθ ac (x) = 1 ˆ Θ g θ (x) dθ Θ Θ Θ Θ A simple derivation shows that the PSF is g ac Θ (x) = g (x) sinc(πν c 1 Θx) Θ = : we recover g θ for θ =. Θ > : this PSF enjoys faster decay in the variable x (a) g θ, θ = (b) g ac Θ, Θ =.5 Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May / 33
12 The dynamic forward problem The dynamic imaging setup consists in the repetition of the static imaging method over time to acquire a collection of images of a medium in motion. There are two time scales: the fast one related to the propagation of the wave is considered instantaneous with respect to the slow one, related to the sequence of the images. We now neglect the time of the propagation of a single wave to obtain static imaging. The time t considered here is related to the slow time scale. At fixed time t, we obtain a static image of the medium n = n(x, t): s(x, t) = (g ac Θ n(, t)) (x). Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May / 33
13 The dynamic inverse problem: Source separation Repeating the process for t [, T ] we obtain the movie s(x, t), which represents the main data we now need to process. Main aim: locating the (possibly very small) blood vessels. Main issue: s(x, t) is highly corrupted by clutter signal, namely, the signal scattered from tissues. Decompose The measurements are n(x, t) = n c (x, t) + n b (x, t) s(x, t) = s c (x, t) + s b (x, t) Inverse problem: determine the spatial support of n b. Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May / 33
14 The SVD separation algorithm (Demené et al., 15) (x i ) i=1,...,mx, (t j ) j=1,...,mt : sampling locations and times. Casorati matrix S C mx mt (m t m x ): The singular value decomposition of S S(i, j) = s(x i, t j ). m t S(i, j) = σ k u k (i) v k (j) k=1 singular vectors: (u1,..., u mx ) and (v 1,..., v mt ) are ONB of C mx and C m t singular values: σ1 σ... σ mt the dynamic data S are expressed as a sum of spatial components uk moving with time profiles v k, with weights σ k. Since the tissue movement has higher spatial coherence than the blood flow, the first factors are expected to contain the clutter signal, and the remainder to provide information about the blood location The blood location may be recovered by looking at the power Doppler m t k=k+1 σ k u k (i) Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May / 33
15 A general multiple scatterer random model Consider N point particles, with positions a k : i.i.d. stochastic processes a k (t), k = 1,..., N. The medium and the measurements are given by n (x, t) = C N N k=1 δ ak (t) (x), s(x, t) = C N N k=1 C > : scattering intensity N : natural normalization factor (central limit theorem) 1 Casorati matrix S N C mx mt : S N (i, j) = s(x i, t j ). g (x a k (t)) Multivariate central limit theorem: S N converges in distribution to a Gaussian complex matrix S C mx mt The distribution of S is entirely determined by g and the law of a k Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May / 33
16 Justification of the SVD method (1D) Using the multiple scatterer random model introduced above, we construct two Casorati matrices S b, S c as limits of particles with the following statistics. Clutter: large support, constant velocities a k (t) = u k + v c t where u k is uniformly distributed in (, L c ). Blood: small support, varying velocities: a k (t) = u k + v b t + σb t where u k is uniformly distributed in (, L b ) (L b L c ) and B t is a Brownian motion. S b and S c may be constructed using the Gaussian limit approximation Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May / 33
17 Justification of the SVD method (1D) # Clutter Blood Noise #1-1.5 Clutter (:5 cm.s!1 ) Clutter (1 cm.s!1 ) Clutter ( cm.s!1 ) #1 1 (a) The clutter model (v c = 1 m s 1 ), the blood model (σ = 1 6 m s 1, v b = 1 m s 1 ) and a white noise model with same variance as the blood #1 1 (b) The clutter model with different velocities. Figure: The distribution of the singular values of the Casorati matrix S in different cases. Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May / 33
18 z Numerical simulations We put one blood vessel in a moving tissue: domain: 5 mm 5 mm Θ = 7 The density of particles for both blood and clutter is, per mm C c = 5C b A single frame of the measurements s(x, t ) is # # x # Need further processing to locate the blood vessel! Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May / 33
19 Numerical simulations: v b > v c (a) Maximum blood velocity: cm s 1 ; mean clutter velocity: 1 cm s 1. 1 Clutter Blood Both Clutter Blood Both Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May 1-5 / 33
20 Numerical simulations: v b = v c Clutter Blood Both (b) Maximum blood velocity: 1 cm s 1 ; mean clutter velocity: 1 cm s 1. Clutter Blood Both Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May / 33
21 Numerical simulations: v b < v c (c) Maximum blood velocity:.5 cm s 1 ; mean clutter velocity: 1 cm s 1. 1 Clutter Blood Both Clutter Blood Both Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May 1-5 / 33
22 Super-resolution with ultrafast ultrasound Errico et. al, Nature 57, 99 5, 15 Use of randomly activated microbubbles as in Super-resolved Fluorescence Microscopy (Nobel Prize in Chemistry 1) Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May 1-5 / 33
23 Superresolution: static vs dynamic Current method a b: SVD b c: identify center of PSF, if well-separated Track bubbles to obtain velocities Drawbacks: slow discard a lot of data Errico et. al, Nature 57, 99 5, 15 New method dynamical superresolution in time and space based on l 1 minimization obtain locations and velocities in one step Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May / 33
24 The static super-resolution problem x 1,..., x N [, 1] d : locations of N particles w 1,..., w N C: weights The corresponding unknown Radon measure is t µ = N w i δ xi M([, 1] d ) i=1 For simplicity, we consider only low frequency Fourier measurements ˆ (Fµ) l = e πil x dµ(x), l { f c,..., f c } d, [,1] d where f c N is the highest measurable frequency. IP: reconstruct µ (i.e. the locations x i with infinite resolution) from (Fµ) l The unknown µ may be recovered via the convex problem min ν TV ν M([,1] d ) subject to Fν = Fµ, provided that the particles are well separated (Candes, Fernandez-Granda, CPAM 13) or d = 1 and w i (De Castro, Gamboa, JMAA, 1) Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May / 33
25 The dynamical super-resolution problem t t 1 t t1 t Time steps t k = kτ where τ > and k = K,..., K Dynamic spikes with positions x i (t k ) at time k: µ tk = N w i δ xi(t k ) i=1 For each t k, consider the static measurements Fµ tk, and so the full data are (Fµ tk ) l, l { f c,..., f c } d, k { K,..., K}. Aim: recover positions and velocities of the particles x i Static approach: at each time step t k, perform static reconstruction, and then use tracking to obtain the velocities We wish to design a fully dynamical reconstruction algorithm Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May / 33
26 Lifting the particles to the phase space Linear approximation of (local) movement of particles: µ tk = N w i δ xi+t k v i, where x i is the position at time t = and v i is the velocity. Lifting the unknown measure to N ω = w i δ (xi,v i) i=1 i=1 allows to separate the particles in the phase space: v 1 x 1 x v x 1 x Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May / 33
27 The measurements The measurements become (Gω) l,k = ˆ e πi(x+kτv) l dω(x, v) = ˆ e πi(x,v) (l,kτl) dω(x, v). In other words, we measure a d-dimensional Fourier transform restricted to { (l, kτl) : l Z d, l f c, k = K,..., K }, Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May / 33
28 The dynamical reconstruction Unknown measure: Measurements: N ω = w i δ (xi,v i) i=1 ˆ (Gω) l,k = e πi(x,v) (l,kτl) dω(x, v). t t 1 t t1 t We seek to reconstruct positions and velocities simultaneously by minimizing min λ TV subject to Gλ = Gω λ Theoretical results (see also (Dossal, Duval, Poon, SINUM 17)): exact reconstruction in the noiseless case, under suitable assumptions on (x i, v i) through a distance dyn in the phase space stable reconsctruction in presence of noise Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May / 33
29 Dynamical vs static reconstructions The results of 18, simulations in 1D with randomly generated particles are: 1. Correct recontruction rate dyn dynamic static static3 1/f c Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May / 33
30 An example We simulate blinking microbubbles in two close vessels with different velocities (a) frame (b) frame (c) frame 8 (d) B-mode image Super resolved reconstruction of the simulated particles The colors represent the velocity in the vertical direction Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May / 33
31 Conclusions I Mathematical modeling of ultrafast ultrasonography I Justification of the SVD algorithm by means of random particles I Design of a fully dynamical reconstruction method for super-resolution images of blood vessels Future directions I I I I Higher order approximation of the particles movement Efficient numerical implementation Simulations with real data Giovanni S Alberti (University of Genoa) Ultrafast ultrasound imaging IPMS 18, May / 33
arxiv: v2 [math.na] 1 Oct 2016
Mathematical Analysis of Ultrafast Ultrasound Imaging Giovanni S. Alberti Habib Ammari Francisco Romero Timothée Wintz September, arxiv:.v [math.na] Oct Abstract This paper provides a mathematical analysis
More informationCombining Sparsity with Physically-Meaningful Constraints in Sparse Parameter Estimation
UIUC CSL Mar. 24 Combining Sparsity with Physically-Meaningful Constraints in Sparse Parameter Estimation Yuejie Chi Department of ECE and BMI Ohio State University Joint work with Yuxin Chen (Stanford).
More informationSimulation of Contrast Agent Enhanced Ultrasound Imaging based on Field II
Simulation of Contrast Agent Enhanced Ultrasound Imaging based on Field II Tobias Gehrke, Heinrich M. Overhoff Medical Engineering Laboratory, University of Applied Sciences Gelsenkirchen tobias.gehrke@fh-gelsenkirchen.de
More informationINVERSION ASSUMING WEAK SCATTERING
INVERSION ASSUMING WEAK SCATTERING Angeliki Xenaki a,b, Peter Gerstoft b and Klaus Mosegaard a a Department of Applied Mathematics and Computer Science, Technical University of Denmark, 2800 Kgs. Lyngby,
More informationSparse Recovery Beyond Compressed Sensing
Sparse Recovery Beyond Compressed Sensing Carlos Fernandez-Granda www.cims.nyu.edu/~cfgranda Applied Math Colloquium, MIT 4/30/2018 Acknowledgements Project funded by NSF award DMS-1616340 Separable Nonlinear
More informationThe Physics of Doppler Ultrasound. HET408 Medical Imaging
The Physics of Doppler Ultrasound HET408 Medical Imaging 1 The Doppler Principle The basis of Doppler ultrasonography is the fact that reflected/scattered ultrasonic waves from a moving interface will
More informationWAVE PROPAGATION AND SCATTERING IN RANDOM MEDIA
WAVE PROPAGATION AND SCATTERING IN RANDOM MEDIA AKIRA ISHIMARU UNIVERSITY of WASHINGTON IEEE Antennas & Propagation Society, Sponsor IEEE PRESS The Institute of Electrical and Electronics Engineers, Inc.
More informationECE 8201: Low-dimensional Signal Models for High-dimensional Data Analysis
ECE 8201: Low-dimensional Signal Models for High-dimensional Data Analysis Lecture 7: Matrix completion Yuejie Chi The Ohio State University Page 1 Reference Guaranteed Minimum-Rank Solutions of Linear
More informationMicrowave-induced thermoacoustic tomography using multi-sector scanning
Microwave-induced thermoacoustic tomography using multi-sector scanning Minghua Xu, Geng Ku, and Lihong V. Wang a) Optical Imaging Laboratory, Biomedical Engineering Program, Texas A&M University, 3120
More informationFast and Robust Phase Retrieval
Fast and Robust Phase Retrieval Aditya Viswanathan aditya@math.msu.edu CCAM Lunch Seminar Purdue University April 18 2014 0 / 27 Joint work with Yang Wang Mark Iwen Research supported in part by National
More informationDOPPLER-BASED ULTRASONIC BLOOD VELOCITY ESTIMATION
DOPPLER-BASED ULTRASONIC BLOOD VELOCITY ESTIMATION BY TOROS ARIKAN THESIS Submitted in partial fulfillment of the requirements for the degree of Bachelor of Science in Electrical and Computer Engineering
More informationInverse Transport Problems and Applications. II. Optical Tomography and Clear Layers. Guillaume Bal
Inverse Transport Problems and Applications II. Optical Tomography and Clear Layers Guillaume Bal Department of Applied Physics & Applied Mathematics Columbia University http://www.columbia.edu/ gb23 gb23@columbia.edu
More informationA method of the forward problem for magneto-acousto-electrical tomography
Technology and Health Care 24 (216) S733 S738 DOI 1.3233/THC-16122 IOS Press S733 A method of the forward problem for magneto-acousto-electrical tomography Jingxiang Lv a,b,c, Guoqiang Liu a,, Xinli Wang
More informationOverview. Optimization-Based Data Analysis. Carlos Fernandez-Granda
Overview Optimization-Based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_spring16 Carlos Fernandez-Granda 1/25/2016 Sparsity Denoising Regression Inverse problems Low-rank models Matrix completion
More informationL-statistics based Modification of Reconstruction Algorithms for Compressive Sensing in the Presence of Impulse Noise
L-statistics based Modification of Reconstruction Algorithms for Compressive Sensing in the Presence of Impulse Noise Srdjan Stanković, Irena Orović and Moeness Amin 1 Abstract- A modification of standard
More informationConstructing Approximation Kernels for Non-Harmonic Fourier Data
Constructing Approximation Kernels for Non-Harmonic Fourier Data Aditya Viswanathan aditya.v@caltech.edu California Institute of Technology SIAM Annual Meeting 2013 July 10 2013 0 / 19 Joint work with
More informationImaging with Ambient Noise II
Imaging with Ambient Noise II George Papanicolaou Stanford University Michigan State University, Department of Mathematics Richard E. Phillips Lecture Series April 21, 2009 With J. Garnier, University
More informationClassical behavior of magnetic dipole vector. P. J. Grandinetti
Classical behavior of magnetic dipole vector Z μ Y X Z μ Y X Quantum behavior of magnetic dipole vector Random sample of spin 1/2 nuclei measure μ z μ z = + γ h/2 group μ z = γ h/2 group Quantum behavior
More informationAn example of Bayesian reasoning Consider the one-dimensional deconvolution problem with various degrees of prior information.
An example of Bayesian reasoning Consider the one-dimensional deconvolution problem with various degrees of prior information. Model: where g(t) = a(t s)f(s)ds + e(t), a(t) t = (rapidly). The problem,
More informationPixel-based Beamforming for Ultrasound Imaging
Pixel-based Beamforming for Ultrasound Imaging Richard W. Prager and Nghia Q. Nguyen Department of Engineering Outline v Introduction of Ultrasound Imaging v Image Formation and Beamforming v New Time-delay
More informationA STATE-SPACE APPROACH FOR THE ANALYSIS OF WAVE AND DIFFUSION FIELDS
ICASSP 2015 A STATE-SPACE APPROACH FOR THE ANALYSIS OF WAVE AND DIFFUSION FIELDS Stefano Maranò Donat Fäh Hans-Andrea Loeliger ETH Zurich, Swiss Seismological Service, 8092 Zürich ETH Zurich, Dept. Information
More informationROBUST BLIND SPIKES DECONVOLUTION. Yuejie Chi. Department of ECE and Department of BMI The Ohio State University, Columbus, Ohio 43210
ROBUST BLIND SPIKES DECONVOLUTION Yuejie Chi Department of ECE and Department of BMI The Ohio State University, Columbus, Ohio 4 ABSTRACT Blind spikes deconvolution, or blind super-resolution, deals with
More informationVirtual Array Processing for Active Radar and Sonar Sensing
SCHARF AND PEZESHKI: VIRTUAL ARRAY PROCESSING FOR ACTIVE SENSING Virtual Array Processing for Active Radar and Sonar Sensing Louis L. Scharf and Ali Pezeshki Abstract In this paper, we describe how an
More informationFast Angular Synchronization for Phase Retrieval via Incomplete Information
Fast Angular Synchronization for Phase Retrieval via Incomplete Information Aditya Viswanathan a and Mark Iwen b a Department of Mathematics, Michigan State University; b Department of Mathematics & Department
More informationCOMPRESSIVE SENSING WITH AN OVERCOMPLETE DICTIONARY FOR HIGH-RESOLUTION DFT ANALYSIS. Guglielmo Frigo and Claudio Narduzzi
COMPRESSIVE SESIG WITH A OVERCOMPLETE DICTIOARY FOR HIGH-RESOLUTIO DFT AALYSIS Guglielmo Frigo and Claudio arduzzi University of Padua Department of Information Engineering DEI via G. Gradenigo 6/b, I
More informationMEDT8007 Simulation methods in ultrasound imaging
1 MEDT8007 Simulation methods in ultrasound imaging Dynamic objects Lasse Løvstakken Department of Circulation and Medical Imaging Norwegian University of Science and Technology 2 Lecture outline Simple
More informationBasic System. Basic System. Sonosite 180. Acuson Sequoia. Echo occurs at t=2z/c where c is approximately 1500 m/s or 1.5 mm/µs
Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2007 Ultrasound Lecture Sonosite 80 From Suetens 2002 Acuson Sequoia Basic System Basic System Echo occurs at t=2/c where c is approximately
More informationDistributed Sensor Imaging
Distributed Sensor Imaging George Papanicolaou Department of Mathematics - Stanford University URL: http://georgep.stanford.edu In collaboration with: Gregoire Derveaux, INRIA, France Chrysoula Tsogka,
More information+ + ( + ) = Linear recurrent networks. Simpler, much more amenable to analytic treatment E.g. by choosing
Linear recurrent networks Simpler, much more amenable to analytic treatment E.g. by choosing + ( + ) = Firing rates can be negative Approximates dynamics around fixed point Approximation often reasonable
More informationRecurrent scattering and memory effect at the Anderson transition
Recurrent scattering and memory effect at the Anderson transition Alexandre Aubry Institut Langevin CNRS UMR 7587, ESPCI ParisTech Paris, France Collaborators: John Page, Laura Cobus (University of Manitoba,
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 informationA Non-sparse Tutorial on Sparse FFTs
A Non-sparse Tutorial on Sparse FFTs Mark Iwen Michigan State University April 8, 214 M.A. Iwen (MSU) Fast Sparse FFTs April 8, 214 1 / 34 Problem Setup Recover f : [, 2π] C consisting of k trigonometric
More informationSuper-Resolution Ultrasound Imaging of Vascular Structures with High Temporal Resolution
> REPLACE THIS LINE WITH YOUR PAPER IDENTIFICATION NUMBER (DOUBLE-CLICK HERE TO EDIT) < 1 Super-Resolution Ultrasound Imaging of Vascular Structures with High Temporal Resolution Avinoam Bar-Zion, Charles
More informationLinear Response and Onsager Reciprocal Relations
Linear Response and Onsager Reciprocal Relations Amir Bar January 1, 013 Based on Kittel, Elementary statistical physics, chapters 33-34; Kubo,Toda and Hashitsume, Statistical Physics II, chapter 1; and
More information31545 Medical Imaging systems
Simulation of ultrasound systems and non-linear imaging 545 Medical Imaging systems Lecture 9: Simulation of ultrasound systems and non-linear imaging Jørgen Arendt Jensen Department of Electrical Engineering
More informationSparse Sensing in Colocated MIMO Radar: A Matrix Completion Approach
Sparse Sensing in Colocated MIMO Radar: A Matrix Completion Approach Athina P. Petropulu Department of Electrical and Computer Engineering Rutgers, the State University of New Jersey Acknowledgments Shunqiao
More informationSingular value decomposition. If only the first p singular values are nonzero we write. U T o U p =0
Singular value decomposition If only the first p singular values are nonzero we write G =[U p U o ] " Sp 0 0 0 # [V p V o ] T U p represents the first p columns of U U o represents the last N-p columns
More informationIEEE Ultrasonic symposium 2002
IEEE Ultrasonic symposium 2002 Short Course 6: Flow Measurements Hans Torp Department of Circulation and Medical Imaging TU, orway Internet-site for short course: http://www.ifbt.ntnu.no/~hanst/flowmeas02/index.html
More informationStructure of Biological Materials
ELEC ENG 3BA3: Structure of Biological Materials Notes for Lecture #19 Monday, November 22, 2010 6.5 Nuclear medicine imaging Nuclear imaging produces images of the distribution of radiopharmaceuticals
More informationRobust Principal Component Analysis
ELE 538B: Mathematics of High-Dimensional Data Robust Principal Component Analysis Yuxin Chen Princeton University, Fall 2018 Disentangling sparse and low-rank matrices Suppose we are given a matrix M
More informationCoherent imaging without phases
Coherent imaging without phases Miguel Moscoso Joint work with Alexei Novikov Chrysoula Tsogka and George Papanicolaou Waves and Imaging in Random Media, September 2017 Outline 1 The phase retrieval problem
More informationBME I5000: Biomedical Imaging
BME I5000: Biomedical Imaging Lecture 9 Magnetic Resonance Imaging (imaging) Lucas C. Parra, parra@ccny.cuny.edu Blackboard: http://cityonline.ccny.cuny.edu/ 1 Schedule 1. Introduction, Spatial Resolution,
More informationNotes on the point spread function and resolution for projection lens/corf. 22 April 2009 Dr. Raymond Browning
Notes on the point spread function and resolution for projection lens/corf Abstract April 009 Dr. Raymond Browning R. Browning Consultants, 4 John Street, Shoreham, NY 786 Tel: (63) 8 348 This is a calculation
More informationAsymptotic Behavior of Waves in a Nonuniform Medium
Available at http://pvamuedu/aam Appl Appl Math ISSN: 1932-9466 Vol 12, Issue 1 June 217, pp 217 229 Applications Applied Mathematics: An International Journal AAM Asymptotic Behavior of Waves in a Nonuniform
More informationSuper-resolution by means of Beurling minimal extrapolation
Super-resolution by means of Beurling minimal extrapolation Norbert Wiener Center Department of Mathematics University of Maryland, College Park http://www.norbertwiener.umd.edu Acknowledgements ARO W911NF-16-1-0008
More informationCreative Commons: Attribution 3.0 Hong Kong License
Title Single-ensemble-based eigen-processing methods for color flow imaging-part I. the Hankel-SVD filter Author(s) Yu, ACH; Cobbold, RSC Citation Ieee Transactions On Ultrasonics, Ferroelectrics, And
More informationApproximate Methods for Time-Reversal Processing of Large Seismic Reflection Data Sets. University of California, LLNL
147th Annual Meeting: Acoustical Society of America New York City, May 25, 2004 Approximate Methods for Time-Reversal Processing of Large Seismic Reflection Data Sets Speaker: James G. Berryman University
More informationInverse problem and optimization
Inverse problem and optimization Laurent Condat, Nelly Pustelnik CNRS, Gipsa-lab CNRS, Laboratoire de Physique de l ENS de Lyon Decembre, 15th 2016 Inverse problem and optimization 2/36 Plan 1. Examples
More informationTwo-Dimensional Blood Flow Velocity Estimation Using Ultrasound Speckle Pattern Dependence on Scan Direction and A-Line Acquisition Velocity
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Biological Systems Engineering: Papers and Publications Biological Systems Engineering 5-2013 Two-Dimensional Blood Flow
More informationDetermination of thin elastic inclusions from boundary measurements.
Determination of thin elastic inclusions from boundary measurements. Elena Beretta in collaboration with E. Francini, S. Vessella, E. Kim and J. Lee September 7, 2010 E. Beretta (Università di Roma La
More informationDoppler echocardiography & Magnetic Resonance Imaging. Doppler echocardiography. History: - Langevin developed sonar.
1 Doppler echocardiography & Magnetic Resonance Imaging History: - Langevin developed sonar. - 1940s development of pulse-echo. - 1950s development of mode A and B. - 1957 development of continuous wave
More informationEE 5345 Biomedical Instrumentation Lecture 6: slides
EE 5345 Biomedical Instrumentation Lecture 6: slides 129-147 Carlos E. Davila, Electrical Engineering Dept. Southern Methodist University slides can be viewed at: http:// www.seas.smu.edu/~cd/ee5345.html
More informationSupermodular ordering of Poisson arrays
Supermodular ordering of Poisson arrays Bünyamin Kızıldemir Nicolas Privault Division of Mathematical Sciences School of Physical and Mathematical Sciences Nanyang Technological University 637371 Singapore
More informationFourier Approach to Wave Propagation
Phys 531 Lecture 15 13 October 005 Fourier Approach to Wave Propagation Last time, reviewed Fourier transform Write any function of space/time = sum of harmonic functions e i(k r ωt) Actual waves: harmonic
More informationA DIRECT IMAGING ALGORITHM FOR EXTENDED TARGETS
A DIRECT IMAGING ALGORITHM FOR EXTENDED TARGETS SONGMING HOU, KNUT SOLNA, AND HONGKAI ZHAO Abstract. We present a direct imaging algorithm for extended targets. The algorithm is based on a physical factorization
More informationRegular Variation and Extreme Events for Stochastic Processes
1 Regular Variation and Extreme Events for Stochastic Processes FILIP LINDSKOG Royal Institute of Technology, Stockholm 2005 based on joint work with Henrik Hult www.math.kth.se/ lindskog 2 Extremes for
More informationLecture 1: Pragmatic Introduction to Stochastic Differential Equations
Lecture 1: Pragmatic Introduction to Stochastic Differential Equations Simo Särkkä Aalto University, Finland (visiting at Oxford University, UK) November 13, 2013 Simo Särkkä (Aalto) Lecture 1: Pragmatic
More informationSparse Solutions of an Undetermined Linear System
1 Sparse Solutions of an Undetermined Linear System Maddullah Almerdasy New York University Tandon School of Engineering arxiv:1702.07096v1 [math.oc] 23 Feb 2017 Abstract This work proposes a research
More informationEL-GY 6813/BE-GY 6203 Medical Imaging, Fall 2016 Final Exam
EL-GY 6813/BE-GY 6203 Medical Imaging, Fall 2016 Final Exam (closed book, 1 sheets of notes double sided allowed, no calculator or other electronic devices allowed) 1. Ultrasound Physics (15 pt) A) (9
More informationWavelet Based Image Restoration Using Cross-Band Operators
1 Wavelet Based Image Restoration Using Cross-Band Operators Erez Cohen Electrical Engineering Department Technion - Israel Institute of Technology Supervised by Prof. Israel Cohen 2 Layout Introduction
More informationGaussian processes. Chuong B. Do (updated by Honglak Lee) November 22, 2008
Gaussian processes Chuong B Do (updated by Honglak Lee) November 22, 2008 Many of the classical machine learning algorithms that we talked about during the first half of this course fit the following pattern:
More informationHigh Frame Rate Vector Velocity Estimation using Plane Waves and Transverse Oscillation
Downloaded from orbit.dtu.dk on: Nov 1, 218 High Frame Rate Vector Velocity Estimation using Plane Waves and Transverse Oscillation Jensen, Jonas; Stuart, Matthias Bo; Jensen, Jørgen Arendt Published in:
More informationI. Rayleigh Scattering. EE Lecture 4. II. Dipole interpretation
I. Rayleigh Scattering 1. Rayleigh scattering 2. Dipole interpretation 3. Cross sections 4. Other approximations EE 816 - Lecture 4 Rayleigh scattering is an approximation used to predict scattering from
More informationFeasibility of non-linear simulation for Field II using an angular spectrum approach
Downloaded from orbit.dtu.dk on: Aug 22, 218 Feasibility of non-linear simulation for using an angular spectrum approach Du, Yigang; Jensen, Jørgen Arendt Published in: 28 IEEE Ultrasonics Symposium Link
More informationRecovery of Compactly Supported Functions from Spectrogram Measurements via Lifting
Recovery of Compactly Supported Functions from Spectrogram Measurements via Lifting Mark Iwen markiwen@math.msu.edu 2017 Friday, July 7 th, 2017 Joint work with... Sami Merhi (Michigan State University)
More informationGoing off the grid. Benjamin Recht Department of Computer Sciences University of Wisconsin-Madison
Going off the grid Benjamin Recht Department of Computer Sciences University of Wisconsin-Madison Joint work with Badri Bhaskar Parikshit Shah Gonnguo Tang We live in a continuous world... But we work
More informationPhoto-Acoustic imaging in layered media
Photo-Acoustic imaging in layered media Faouzi TRIKI Université Grenoble-Alpes, France (joint works with Kui Ren, Texas at Austin, USA.) 4 juillet 2016 Photo-Acoustic effect Photo-Acoustic effect (Bell,
More informationAn eigenvalue method using multiple frequency data for inverse scattering problems
An eigenvalue method using multiple frequency data for inverse scattering problems Jiguang Sun Abstract Dirichlet and transmission eigenvalues have important applications in qualitative methods in inverse
More informationStochastic Processes. M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno
Stochastic Processes M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno 1 Outline Stochastic (random) processes. Autocorrelation. Crosscorrelation. Spectral density function.
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 informationAn Overview of Sparsity with Applications to Compression, Restoration, and Inverse Problems
An Overview of Sparsity with Applications to Compression, Restoration, and Inverse Problems Justin Romberg Georgia Tech, School of ECE ENS Winter School January 9, 2012 Lyon, France Applied and Computational
More informationDeep convolutional framelets:
Deep convolutional framelets: application to diffuse optical tomography : learning based approach for inverse scattering problems Jaejun Yoo NAVER Clova ML OUTLINE (bottom up!) I. INTRODUCTION II. EXPERIMENTS
More informationTwo-View Segmentation of Dynamic Scenes from the Multibody Fundamental Matrix
Two-View Segmentation of Dynamic Scenes from the Multibody Fundamental Matrix René Vidal Stefano Soatto Shankar Sastry Department of EECS, UC Berkeley Department of Computer Sciences, UCLA 30 Cory Hall,
More informationSPARSE signal representations have gained popularity in recent
6958 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 57, NO. 10, OCTOBER 2011 Blind Compressed Sensing Sivan Gleichman and Yonina C. Eldar, Senior Member, IEEE Abstract The fundamental principle underlying
More informationSynchronization Transitions in Complex Networks
Synchronization Transitions in Complex Networks Y. Moreno 1,2,3 1 Institute for Biocomputation and Physics of Complex Systems (BIFI) University of Zaragoza, Zaragoza 50018, Spain 2 Department of Theoretical
More informationLight Propagation in Free Space
Intro Light Propagation in Free Space Helmholtz Equation 1-D Propagation Plane waves Plane wave propagation Light Propagation in Free Space 3-D Propagation Spherical Waves Huygen s Principle Each point
More informationStatistical signal processing
Statistical signal processing Short overview of the fundamentals Outline Random variables Random processes Stationarity Ergodicity Spectral analysis Random variable and processes Intuition: A random variable
More informationLecture Notes 5: Multiresolution Analysis
Optimization-based data analysis Fall 2017 Lecture Notes 5: Multiresolution Analysis 1 Frames A frame is a generalization of an orthonormal basis. The inner products between the vectors in a frame and
More informationPolarization analysis on three-component seismic data
Bollettino di Geofisica Teorica ed Applicata Vol., n. -, pp. 9-; Sep.-Dec. Polarization analysis on three-component seismic data P. Capizzi, L. De Luca and M. Vitale Department of Chemistry and Physics
More informationOn the recovery of measures without separation conditions
Norbert Wiener Center Department of Mathematics University of Maryland, College Park http://www.norbertwiener.umd.edu Applied and Computational Mathematics Seminar Georgia Institute of Technology October
More informationA Probability Review
A Probability Review Outline: A probability review Shorthand notation: RV stands for random variable EE 527, Detection and Estimation Theory, # 0b 1 A Probability Review Reading: Go over handouts 2 5 in
More informationAll-Optical Delay with Large Dynamic Range Using Atomic Dispersion
All-Optical Delay with Large Dynamic Range Using Atomic Dispersion M. R. Vanner, R. J. McLean, P. Hannaford and A. M. Akulshin Centre for Atom Optics and Ultrafast Spectroscopy February 2008 Motivation
More informationSupplement (videos)
Supplement (videos) Ruben s tube (sound): http://www.youtube.com/watch?v=gpcquuwqayw Doppler US (diagnostic use): http://www.youtube.com/watch?v=fgxzg-j_hfw http://www.youtube.com/watch?v=upsmenyoju8 High
More informationInfinite Series. 1 Introduction. 2 General discussion on convergence
Infinite Series 1 Introduction I will only cover a few topics in this lecture, choosing to discuss those which I have used over the years. The text covers substantially more material and is available for
More informationPredicting Phase Synchronization for Homoclinic Chaos in a CO 2 Laser
Predicting Phase Synchronization for Homoclinic Chaos in a CO 2 Laser Isao Tokuda, Jürgen Kurths, Enrico Allaria, Riccardo Meucci, Stefano Boccaletti and F. Tito Arecchi Nonlinear Dynamics, Institute of
More informationGaussian Processes. Le Song. Machine Learning II: Advanced Topics CSE 8803ML, Spring 2012
Gaussian Processes Le Song Machine Learning II: Advanced Topics CSE 8803ML, Spring 01 Pictorial view of embedding distribution Transform the entire distribution to expected features Feature space Feature
More informationFourier transforms and convolution
Fourier transforms and convolution (without the agonizing pain) CS/CME/BioE/Biophys/BMI 279 Oct. 26, 2017 Ron Dror 1 Outline Why do we care? Fourier transforms Writing functions as sums of sinusoids The
More informationRecovering overcomplete sparse representations from structured sensing
Recovering overcomplete sparse representations from structured sensing Deanna Needell Claremont McKenna College Feb. 2015 Support: Alfred P. Sloan Foundation and NSF CAREER #1348721. Joint work with Felix
More informationAdaptive interferometric imaging in clutter and optimal illumination
Adaptive interferometric imaging in clutter and optimal illumination Liliana Borcea and George Papanicolaou and Chrysoula Tsogka Computational and Applied Mathematics, MS 134, Rice University, 61 Main
More informationECE 636: Systems identification
ECE 636: Systems identification Lectures 3 4 Random variables/signals (continued) Random/stochastic vectors Random signals and linear systems Random signals in the frequency domain υ ε x S z + y Experimental
More informationPerformance Analysis for Strong Interference Remove of Fast Moving Target in Linear Array Antenna
Performance Analysis for Strong Interference Remove of Fast Moving Target in Linear Array Antenna Kwan Hyeong Lee Dept. Electriacal Electronic & Communicaton, Daejin University, 1007 Ho Guk ro, Pochen,Gyeonggi,
More informationBROADBAND MIMO SONAR SYSTEM: A THEORETICAL AND EXPERIMENTAL APPROACH
BROADBAND MIMO SONAR SYSTM: A THORTICAL AND XPRIMNTAL APPROACH Yan Pailhas a, Yvan Petillot a, Chris Capus a, Keith Brown a a Oceans Systems Lab., School of PS, Heriot Watt University, dinburgh, Scotland,
More informationELE 538B: Sparsity, Structure and Inference. Super-Resolution. Yuxin Chen Princeton University, Spring 2017
ELE 538B: Sparsity, Structure and Inference Super-Resolution Yuxin Chen Princeton University, Spring 2017 Outline Classical methods for parameter estimation Polynomial method: Prony s method Subspace method:
More informationPart 3 - Image Formation
Part 3 - Image Formation Three classes of scattering outcomes Types of electron microscopes Example SEM image: fly nose Example TEM image: muscle Skeletal muscle. Cell and Tissue Ultrastructure Mercer
More informationDS-GA 1002 Lecture notes 10 November 23, Linear models
DS-GA 2 Lecture notes November 23, 2 Linear functions Linear models A linear model encodes the assumption that two quantities are linearly related. Mathematically, this is characterized using linear functions.
More informationThe Linear Sampling Method and the MUSIC Algorithm
CODEN:LUTEDX/(TEAT-7089)/1-6/(2000) The Linear Sampling Method and the MUSIC Algorithm Margaret Cheney Department of Electroscience Electromagnetic Theory Lund Institute of Technology Sweden Margaret Cheney
More informationRecent progress on the explicit inversion of geodesic X-ray transforms
Recent progress on the explicit inversion of geodesic X-ray transforms François Monard Department of Mathematics, University of Washington. Geometric Analysis and PDE seminar University of Cambridge, May
More informationMassive MIMO: Signal Structure, Efficient Processing, and Open Problems II
Massive MIMO: Signal Structure, Efficient Processing, and Open Problems II Mahdi Barzegar Communications and Information Theory Group (CommIT) Technische Universität Berlin Heisenberg Communications and
More informationGauge optimization and duality
1 / 54 Gauge optimization and duality Junfeng Yang Department of Mathematics Nanjing University Joint with Shiqian Ma, CUHK September, 2015 2 / 54 Outline Introduction Duality Lagrange duality Fenchel
More information