Inverse Problem in Quantitative Susceptibility Mapping From Theory to Application. Jae Kyu Choi.
|
|
- Shannon Greene
- 6 years ago
- Views:
Transcription
1 Inverse Problem in Quantitative Susceptibility Mapping From Theory to Application Jae Kyu Choi Department of CSE, Yonsei University Hangzhou Zhe Jiang University Jae Kyu Choi (Yonsei University) Inverse Problem in QSM / 23
2 1 Introduction 2 Theories in Inverse Problem 3 Applications 4 Conclusions Jae Kyu Choi (Yonsei University) Inverse Problem in QSM / 23
3 Introduction Jae Kyu Choi (Yonsei University) Inverse Problem in QSM / 23
4 Quantitative Susceptibility Mapping (QSM) Magnetic susceptibility χ: an intrinsic property of the material relating the magnetization M and the magnetic field H via M = χh. QSM: aims to visualize magnetic susceptibility χ from MR data. MR Magnitude Image Susceptibility Image Applications: Disease diagnosis (Parkinson s, Alzheimer, stroke, blood calcification, etc.) Jae Kyu Choi (Yonsei University) Inverse Problem in QSM / 23
5 Jae Kyu Choi (Yonsei University) Inverse Problem in QSM / 23
6 Inverse Problem of QSM Solve the deconvolution problem for χ: ˆR3 ψ(x) = pv d(x y)χ(y)dy d(x) = 2x2 3 x 2 1 x 2 2 4π x 5 ( ) 1 Ψ(ξ) = }{{} 3 ξ2 3 X (ξ) = D(ξ)X (ξ). ξ 2 }{{} F(ψ)(ξ) F(χ)(ξ) (pv: the principal value of the integral, F: Fourier transform) Data: relative difference field (RDF) ψ (noisy) Integral kernel d: singular (d(rx) = r 3 d(x) for r > 0). = Data ψ (Noisy) Unit dipole field d Susceptibility χ Jae Kyu Choi (Yonsei University) Inverse Problem in QSM / 23
7 Challenging Issue ˆ d(x y)χ(y)dy! ξ32 1 Ψ(ξ) = X (ξ) 3 ξ 2 {z } {z } {z } F(χ)(ξ) F(ψ)(ξ) ψ(x) = pv (IP-I) R3 (IP-F) D(ξ) ill-posed since D(ξ) = 0 in Γ0 = ξ R3 : ξ12 + ξ22 2ξ32 = 0, and this leads to the streaking artifacts. Data ψ Jae Kyu Choi (Yonsei University) Γ0 in Fourier domain Inverse Problem in QSM Reconstructed χ using (IP-F) / 23
8 Theoretical Results Existence and Uniqueness Result Characterization of Streaking Artifacts Tool Microlocal Analysis & Hörmander s Theory Both of the results are established for the first time. Jae Kyu Choi (Yonsei University) Inverse Problem in QSM / 23
9 Source of Error Propagation For a given measurement ψ E, we aim to obtain χ E using (IP-I) or (IP-F): ψ(x) = lim d(x y)χ(y)dy ε 0 ˆ x y >ε ( ) 1 Ψ(ξ) = D(ξ)X (ξ) = 3 ξ2 3 X (ξ) ξ 2 (IP-I) (IP-F). (D : space of distributions, E : space of compactly supported distributions, S : space of tempered distributions) If ψ E satisfies (IP-F) for some χ E, then ψ must lie in E := { u E : û(ξ)/p(ξ) is bounded near Γ 0 }. Here, P(ξ) is the polynomial defined as P(ξ) = 4π2 3 (ξ2 1 + ξ2 2 2ξ2 3 ). Jae Kyu Choi (Yonsei University) Inverse Problem in QSM / 23
10 Theorem (Existence and Uniqueness [J.K.Choi et al ]) If ψ E, we have the unique χ E satisfying (IP-F), and X = F(χ) can be represented as X (ξ) = 4π 2 ξ 2 Ψ(ξ) P(ξ) 9ξ3 4 (Proof follows from Paley-Wiener-Schwartz theorem.) if ξ Γ 0 Ψ ξ 3 (ξ) if ξ Γ 0 \ {0}. ( ) Reference χ ψ E Reconstructed χ using ( ) Jae Kyu Choi (Yonsei University) Inverse Problem in QSM / 23
11 Cause of Streaking Artifacts Streaking artifacts: closely related with the PDE ) 13 2 P(D)χ = ( + χ = ψ x 2 3 (IP-PDE) The solution χ D to (IP-PDE) is expressed as ˆ χ (x) = E ( ψ)(x) = E(x y) yψ(y)dy ψ E (S-PDE) R 3 where E(x) is the fundamental solution of P(D): 3 4π x E(x) = 2 3 2(x2 1 + x2 2 ) if 2(x x 2 2) < x otherwise. E(x) has the singular support along { x R 3 : 2(x x2 2 ) = x2 3}. Jae Kyu Choi (Yonsei University) Inverse Problem in QSM / 23
12 Theorem (Characterization of Artifacts [J.K.Choi et al ]) For ψ E, the wave front set of χ = χ satisfies WF(χ) \ WF(ψ) { (t P(ξ) + x, ξ) : ξ Γ 0 \ {0}, t 0, (x, ξ) WF(ψ) } Moreover, if (x, ξ) WF(χ) \ WF(ψ), then 1 ξ Γ 0 2 for any open interval (a, b) containing 0 such that { (x + t P(ξ), ξ) : t (a, b) } WF(ψ) =, we have { (x + t P(ξ), ξ) : t (a, b) } WF(χ). Simulated ψ E \ E Fundamental solution E χ = E ( ψ) Jae Kyu Choi (Yonsei University) Inverse Problem in QSM / 23
13 Our Proposed Approach Jae Kyu Choi (Yonsei University) Inverse Problem in QSM / 23
14 Description Want to develop a reconstruction method w/o regularization. = Artifacted Reconstructed Image χ + True Susceptibility χ Streaking Artifacts χv Based on the decomposition of a measured data = Measured data ψ Jae Kyu Choi (Yonsei University) + True data ψ = d χ Inverse Problem in QSM Violator ψv / 23
15 Assume that both ψ and ψ are supported in Ω. Since ψ = d χ, it follows that ψ E. However, ψ E \ E due to the non-vanishing compatibility violator ψ v := ψ ψ. This ψ v is the cause of artifacts in the reconstructed image χ: ˆ χ(x) = E(x y) yψ(y)dy. (S-PDE) R 3 Since χ is compactly supported, we have ˆ { χ(x) for x R 3 \ Ω, χ v(x) = E(x y) yψ v(y)dy = R 3 χ(x) χ (x) for x Ω. (1) In theory, it is possible to recover ψ v (with a special form of point sources) from the knowledge of χ = χ v in R 3 \ Ω. Jae Kyu Choi (Yonsei University) Inverse Problem in QSM / 23
16 Observation Assume that ψ E with supp(ψ) Ω can be decomposed into ψ = ψ + ψ v with ψ v being of the form ψ v(x) = N c j δ(x x j ) c j R, (2) j=1 where x 1,, x N supp(ψ). Then the source of the streaking artifacts ψ v E is determined uniquely from { P(D)χv = ψ v in R 3, with a given χ D defined as (1) in R 3 \ Ω. χ v = χ in R 3 \ Ω Once ψ v is obtained from ψ, then we are able to recover the streaking artifacts uniquely. Jae Kyu Choi (Yonsei University) Inverse Problem in QSM / 23
17 Issues in Zero-Cone Related Data Loss Artifacts All the previous methods of QSM rely on the direct inversion of in the frequency domain. Ψ(ξ) = D(ξ)X (ξ) χ can be obtained formally by the formula χ = F 1 (Ψ/D) with Ψ/D properly defined as a limit near Γ 0. If ψ v 0, then there exist ξ Γ 0 such that Ψ(ξ)/D(ξ) =. Hence, for given data available in discrete sense, we should discard the information of Ψ/D on Γ ε defined as { } Γ ε := ξ R 3 : dist(ξ, Γ 0) < ε/2 where 1 ε the distance between the centers of the neighboring voxels in frequency domain 2 dist(ξ, Γ 0): the distance between ξ and Γ 0. Jae Kyu Choi (Yonsei University) Inverse Problem in QSM / 23
18 Most previous methods use regularization with the fidelity of the following approximation: χ F 1( Ψ/D (1 η Γε ) ) (η Γε :the characteristic function of the set Γ ε). The information loss of (Ψ/D)η Γε causes unexpected artifacts in the reconstructed image, different from the streaking artifacts χ v generated by ψ v. In theory, it must be [ χ = lim F 1( Ψ/D(1 η Γε ) ) ] χ v. ε 0 However, in the practical computation with discrete images, it is inevitable to lose the information of Ψ/D near the zero cone Γ 0. Even though we choose ε > 0 so small that Γ ε consists of all voxels containing ξ Γ 0 only, this information loss causes the unexpected artifacts, which makes F 1( Ψ/D(1 η Γε ) ) ψ different from χ. Jae Kyu Choi (Yonsei University) Inverse Problem in QSM / 23
19 = χ = χ + χv + χ = F 1 Ψ/D(1 ηγ ) ε Jae Kyu Choi (Yonsei University) χv + F 1 Ψ/D(1 ηγ ) χv ε Inverse Problem in QSM χv / 23
20 The artifacts due to the missing data near Γ 0 are closely related to the singular fiber of the compatible part ψ (or Σ x(ψ ).) We present the key observation about the cause of the artifacts due to the singular fiber of the compatible part ψ. Observation Let ψ = ψ + ψ v be the decomposition of ψ with the non-vanishing ψ v. With the mixed use of discrete and continuous variables, let ε be the distance between the centers of the neighboring voxels in the frequency domain. Define χ ε by χ ε = F 1 (Ψ/D(1 η Γε )). (3) The artifacts of χ ε χ caused by the loss of Ψ/D on Γ ε can be explained by using the wave front set. More precisely, the wave front set of χ ε χ is explained by the set (in the continuous concept) { (t P(ξ) + x, ξ) : ξ Σx(ψ ) (Γ 0 \ {0}), t R }. Jae Kyu Choi (Yonsei University) Inverse Problem in QSM / 23
21 ψ = ψ + ψv ψ ψv χε = χ + χv + χ χ + χ χv Jae Kyu Choi (Yonsei University) Inverse Problem in QSM / 23
22 Conclusions Jae Kyu Choi (Yonsei University) Inverse Problem in QSM / 23
23 Conclusion We established the theoretical ground for the inverse problem of QSM; the compatibility condition of the data ψ E is founded so that the inverse problem can be solvable uniquely. Any data that violates the condition will cause streaking artifacts due to Γ 0, which can be analyzed using the wave front set. For the first time try to find directly the source of the streaking artifacts. In theory, it is possible to extract the source of the streaking artifacts using their characteristics (e.g., wave front set.) In addition, observed that the non-vanishing ψ v (with degeneracy of D on Γ 0 in the frequency domain) produces two types of streaking artifact in the reconstructed image: one caused by ψ v itself and another by the loss of information of X on Γ 0. Future works may consider Proving whether it is possible to pin-point the source of streaking artifacts ψ v when it is compactly supported in Ω. (Our presented case is based on the analyticity of E. However, for general ψ v E, we may not be able to adopt the analyticity of E any more.) Developing a robust numerical method to extract ψ v out of ψ. (Even in our case, the numerical implementation may be unstable because our theory is currently based on the analyticity of E.) Jae Kyu Choi (Yonsei University) Inverse Problem in QSM / 23
Oscillatory integrals
Chapter Oscillatory integrals. Fourier transform on S The Fourier transform is a fundamental tool in microlocal analysis and its application to the theory of PDEs and inverse problems. In this first section
More informationIntroduction to Microlocal Analysis
Introduction to Microlocal Analysis First lecture: Basics Dorothea Bahns (Göttingen) Third Summer School on Dynamical Approaches in Spectral Geometry Microlocal Methods in Global Analysis University of
More informationLECTURE 5: THE METHOD OF STATIONARY PHASE
LECTURE 5: THE METHOD OF STATIONARY PHASE Some notions.. A crash course on Fourier transform For j =,, n, j = x j. D j = i j. For any multi-index α = (α,, α n ) N n. α = α + + α n. α! = α! α n!. x α =
More informationWe denote the space of distributions on Ω by D ( Ω) 2.
Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study
More informationA Walking Tour of Microlocal Analysis
A Walking Tour of Microlocal Analysis Jeff Schonert August 10, 2006 Abstract We summarize some of the basic principles of microlocal analysis and their applications. After reviewing distributions, we then
More informationMicro-local analysis in Fourier Lebesgue and modulation spaces.
Micro-local analysis in Fourier Lebesgue and modulation spaces. Stevan Pilipović University of Novi Sad Nagoya, September 30, 2009 (Novi Sad) Nagoya, September 30, 2009 1 / 52 Introduction We introduce
More information2 Infinite products and existence of compactly supported φ
415 Wavelets 1 Infinite products and existence of compactly supported φ Infinite products.1 Infinite products a n need to be defined via limits. But we do not simply say that a n = lim a n N whenever the
More informationLECTURES ON MICROLOCAL CHARACTERIZATIONS IN LIMITED-ANGLE
LECTURES ON MICROLOCAL CHARACTERIZATIONS IN LIMITED-ANGLE TOMOGRAPHY Jürgen Frikel 4 LECTURES 1 Nov. 11: Introduction to the mathematics of computerized tomography 2 Today: Introduction to the basic concepts
More informationMicrolocal Analysis : a short introduction
Microlocal Analysis : a short introduction Plamen Stefanov Purdue University Mini Course, Fields Institute, 2012 Plamen Stefanov (Purdue University ) Microlocal Analysis : a short introduction 1 / 25 Introduction
More information44 CHAPTER 3. CAVITY SCATTERING
44 CHAPTER 3. CAVITY SCATTERING For the TE polarization case, the incident wave has the electric field parallel to the x 3 -axis, which is also the infinite axis of the aperture. Since both the incident
More informationMICROLOCAL ANALYSIS METHODS
MICROLOCAL ANALYSIS METHODS PLAMEN STEFANOV One of the fundamental ideas of classical analysis is a thorough study of functions near a point, i.e., locally. Microlocal analysis, loosely speaking, is analysis
More informationWeak universality of the KPZ equation
Weak universality of the KPZ equation M. Hairer, joint with J. Quastel University of Warwick Buenos Aires, July 29, 2014 KPZ equation Introduced by Kardar, Parisi and Zhang in 1986. Stochastic partial
More informationWavelets and applications
Chapter 3 Wavelets and applications 3. Multiresolution analysis 3.. The limits of Fourier analysis Although along this chapter the underlying Hilbert space will be L 2 (R), we start with a completely explicit
More informationISSN: [Engida* et al., 6(4): April, 2017] Impact Factor: 4.116
IJESRT INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH TECHNOLOGY ELLIPTICITY, HYPOELLIPTICITY AND PARTIAL HYPOELLIPTICITY OF DIFFERENTIAL OPERATORS Temesgen Engida*, Dr.Vishwajeet S. Goswami
More informationAn inverse scattering problem in random media
An inverse scattering problem in random media Pedro Caro Joint work with: Tapio Helin & Matti Lassas Computational and Analytic Problems in Spectral Theory June 8, 2016 Outline Introduction and motivation
More informationRegularity Structures
Regularity Structures Martin Hairer University of Warwick Seoul, August 14, 2014 What are regularity structures? Algebraic structures providing skeleton for analytical models mimicking properties of Taylor
More informationWe have to prove now that (3.38) defines an orthonormal wavelet. It belongs to W 0 by Lemma and (3.55) with j = 1. We can write any f W 1 as
88 CHAPTER 3. WAVELETS AND APPLICATIONS We have to prove now that (3.38) defines an orthonormal wavelet. It belongs to W 0 by Lemma 3..7 and (3.55) with j =. We can write any f W as (3.58) f(ξ) = p(2ξ)ν(2ξ)
More informationImage amplitudes in reverse time migration/inversion
Image amplitudes in reverse time migration/inversion Chris Stolk 1, Tim Op t Root 2, Maarten de Hoop 3 1 University of Amsterdam 2 University of Twente 3 Purdue University TRIP Seminar, October 6, 2010
More information8 Singular Integral Operators and L p -Regularity Theory
8 Singular Integral Operators and L p -Regularity Theory 8. Motivation See hand-written notes! 8.2 Mikhlin Multiplier Theorem Recall that the Fourier transformation F and the inverse Fourier transformation
More informationThe wave front set of a distribution
The wave front set of a distribution The Fourier transform of a smooth compactly supported function u(x) decays faster than any negative power of the dual variable ξ; that is for every number N there exists
More informationMicrolocal Methods in X-ray Tomography
Microlocal Methods in X-ray Tomography Plamen Stefanov Purdue University Lecture I: Euclidean X-ray tomography Mini Course, Fields Institute, 2012 Plamen Stefanov (Purdue University ) Microlocal Methods
More informationEXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018
EXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018 While these notes are under construction, I expect there will be many typos. The main reference for this is volume 1 of Hörmander, The analysis of liner
More informationFourier Transform & Sobolev Spaces
Fourier Transform & Sobolev Spaces Michael Reiter, Arthur Schuster Summer Term 2008 Abstract We introduce the concept of weak derivative that allows us to define new interesting Hilbert spaces the Sobolev
More informationGaussian beams and the propagation of singularities
Gaussian beams and the propagation of singularities Dominic Dold, Karen Habermann, Ellen Powell March 17, 2014 Contents 1 Introduction 3 2 Geometrical optics 5 2.1 Construction of the approximation...............................
More informationNOTES FOR CARDIFF LECTURES ON MICROLOCAL ANALYSIS
NOTES FOR CARDIFF LECTURES ON MICROLOCAL ANALYSIS JARED WUNSCH Note that these lectures overlap with Alex s to a degree, to ensure a smooth handoff between lecturers! Our notation is mostly, but not completely,
More informationFeshbach-Schur RG for the Anderson Model
Feshbach-Schur RG for the Anderson Model John Z. Imbrie University of Virginia Isaac Newton Institute October 26, 2018 Overview Consider the localization problem for the Anderson model of a quantum particle
More informationYAO LIU. f(x t) + f(x + t)
DUNKL WAVE EQUATION & REPRESENTATION THEORY OF SL() YAO LIU The classical wave equation u tt = u in n + 1 dimensions, and its various modifications, have been studied for centuries, and one of the most
More informationCan There Be a General Theory of Fourier Integral Operators?
Can There Be a General Theory of Fourier Integral Operators? Allan Greenleaf University of Rochester Conference on Inverse Problems in Honor of Gunther Uhlmann UC, Irvine June 21, 2012 How I started working
More informationTopics in Harmonic Analysis Lecture 6: Pseudodifferential calculus and almost orthogonality
Topics in Harmonic Analysis Lecture 6: Pseudodifferential calculus and almost orthogonality Po-Lam Yung The Chinese University of Hong Kong Introduction While multiplier operators are very useful in studying
More informationWEYL S LEMMA, ONE OF MANY. Daniel W. Stroock
WEYL S LEMMA, ONE OF MANY Daniel W Stroock Abstract This note is a brief, and somewhat biased, account of the evolution of what people working in PDE s call Weyl s Lemma about the regularity of solutions
More informationFourier analysis, measures, and distributions. Alan Haynes
Fourier analysis, measures, and distributions Alan Haynes 1 Mathematics of diffraction Physical diffraction As a physical phenomenon, diffraction refers to interference of waves passing through some medium
More informationQuantitative Susceptibility Mapping and Susceptibility Tensor Imaging. Magnetization and Susceptibility
Quantitative Susceptibility Mapping and Susceptibility Tensor Imaging 1, Chunlei Liu, Ph.D. 1 Brain Imaging and Analysis Center Department of Radiology Duke University, Durham, NC, USA 1 Magnetization
More informationFourier transform of tempered distributions
Fourier transform of tempered distributions 1 Test functions and distributions As we have seen before, many functions are not classical in the sense that they cannot be evaluated at any point. For eample,
More informationEE 364B Project Final Report: SCP Solver for a Nonconvex Quantitative Susceptibility Mapping Formulation
EE 364B Project Final Report: SCP Solver for a Nonconvex Quantitative Susceptibility Mapping Formulation Grant Yang Thanchanok Teeraratkul June 4, 2014 1 Background Changes in iron concentration in brain
More informationELLIPTIC REGULARITY AND SOLVABILITY FOR PARTIAL DIFFERENTIAL EQUATIONS WITH COLOMBEAU COEFFICIENTS
Electronic Journal of Differential Equations, Vol. 2004(2004), No. 14, pp. 1 30. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) ELLIPTIC
More informationStability estimates in stationary inverse transport
Stability estimates in stationary inverse transport Guillaume Bal and Alexandre Jollivet April 8, 28 Abstract We study the stability of the reconstruction of the scattering and absorption coefficients
More informationSeismic inverse scattering by reverse time migration
Seismic inverse scattering by reverse time migration Chris Stolk 1, Tim Op t Root 2, Maarten de Hoop 3 1 University of Amsterdam 2 University of Twente 3 Purdue University MSRI Inverse Problems Seminar,
More informationBorn-Oppenheimer Approximation in a singular system Prague. O.T.Turgut Bogazici Univ.
Born-Oppenheimer Approximation in a singular system Prague O.T.Turgut Bogazici Univ. June 5, 206 How can we use Born-Oppenheimer in a singular system (if we can)? What is the usual approach: Let us consider
More informationAppliedMathematics. Iterative Methods for Photoacoustic Tomography with Variable Sound Speed. Markus Haltmeier, Linh V. Nguyen
Nr. 29 7. December 26 Leopold-Franzens-Universität Innsbruck Preprint-Series: Department of Mathematics - Applied Mathematics Iterative Methods for Photoacoustic Tomography with Variable Sound Speed Markus
More informationA new class of pseudodifferential operators with mixed homogenities
A new class of pseudodifferential operators with mixed homogenities Po-Lam Yung University of Oxford Jan 20, 2014 Introduction Given a smooth distribution of hyperplanes on R N (or more generally on a
More informationContributors and Resources
for Introduction to Numerical Methods for d-bar Problems Jennifer Mueller, presented by Andreas Stahel Colorado State University Bern University of Applied Sciences, Switzerland Lexington, Kentucky, May
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 15 Heat with a source So far we considered homogeneous wave and heat equations and the associated initial value problems on the whole line, as
More informationSOLUTIONS TO HOMEWORK ASSIGNMENT 4
SOLUTIONS TO HOMEWOK ASSIGNMENT 4 Exercise. A criterion for the image under the Hilbert transform to belong to L Let φ S be given. Show that Hφ L if and only if φx dx = 0. Solution: Suppose first that
More informationWhen is the single-scattering approximation valid? Allan Greenleaf
When is the single-scattering approximation valid? Allan Greenleaf University of Rochester, USA Mathematical and Computational Aspects of Radar Imaging ICERM October 17, 2017 Partially supported by DMS-1362271,
More informationQuantum Mechanics Solutions. λ i λ j v j v j v i v i.
Quantum Mechanics Solutions 1. (a) If H has an orthonormal basis consisting of the eigenvectors { v i } of A with eigenvalues λ i C, then A can be written in terms of its spectral decomposition as A =
More informationarxiv: v3 [math-ph] 8 Mar 2009
On determinism and well-posedness in multiple time dimensions arxiv:08.00v3 [math-ph] 8 Mar 009 By Walter Craig and Steven Weinstein Department of Mathematics and Statistics, McMaster University, Hamilton,
More informationAlberto Bressan. Department of Mathematics, Penn State University
Non-cooperative Differential Games A Homotopy Approach Alberto Bressan Department of Mathematics, Penn State University 1 Differential Games d dt x(t) = G(x(t), u 1(t), u 2 (t)), x(0) = y, u i (t) U i
More informationStrichartz Estimates for the Schrödinger Equation in Exterior Domains
Strichartz Estimates for the Schrödinger Equation in University of New Mexico May 14, 2010 Joint work with: Hart Smith (University of Washington) Christopher Sogge (Johns Hopkins University) The Schrödinger
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationThe Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition
The Dirichlet boundary problems for second order parabolic operators satisfying a Martin Dindos Sukjung Hwang University of Edinburgh Satellite Conference in Harmonic Analysis Chosun University, Gwangju,
More informationRESULTS ON FOURIER MULTIPLIERS
RESULTS ON FOURIER MULTIPLIERS ERIC THOMA Abstract. The problem of giving necessary and sufficient conditions for Fourier multipliers to be bounded on L p spaces does not have a satisfactory answer for
More informationUltra-Hyperbolic Equations
Ultra-Hyperbolic Equations Minh Truong, Ph.D. Fontbonne University December 16, 2016 MSI Tech Support (Institute) Slides - Beamer 09/2014 1 / 28 Abstract We obtained a classical solution in Fourier space
More informationDiffraction by Edges. András Vasy (with Richard Melrose and Jared Wunsch)
Diffraction by Edges András Vasy (with Richard Melrose and Jared Wunsch) Cambridge, July 2006 Consider the wave equation Pu = 0, Pu = D 2 t u gu, on manifolds with corners M; here g 0 the Laplacian, D
More informationEfficient Solution Methods for Inverse Problems with Application to Tomography Radon Transform and Friends
Efficient Solution Methods for Inverse Problems with Application to Tomography Radon Transform and Friends Alfred K. Louis Institut für Angewandte Mathematik Universität des Saarlandes 66041 Saarbrücken
More informationEulerian Vortex Motion in Two and Three Dimensions
Eulerian Vortex Motion in Two and Three Dimensions Igor Yanovsky June 2006 1 Objective An Eulerian approach is used to solve the motion of an incompressible fluid, in two and three dimensions, in which
More informationProceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005
Proceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005 SOME INVERSE SCATTERING PROBLEMS FOR TWO-DIMENSIONAL SCHRÖDINGER
More informationInverse Transport Theory of Photoacoustics
Inverse Transport Theory of Photoacoustics Guillaume Bal Department of Applied Physics & Applied Mathematics, Columbia University, New York, NY 17 E-mail: gb3@columbia.edu Alexandre Jollivet Department
More information1. Plurisubharmonic functions and currents The first part of the homework studies some properties of PSH functions.
MATH 263: PROBLEM SET 2: PSH FUNCTIONS, HORMANDER S ESTIMATES AND VANISHING THEOREMS 1. Plurisubharmonic functions and currents The first part of the homework studies some properties of PSH functions.
More informationInverse Theory Methods in Experimental Physics
Inverse Theory Methods in Experimental Physics Edward Sternin PHYS 5P10: 2018-02-26 Edward Sternin Inverse Theory Methods PHYS 5P10: 2018-02-26 1 / 29 1 Introduction Indirectly observed data A new research
More informationAN INVERSE PROBLEM FOR THE WAVE EQUATION WITH A TIME DEPENDENT COEFFICIENT
AN INVERSE PROBLEM FOR THE WAVE EQUATION WITH A TIME DEPENDENT COEFFICIENT Rakesh Department of Mathematics University of Delaware Newark, DE 19716 A.G.Ramm Department of Mathematics Kansas State University
More informationLECTURES ON MICROLOCAL CHARACTERIZATIONS IN LIMITED-ANGLE
LECTURES ON MICROLOCAL CHARACTERIZATIONS IN LIMITED-ANGLE TOMOGRAPHY Jürgen Frikel 4 LECTURES 1 Today: Introduction to the mathematics of computerized tomography 2 Mathematics of computerized tomography
More informationMagnetic wells in dimension three
Magnetic wells in dimension three Yuri A. Kordyukov joint with Bernard Helffer & Nicolas Raymond & San Vũ Ngọc Magnetic Fields and Semiclassical Analysis Rennes, May 21, 2015 Yuri A. Kordyukov (Ufa) Magnetic
More informationMicrolocal analysis and inverse problems Lecture 4 : Uniqueness results in admissible geometries
Microlocal analysis and inverse problems Lecture 4 : Uniqueness results in admissible geometries David Dos Santos Ferreira LAGA Université de Paris 13 Wednesday May 18 Instituto de Ciencias Matemáticas,
More informationTheory and applications of time-frequency analysis
Theory and applications of time-frequency analysis Ville Turunen (ville.turunen@aalto.fi) Abstract: When and how often something happens in a signal? By properly quantizing these questions, we obtain the
More informationMiroslav Engliš. (Ω) possesses a reproducing kernel the Bergman kernel B(x, y) of Ω; namely, B(, y) L 2 hol
BOUNDARY SINGULARITY OF POISSON AND HARMONIC BERGMAN KERNELS Miroslav Engliš Abstract. We give a complete description of the boundary behaviour of the Poisson kernel and the harmonic Bergman kernel of
More informationSEMICLASSICAL LAGRANGIAN DISTRIBUTIONS
SEMICLASSICAL LAGRANGIAN DISTRIBUTIONS SEMYON DYATLOV Abstract. These are (rather sloppy) notes for the talk given in UC Berkeley in March 2012, attempting to explain the global theory of semiclassical
More informationx λ ϕ(x)dx x λ ϕ(x)dx = xλ+1 λ + 1 ϕ(x) u = xλ+1 λ + 1 dv = ϕ (x)dx (x))dx = ϕ(x)
1. Distributions A distribution is a linear functional L : D R (or C) where D = C c (R n ) which is linear and continuous. 1.1. Topology on D. Let K R n be a compact set. Define D K := {f D : f outside
More informationChapter One. The Calderón-Zygmund Theory I: Ellipticity
Chapter One The Calderón-Zygmund Theory I: Ellipticity Our story begins with a classical situation: convolution with homogeneous, Calderón- Zygmund ( kernels on R n. Let S n 1 R n denote the unit sphere
More informationHyperbolic inverse problems and exact controllability
Hyperbolic inverse problems and exact controllability Lauri Oksanen University College London An inverse initial source problem Let M R n be a compact domain with smooth strictly convex boundary, and let
More informationLeast squares: introduction to the network adjustment
Least squares: introduction to the network adjustment Experimental evidence and consequences Observations of the same quantity that have been performed at the highest possible accuracy provide different
More informationUnbiased Risk Estimation as Parameter Choice Rule for Filter-based Regularization Methods
Unbiased Risk Estimation as Parameter Choice Rule for Filter-based Regularization Methods Frank Werner 1 Statistical Inverse Problems in Biophysics Group Max Planck Institute for Biophysical Chemistry,
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Error Estimates for Filtered Back Projection Matthias Beckmann and Armin Iske Nr. 2015-03 January 2015 Error Estimates for Filtered Back Projection Matthias
More informationHARMONIC ANALYSIS. Date:
HARMONIC ANALYSIS Contents. Introduction 2. Hardy-Littlewood maximal function 3. Approximation by convolution 4. Muckenhaupt weights 4.. Calderón-Zygmund decomposition 5. Fourier transform 6. BMO (bounded
More informationLecture 2: Some basic principles of the b-calculus
Lecture 2: Some basic principles of the b-calculus Daniel Grieser (Carl von Ossietzky Universität Oldenburg) September 20, 2012 Summer School Singular Analysis Daniel Grieser (Oldenburg) Lecture 2: Some
More informationS chauder Theory. x 2. = log( x 1 + x 2 ) + 1 ( x 1 + x 2 ) 2. ( 5) x 1 + x 2 x 1 + x 2. 2 = 2 x 1. x 1 x 2. 1 x 1.
Sep. 1 9 Intuitively, the solution u to the Poisson equation S chauder Theory u = f 1 should have better regularity than the right hand side f. In particular one expects u to be twice more differentiable
More informationInverse problems for hyperbolic PDEs
Inverse problems for hyperbolic PDEs Lauri Oksanen University College London Example: inverse problem for the wave equation Let c be a smooth function on Ω R n and consider the wave equation t 2 u c 2
More informationMATH 220 solution to homework 5
MATH 220 solution to homework 5 Problem. (i Define E(t = k(t + p(t = then E (t = 2 = 2 = 2 u t u tt + u x u xt dx u 2 t + u 2 x dx, u t u xx + u x u xt dx x [u tu x ] dx. Because f and g are compactly
More informationFree boundaries in fractional filtration equations
Free boundaries in fractional filtration equations Fernando Quirós Universidad Autónoma de Madrid Joint work with Arturo de Pablo, Ana Rodríguez and Juan Luis Vázquez International Conference on Free Boundary
More informationMathematical Introduction
Chapter 1 Mathematical Introduction HW #1: 164, 165, 166, 181, 182, 183, 1811, 1812, 114 11 Linear Vector Spaces: Basics 111 Field A collection F of elements a,b etc (also called numbers or scalars) with
More informationSingular integrals on NA groups
Singular integrals on NA groups joint work with Alessio Martini and Alessandro Ottazzi Maria Vallarino Dipartimento di Scienze Matematiche Politecnico di Torino LMS Midlands Regional Meeting and Workshop
More informationNotes. 1 Fourier transform and L p spaces. March 9, For a function in f L 1 (R n ) define the Fourier transform. ˆf(ξ) = f(x)e 2πi x,ξ dx.
Notes March 9, 27 1 Fourier transform and L p spaces For a function in f L 1 (R n ) define the Fourier transform ˆf(ξ) = f(x)e 2πi x,ξ dx. Properties R n 1. f g = ˆfĝ 2. δλ (f)(ξ) = ˆf(λξ), where δ λ f(x)
More informationQuantum Mechanics I Physics 5701
Quantum Mechanics I Physics 5701 Z. E. Meziani 02/23//2017 Physics 5701 Lecture Outline 1 General Formulation of Quantum Mechanics 2 Measurement of physical quantities and observables 3 Representations
More informationThe Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition
The Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition Sukjung Hwang CMAC, Yonsei University Collaboration with M. Dindos and M. Mitrea The 1st Meeting of
More informationHere we used the multiindex notation:
Mathematics Department Stanford University Math 51H Distributions Distributions first arose in solving partial differential equations by duality arguments; a later related benefit was that one could always
More informationPart 2 Introduction to Microlocal Analysis
Part 2 Introduction to Microlocal Analysis Birsen Yazıcı & Venky Krishnan Rensselaer Polytechnic Institute Electrical, Computer and Systems Engineering August 2 nd, 2010 Outline PART II Pseudodifferential
More informationHarmonic Analysis Homework 5
Harmonic Analysis Homework 5 Bruno Poggi Department of Mathematics, University of Minnesota November 4, 6 Notation Throughout, B, r is the ball of radius r with center in the understood metric space usually
More informationQuintic deficient spline wavelets
Quintic deficient spline wavelets F. Bastin and P. Laubin January 19, 4 Abstract We show explicitely how to construct scaling functions and wavelets which are quintic deficient splines with compact support
More informationAn Empirical Chaos Expansion Method for Uncertainty Quantification
An Empirical Chaos Expansion Method for Uncertainty Quantification Melvin Leok and Gautam Wilkins Abstract. Uncertainty quantification seeks to provide a quantitative means to understand complex systems
More informationInverse problems in statistics
Inverse problems in statistics Laurent Cavalier (Université Aix-Marseille 1, France) YES, Eurandom, 10 October 2011 p. 1/27 Table of contents YES, Eurandom, 10 October 2011 p. 2/27 Table of contents 1)
More informationWell-Posedness and Adiabatic Limit for Quantum Zakharov System
Well-Posedness and Adiabatic Limit for Quantum Zakharov System Yung-Fu Fang (joint work with Tsai-Jung Chen, Jun-Ichi Segata, Hsi-Wei Shih, Kuan-Hsiang Wang, Tsung-fang Wu) Department of Mathematics National
More informationCONVEXIFICATION OF COEFFICIENT INVERSE PROBLEMS
CONVEXIFICATION OF COEFFICIENT INVERSE PROBLEMS Michael V. Klibanov Department of Mathematics and Statistics University of North Carolina at Charlotte Charlotte, NC 28223, USA THE MOST RECENT RESULT PDE-based
More informationf (t) K(t, u) d t. f (t) K 1 (t, u) d u. Integral Transform Inverse Fourier Transform
Integral Transforms Massoud Malek An integral transform maps an equation from its original domain into another domain, where it might be manipulated and solved much more easily than in the original domain.
More informationAlgebras of singular integral operators with kernels controlled by multiple norms
Algebras of singular integral operators with kernels controlled by multiple norms Alexander Nagel Conference in Harmonic Analysis in Honor of Michael Christ This is a report on joint work with Fulvio Ricci,
More informationSINGULARITY DETECTION AND PROCESSING WITH WAVELETS
SINGULARITY DETECTION AND PROCESSING WITH WAVELETS Stephane Mallat and Wen Liang Hwang Courant Institute of Mathematical Sciences New York University, New York, NY 10012 Technical Report March 1991 Abstract
More informationAn Adaptive Pseudo-Wavelet Approach for Solving Nonlinear Partial Differential Equations
An Adaptive Pseudo-Wavelet Approach for Solving Nonlinear Partial Differential Equations Gregory Beylkin and James M. Keiser Wavelet Analysis and Applications, v.6, 1997, Academic Press. Contents 1 Introduction
More informationYongdeok Kim and Seki Kim
J. Korean Math. Soc. 39 (00), No. 3, pp. 363 376 STABLE LOW ORDER NONCONFORMING QUADRILATERAL FINITE ELEMENTS FOR THE STOKES PROBLEM Yongdeok Kim and Seki Kim Abstract. Stability result is obtained for
More informationRANDOM PROPERTIES BENOIT PAUSADER
RANDOM PROPERTIES BENOIT PAUSADER. Quasilinear problems In general, one consider the following trichotomy for nonlinear PDEs: A semilinear problem is a problem where the highest-order terms appears linearly
More informationNotes on Quantum Mechanics
Notes on Quantum Mechanics Kevin S. Huang Contents 1 The Wave Function 1 1.1 The Schrodinger Equation............................ 1 1. Probability.................................... 1.3 Normalization...................................
More information2 Lebesgue integration
2 Lebesgue integration 1. Let (, A, µ) be a measure space. We will always assume that µ is complete, otherwise we first take its completion. The example to have in mind is the Lebesgue measure on R n,
More informationThe double layer potential
The double layer potential In this project, our goal is to explain how the Dirichlet problem for a linear elliptic partial differential equation can be converted into an integral equation by representing
More information