Feature Reconstruction in Tomography
|
|
- Brittney Stewart
- 5 years ago
- Views:
Transcription
1 Feature Reconstruction in Tomography Alfred K. Louis Institut für Angewandte Mathematik Universität des Saarlandes Saarbrücken Wien, July 20, 2009
2 Images
3 Typical Procedure
4 Typical Procedure Image Reconstruction
5 Typical Procedure Image Reconstruction Image Enhancement
6 Typical Procedure Image Reconstruction Image Enhancement Disadvantage: The two operations are non adjusted and possibly too time-consuming
7 Requirements for Repeated Solution of Same Problem with Different Data
8 Requirements for Repeated Solution of Same Problem with Different Data Combination of reconstruction and image analysis in ONE step
9 Requirements for Repeated Solution of Same Problem with Different Data Combination of reconstruction and image analysis in ONE step Precompute solution operator for efficient evaluation of the same problem with different data sets
10 Requirements for Repeated Solution of Same Problem with Different Data Combination of reconstruction and image analysis in ONE step Precompute solution operator for efficient evaluation of the same problem with different data sets Use properties of the operator for fast algorithms
11 Content 1 Formulation of Problem 2 Inverse Problems and Regularization 3 Approximate Inverse for Combining Reconstruction and Analysis 4 Radon Transform and Edge Detection 5 Radon Transform and Diffusion 6 Radon Transform and Wavelet Decomposition 7 Nonlinear Problems 8 Future
12 Formulation of Problem Content 1 Formulation of Problem 2 Inverse Problems and Regularization 3 Approximate Inverse for Combining Reconstruction and Analysis 4 Radon Transform and Edge Detection 5 Radon Transform and Diffusion 6 Radon Transform and Wavelet Decomposition 7 Nonlinear Problems 8 Future
13 Formulation of Problem Target Combine the two steps in ONE algorithm Instead of Solve reconstruction part Af = g Evaluate the result Lf Compute directly LA g
14 Formulation of Problem Feature Determination Image Enhancement Image Registration Classification
15 Formulation of Problem Examples of Existing Methods A = R, L = I 1 Lambda CT 2D : Smith, Vainberg Lambda 3D: L-Maass, Katsevich Lambda with SPECT, ET: Quinto, Öktem Lambda with SONAR : L., Quinto
16 Formulation of Problem Image Analysis compute enhanced image Lf operate on Lf
17 Formulation of Problem Image Analysis compute enhanced image Lf operate on Lf Example for Calculation Edge Detection L kβ f = ( ) G β f x k
18 Formulation of Problem Image Analysis compute enhanced image Lf operate on Lf Example for Calculation Edge Detection L kβ f = ( ) G β f x k L. 96: Compute directly derivative of solution Schuster, 2000 : Compute divergence of vector fields
19 Formulation of Problem Applications for X-Ray CT Detect defects like blowholes In-line inspection 3D Dimensioning ( dimensionelles Messen )
20 Formulation of Problem X-Ray CT
21 Formulation of Problem X-Ray Tomography Assumptions: X-Rays travel on straight lines Attenuation proportional to path length Attenuation proportional to number of photons Radon Transform I = I t f g L = ln(i 0 /I L ) = f dt L
22 Formulation of Problem X-Ray Tomography Assumptions: X-Rays travel on straight lines Attenuation proportional to path length Attenuation proportional to number of photons Radon Transform I = I t f g L = ln(i 0 /I L ) = f dt L
23 Inverse Problems Content 1 Formulation of Problem 2 Inverse Problems and Regularization 3 Approximate Inverse for Combining Reconstruction and Analysis 4 Radon Transform and Edge Detection 5 Radon Transform and Diffusion 6 Radon Transform and Wavelet Decomposition 7 Nonlinear Problems 8 Future
24 Inverse Problems Inverse Problems OBSERVATION g SEARCHED-FOR DISTRIBUTION f Af = g A L(X, Y) Integral or Differential - Operator Difficulties: Solution does not exist Solution is not unique Solution does not depend continuously on g Problem ill - posed ( mal posé )
25 Inverse Problems Linear Regularization Y 1 A + Mγ A : X Y M γ X 1 A + S γ = M γ A + = A + Mγ
26 Inverse Problems General Purpose Regularization Methods R γ g = µ σ 1 µ g, u µ Y v µ
27 Inverse Problems General Purpose Regularization Methods R γ g = µ σ 1 µ F γ (σ µ ) g, u µ Y v µ
28 Inverse Problems General Purpose Regularization Methods R γ g = µ σ 1 µ F γ (σ µ ) g, u µ Y v µ Truncated SVD F γ (t) = { 1 : t γ 0 : t < γ
29 Inverse Problems General Purpose Regularization Methods R γ g = µ σ 1 µ F γ (σ µ ) g, u µ Y v µ Truncated SVD F γ (t) = { 1 : t γ 0 : t < γ Tikhonov - Phillips - Regularization ( Levenberg-Marquardt ) F γ (t) = t 2 /(t 2 + γ 2 )
30 Inverse Problems General Purpose Regularization Methods R γ g = µ σ 1 µ F γ (σ µ ) g, u µ Y v µ Truncated SVD F γ (t) = { 1 : t γ 0 : t < γ Tikhonov - Phillips - Regularization ( Levenberg-Marquardt ) F γ (t) = t 2 /(t 2 + γ 2 ) arg min{ Af g 2 Y + γ 2 f 2 X : f X }
31 Inverse Problems General Purpose Regularization Methods R γ g = µ σ 1 µ F γ (σ µ ) g, u µ Y v µ Truncated SVD F γ (t) = { 1 : t γ 0 : t < γ Tikhonov - Phillips - Regularization ( Levenberg-Marquardt ) F γ (t) = t 2 /(t 2 + γ 2 ) arg min{ Af g 2 Y + γ 2 f 2 X : f X } Iterative Methods (Landweber, CG, etc.)
32 Inverse Problems Regularization by Linear Functionals Likht, 1967 Backus-Gilbert, 1967 L.-Maass, 1990 L, 1996
33 Inverse Problems Regularization by Linear Functionals Likht, 1967 Backus-Gilbert, 1967 L.-Maass, 1990 L, 1996 Approximate Inverse: choose mollifier e γ (x, ) δ x and solve A ψ γ (x, ) = e γ (x, ) put f γ (x) := E γ f (x) := f, e γ (x, ) = g, ψ γ (x, )
34 Approximate Inverse for Combining Reconstruction and Analysis Content 1 Formulation of Problem 2 Inverse Problems and Regularization 3 Approximate Inverse for Combining Reconstruction and Analysis 4 Radon Transform and Edge Detection 5 Radon Transform and Diffusion 6 Radon Transform and Wavelet Decomposition 7 Nonlinear Problems 8 Future
35 Approximate Inverse for Combining Reconstruction and Analysis Approximate Inverse L.: SIAM J. Imaging Sciences , 2008 Aim Solve Af = g and compute Lf in one step Smoothed version (Lf ) γ = Lf, e γ with mollifier e γ
36 Approximate Inverse for Combining Reconstruction and Analysis Approximate Inverse L.: SIAM J. Imaging Sciences , 2008 Aim Solve Af = g and compute Lf in one step Smoothed version (Lf ) γ = Lf, e γ with mollifier e γ Assume a solution exists of A ψ γ (x, ) = L e γ (x, )
37 Approximate Inverse for Combining Reconstruction and Analysis Approximate Inverse L.: SIAM J. Imaging Sciences , 2008 Aim Solve Af = g and compute Lf in one step Smoothed version (Lf ) γ = Lf, e γ with mollifier e γ Assume a solution exists of A ψ γ (x, ) = L e γ (x, ) Then (Lf ) γ (x) = g, ψ γ (x, )
38 Approximate Inverse for Combining Reconstruction and Analysis Proof (Lf ) γ (x) = Lf, e γ (x, ) = f, L e γ (x, ) = f, A ψ γ (x, ) = g, ψ γ (x, )
39 Approximate Inverse for Combining Reconstruction and Analysis Proof (Lf ) γ (x) = Lf, e γ (x, ) = f, L e γ (x, ) = f, A ψ γ (x, ) = g, ψ γ (x, ) Conditions on e γ such that S γ is a regularization for determining Lf are given in the above mentioned paper.
40 Approximate Inverse for Combining Reconstruction and Analysis Example 1: L = I Consider with Then Mollifier A : L 2 (0, 1) L 2 (0, 1) Af (x) = x 0 f = g f (y)dy e γ (x, y) = 1 2γ χ [x γ,x+γ](y)
41 Approximate Inverse for Combining Reconstruction and Analysis Example 1 continued Auxiliary Equation: A ψ γ (x, y) = 1 y ψ γ (x, t)dt = e γ (x, y) leading to and ψ γ (x, y) = 1 2γ (δ x+γ δ x γ )(y) S γ g(x) = g(x + γ) g(x γ) 2γ
42 Approximate Inverse for Combining Reconstruction and Analysis Example 2: L = d dx Consider with Then Mollifier A : L 2 (0, 1) L 2 (0, 1) Af (x) = x 0 Lf = g f (y)dy e γ (x, y) = { y (x γ) γ x+γ y γ for x γ y x for x y x + γ
43 Approximate Inverse for Combining Reconstruction and Analysis Example 2 continued Auxiliary Equation A ψ γ (x, y) = leading to 1 y ψ γ (x, t)dt = L e γ (x, y) = y e γ(x, y) ψ γ (x, y) = 1 γ 2 (δ x+γ 2δ x + δ x γ )(y) and S γ g(x) = g(x + γ) 2g(x) + g(x γ) γ 2
44 Approximate Inverse for Combining Reconstruction and Analysis Existence Consider with and with Then formally and regularization via A : H s H s+α Af H s+α c 1 f H s L : D(L) L 2 L 2 Lf H s t c 2 f H s LA : L 2 H t α E γ : H t α L 2
45 Approximate Inverse for Combining Reconstruction and Analysis Possibilities
46 Approximate Inverse for Combining Reconstruction and Analysis Possibilities t < 0: addtional regularization needed ( Ex.: L differentiation )
47 Approximate Inverse for Combining Reconstruction and Analysis Possibilities t < 0: addtional regularization needed ( Ex.: L differentiation ) t = 0: only regularization of A needed ( Ex.: Wavelet decomposition )
48 Approximate Inverse for Combining Reconstruction and Analysis Possibilities t < 0: addtional regularization needed ( Ex.: L differentiation ) t = 0: only regularization of A needed ( Ex.: Wavelet decomposition ) t > α: no regularization at all ( Ex.: diffusion )
49 Approximate Inverse for Combining Reconstruction and Analysis Invariances Theorem (L, 1997/2007) Let the operator T 1, T 2, T 3 be such that L T 1 = T 2 L T 2 A = A T 3 and solve for a reference mollifier E γ the equation R Ψ γ = L E γ Then the general reconstruction kernel for the general mollifier e γ = T 1 E γ is ψ γ = T 3 Ψ γ and f γ = g, T 3 Ψ γ
50 Approximate Inverse for Combining Reconstruction and Analysis Possible Calculation of Reconstruction Kernel If the reconstruction of f is achieved as E γ f (x) = ψ γ (x, y)g(y)dy then the reconstruction kernel for calculating L β f can be computed as ψ βγ (, y) = L β ψ γ (, y)
51 Radon Transform and Edge Detection Content 1 Formulation of Problem 2 Inverse Problems and Regularization 3 Approximate Inverse for Combining Reconstruction and Analysis 4 Radon Transform and Edge Detection 5 Radon Transform and Diffusion 6 Radon Transform and Wavelet Decomposition 7 Nonlinear Problems 8 Future
52 Radon Transform and Edge Detection Radon Transform Rf (θ, s) = f (x)δ(s x θ)dx R 2 θ S 1 and s R. Rf (θ, σ) = (2π) 1/2ˆf (σθ) R 1 = 1 4π R I 1 where R is the adjoint operator from L 2 to L 2 known as backprojection R g(x) = g(θ, x θ)dθ S 1 and the Riesz potential I 1 is defined with the Fourier transform I 1 g(θ, σ) = σ ĝ(θ, σ)
53 Radon Transform and Edge Detection Radon Transform and Derivatives The Radon transform of a derivative is R f (θ, s) = θ k Rf (θ, s) x k s
54 Radon Transform and Edge Detection Invariances Let for x R 2 the shift operators T2 x f (y) = f (y x) and T x θ 3 g(θ, s) = g(θ, s x θ) then R θ T x 2 = T x θ 3 R θ Let U be a unitary 2 2 matrix and D2 U f (y) = f (Uy). then where D3 U g(θ, s) = g(uθ, s). RD U 2 = DU 3 R (T R) = R T
55 Radon Transform and Edge Detection Reconstruction Kernel Theorem Let the mollifier be given as e γ (x, y) = E γ ( x y ). Then the reconstruction kernel for L k = / x k is ψ k,γ = θ k Ψ γ ( s x θ ) with Ψ γ = 1 4π s I 1 RE γ
56 Radon Transform and Edge Detection Reconstruction Kernel Theorem Let the mollifier be given as e γ (x, y) = E γ ( x y ). Then the reconstruction kernel for L k = / x k is ψ k,γ = θ k Ψ γ ( s x θ ) with Ψ γ = 1 4π s I 1 RE γ Same costs for calculating f x k as for calculating f.
57 Radon Transform and Edge Detection Proof R ψ kγ = L k e γ = R 1 RL k e γ = 1 4π R I 1 RL k e γ Hence ψ kγ = 1 4π I 1 RL k e γ L k = L k and differentiation rule completes the proof.
58 Radon Transform and Edge Detection Example Let E βγ be given as Ê β γ(ξ) = (2π) 1 sinc ξ π β Then the reconstruction kernel for β = γ is sinc ξ π 2γ χ [ γ,γ]( ξ ) }{{} Shepp Logan kernel ψ π/h (lh) = 1 8l π 2 h 3 ( ) 3 + 4l 2 2, l Z 64l 2
59 Radon Transform and Edge Detection Kernel with β = γ and β = 2γ
60 Radon Transform and Edge Detection Exact Data Shepp - Logan phantom with original densities
61 Radon Transform and Edge Detection Noisy Data
62 Radon Transform and Edge Detection Differentiation of Image vs. New Method
63 Radon Transform and Edge Detection New Differentiation with respect to x and y.
64 Radon Transform and Edge Detection Lambda and wrong Parameter β
65 Radon Transform and Edge Detection Differentiation of Image vs. New Method Sum of absolute values of the two derivatives Apply support theorem of Boman and Quinto.
66 Radon Transform and Edge Detection Reconstruction from Measured Data, Fan Beam Geometry
67 Radon Transform and Edge Detection Applications 3D CT ( obvious with Feldkamp algorithm ) Phase Contrast Phase Contrast with Eikonal ( High Energy ) Approximation Electron Microscopy Adapt for cone beam tomography and inversion formula D 1 = cd TM Γ T where T is a differential and M Γ a trajectory dependent multiplication operator
68 Radon Transform and Diffusion Content 1 Formulation of Problem 2 Inverse Problems and Regularization 3 Approximate Inverse for Combining Reconstruction and Analysis 4 Radon Transform and Edge Detection 5 Radon Transform and Diffusion 6 Radon Transform and Wavelet Decomposition 7 Nonlinear Problems 8 Future
69 Radon Transform and Diffusion Diffusion Processes Smoothing of noisy images f by diffusion t u = u u(0, x) = f (x) Problem Edges are smeared out Perona-Malik: Nonlinear Diffusion u t = u u
70 Radon Transform and Diffusion Example Shepp-Logan phantom, 12% noise
71 Radon Transform and Diffusion Example Shepp-Logan phantom, 12% noise Normal reconstruction (left) and smoothing by linear diffusion, T = 0.1
72 Radon Transform and Diffusion Example Shepp-Logan phantom, 12% noise Normal reconstruction (left) and smoothing by linear diffusion, T = 0.1
73 Radon Transform and Diffusion Fundamental Solution L.: Inverse Problems and Imaging, to appear Fundamental Solution of the heat equation G(t, x) = 1 (4πt) n/2 exp( x 2 /4t) The solution of the heat equation for initial condition u(0, x) = f (x) can be calculated as L T u(0, ) = u(t, ) = G(T, ) f
74 Radon Transform and Diffusion Combining Reconstruction and Diffusion Let e γ (x, ) = δ x and L = L T Then the reconstruction kernel for time t = T has to fulfil and is of convolution type. R ψ T (x, y) = L T δ x(y) = G(T, x y)
75 Radon Transform and Diffusion Reconstruction Kernel Let Then and ψ T (s) = 1 2π 0 E γ (x) = δ x σ exp( T σ 2 ) cos(sσ)dσ ψ T (s) = 1 ( 4π 2 1 s D ( s T T 2 ) ) T with the Dawson integral s D(s) = exp( s 2 ) exp(t 2 )dt 0
76 Radon Transform and Diffusion Numerical Tests with noisy Data, 12 % noise
77 Radon Transform and Diffusion Numerical Tests with noisy Data, 12 % noise Smoothing by linear diffusion ( left ) and combined reconstruction (right).
78 Radon Transform and Diffusion Reconstruction of Derivatives for Noisy Data Reconstruction Kernel for Derivative in Direction x k ψ T,k (s) = θ ( k s 4π 2 T 3/2 2 T + ( 1 s2 2T ) D( s 2 T ) )
79 Radon Transform and Diffusion Noise and Differentiation, 6% Noise
80 Radon Transform and Diffusion Noise and Differentiation, 6% Noise Smoothing by linear diffusion ( left ) and combined reconstruction (right).
81 Radon Transform and Diffusion Noise and Differentiation, 6% and 12% Noise Combined reconstruction of sum of absolute values of derivatives for 6% and 12%noise.
82 Radon Transform and Diffusion Noise and Differentiation, 6% and 12% Noise Combined reconstruction of sum of absolute values of derivatives for 6% and 12%noise. No useful results with 12% noise with separate calculation of reconstruction and smoothed derivative.
83 Radon Transform and Wavelet Decomposition Content 1 Formulation of Problem 2 Inverse Problems and Regularization 3 Approximate Inverse for Combining Reconstruction and Analysis 4 Radon Transform and Edge Detection 5 Radon Transform and Diffusion 6 Radon Transform and Wavelet Decomposition 7 Nonlinear Problems 8 Future
84 Radon Transform and Wavelet Decomposition Wavelet Components Combination of reconstruction and wavelet decomposition: L, Maass, Rieder, 1994 Image reconstruction and wavelets : Holschneider, 1991 Berenstein, Walnut, 1996 Bhatia, Karl, Willsky, 1996 Bonnet, Peyrin, Turjman, 2002 Wang et al, 2004
85 Radon Transform and Wavelet Decomposition Inversion with Wavelets L.,MAASS,RIEDER 94 Represent the solution of Rf = g as f = k c M k ϕ Mk + m dl m ψ ml l
86 Radon Transform and Wavelet Decomposition Inversion with Wavelets L.,MAASS,RIEDER 94 Represent the solution of Rf = g as f = k c M k ϕ Mk + m dl m ψ ml l Precompute v Mk, w ml as R v Mk = ϕ Mk R w ml = ψ ml
87 Radon Transform and Wavelet Decomposition Inversion with Wavelets L.,MAASS,RIEDER 94 Represent the solution of Rf = g as f = k c M k ϕ Mk + m dl m ψ ml l Precompute v Mk, w ml as R v Mk = ϕ Mk Then R w ml = ψ ml c M k = g, v Mk d m l = g, w ml
88 Radon Transform and Wavelet Decomposition This is joint work with Steven Oeckl Department Process Integrated Inspection Systems Fraunhofer IIS
89 Reconstruction result: Three-dimensional registration of the object Main inspection task Dimensional measurement Detection of blow holes and porosity Feature Reconstruction in Tomography Radon Transform and Wavelet Decomposition New Application, first in NDT In-line inspection in production process
90 Radon Transform and Wavelet Decomposition Typical Scanning Geometry in NDT
91 Radon Transform and Wavelet Decomposition Wavelet Coefficient for 2D Parallel Geoemtry f, ψ jk = 1 I 1 g(θ, s)r θ ψ jk (s)dsdθ 4π S 1 IR = 1 ( I 1 ) g(θ, ) R θ ψ j0 ( D j k, θ )dθ 4π S 1
92 Radon Transform and Wavelet Decomposition 2D Shepp-Logan Phantom, Fan Beam
93 Radon Transform and Wavelet Decomposition Cone Beam: decentralized slice
94 Radon Transform and Wavelet Decomposition Cone Beam Local Reconstruction
95 Nonlinear Problems Content 1 Formulation of Problem 2 Inverse Problems and Regularization 3 Approximate Inverse for Combining Reconstruction and Analysis 4 Radon Transform and Edge Detection 5 Radon Transform and Diffusion 6 Radon Transform and Wavelet Decomposition 7 Nonlinear Problems 8 Future
96 Nonlinear Problems Considered Nonlinearity where Af = A l f l=1 A 1 f (x) = k 1 (x, y)f (y)dy A 2 f (x) = k 2 (x, y 1, y 2 )f (y 1 )f (y 2 )dy Considered by : Snieder for Backus-Gilbert Variants, 1991 L: Approximate Inverse, 1995
97 Nonlinear Problems Ansatz for Inversion, Presented for 2 Terms Put with and S γ g = S 1 g + S 2 g S 1 g(x) = g, ψ 1 (x, ) S 2 g(x) = g, Ψ 2 (x,, ), g Replace g by Af and omit higher order terms. Then S γ Af = S 1 A 1 f + S }{{} 1 A 2 f + S 2 A 1 f }{{} LE γf 0
98 Nonlinear Problems Determination of the Square Term Consider A : L 2 (Ω) R N Minimizing the defect leads to the equation N A 1 A 1 Ψ γ(x)a 1 A 1 = ψ γ,n (x)b n where B n = Ω Ω Note: B n is independent of x! n=1 k 1 (y 1 )k 2n (y 1, y 2 )k 1 (y 2 )dy 1 dy 2
99 Nonlinear Problems Example: A : L 2 R N Vibrating String u (x) + ρ(x)ω 2 u(x) = 0, u(0) = u(1) = 0 ρ(x) = 1 + f (x) Data g n = ω2 n (ω 0 n) 2 (ω 0 n) 2,, n = 1,..., N k 1,n (y) = 2 sin 2 (nπy) k 2,n (y 1, y 2 ) = 4 sin 2 (nπy 1 ) sin 2 (nπy 2 ) +4 n m n 2 n 2 m 2 sin(nπy 1) sin(mπy 1 ) sin(nπy 2 ) sin(mπy 2 )
100 Nonlinear Problems Numerical Test, L = d dx Linear and quadratic approximation of function (left) and derivative (right).
101 Future Content 1 Formulation of Problem 2 Inverse Problems and Regularization 3 Approximate Inverse for Combining Reconstruction and Analysis 4 Radon Transform and Edge Detection 5 Radon Transform and Diffusion 6 Radon Transform and Wavelet Decomposition 7 Nonlinear Problems 8 Future
102 Future Work in Progress Weakly nonlinear operators A Electron microscopy ( Kohr ) Inverse problems for Maxwell s equation ( A. Lakhal) Spherical Radon Transform ( Riplinger ) System Biology ( Groh ) Gait Analysis ( Bechtel, Johann ) Mathematical Finance ( M. Lakhal ) Semi - Discrete Problems ( Krebs ) 3D X-ray CT
103 Future Open Problems Nonlinear Problems, esp. nonlinear in L Time dependent problems other scanning geometries other operators L
Efficient 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 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 information1 Radon Transform and X-Ray CT
Radon Transform and X-Ray CT In this section we look at an important class of imaging problems, the inverse Radon transform. We have seen that X-ray CT (Computerized Tomography) involves the reconstruction
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 informationThe Radon transform. Chris Stolk. December 18, 2014
The Radon transform Chris Stolk December 18, 2014 1 Introduction In two dimensions the Radon transform is an integral transform that maps a function to its integrals over lines. Let θ S 1 and s R then
More informationIterative regularization of nonlinear ill-posed problems in Banach space
Iterative regularization of nonlinear ill-posed problems in Banach space Barbara Kaltenbacher, University of Klagenfurt joint work with Bernd Hofmann, Technical University of Chemnitz, Frank Schöpfer and
More information2D X-Ray Tomographic Reconstruction From Few Projections
2D X-Ray Tomographic Reconstruction From Few Projections Application of Compressed Sensing Theory CEA-LID, Thalès, UJF 6 octobre 2009 Introduction Plan 1 Introduction 2 Overview of Compressed Sensing Theory
More informationNumerical Methods for geodesic X-ray transforms and applications to open theoretical questions
Numerical Methods for geodesic X-ray transforms and applications to open theoretical questions François Monard Department of Mathematics, University of Washington. Nov. 13, 2014 UW Numerical Analysis Research
More informationSaturation Rates for Filtered Back Projection
Saturation Rates for Filtered Back Projection Matthias Beckmann, Armin Iske Fachbereich Mathematik, Universität Hamburg 3rd IM-Workshop on Applied Approximation, Signals and Images Bernried, 21st February
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 informationROI Reconstruction in CT
Reconstruction in CT From Tomo reconstruction in the 21st centery, IEEE Sig. Proc. Magazine (R.Clackdoyle M.Defrise) L. Desbat TIMC-IMAG September 10, 2013 L. Desbat Reconstruction in CT Outline 1 CT Radiology
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 informationIterative Methods for Inverse & Ill-posed Problems
Iterative Methods for Inverse & Ill-posed Problems Gerd Teschke Konrad Zuse Institute Research Group: Inverse Problems in Science and Technology http://www.zib.de/ag InverseProblems ESI, Wien, December
More informationLecture 9 February 2, 2016
MATH 262/CME 372: Applied Fourier Analysis and Winter 26 Elements of Modern Signal Processing Lecture 9 February 2, 26 Prof. Emmanuel Candes Scribe: Carlos A. Sing-Long, Edited by E. Bates Outline Agenda:
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 informationParameter Identification in Partial Differential Equations
Parameter Identification in Partial Differential Equations Differentiation of data Not strictly a parameter identification problem, but good motivation. Appears often as a subproblem. Given noisy observation
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 informationInverse problems in statistics
Inverse problems in statistics Laurent Cavalier (Université Aix-Marseille 1, France) Yale, May 2 2011 p. 1/35 Introduction There exist many fields where inverse problems appear Astronomy (Hubble satellite).
More informationEE 4372 Tomography. Carlos E. Davila, Dept. of Electrical Engineering Southern Methodist University
EE 4372 Tomography Carlos E. Davila, Dept. of Electrical Engineering Southern Methodist University EE 4372, SMU Department of Electrical Engineering 86 Tomography: Background 1-D Fourier Transform: F(
More informationAn Application of Abel s Method to the Inverse Radon Transform
J. Inverse Ill-Posed Probl.? (????), 1 14 DOI 1.1515/jip-????-??? de Gruyter???? An Application of Abel s Method to the Inverse Radon Transform Shavkat Alimov, Joseph David, Alexander Nolte and Julie Sherman
More informationApplication of local operators for numerical reconstruction of the singular support of a vector field by its known ray transforms
Application of local operators for numerical reconstruction of the singular support of a vector field by its known ray transforms E Yu Derevtsov 1, V V Pickalov 2 and T Schuster 3 1 Sobolev Institute of
More informationMumford-Shah Level-Set Ideas for Problems in Medical Imaging Paris Dauphine 2005
1/26 Mumford-Shah Level-Set Ideas for Problems in Medical Imaging Paris Dauphine 2005 Wolfgang Ring, Marc Droske, Ronny Ramlau Institut für Mathematik, Universität Graz General Mumford-Shah like 2/26 functionals
More informationRestoration of Missing Data in Limited Angle Tomography Based on Helgason- Ludwig Consistency Conditions
Restoration of Missing Data in Limited Angle Tomography Based on Helgason- Ludwig Consistency Conditions Yixing Huang, Oliver Taubmann, Xiaolin Huang, Guenter Lauritsch, Andreas Maier 26.01. 2017 Pattern
More informationIntro to harmonic analysis on groups Risi Kondor
Risi Kondor Any (sufficiently smooth) function f on the unit circle (equivalently, any 2π periodic f ) can be decomposed into a sum of sinusoidal waves f(x) = k= c n e ikx c n = 1 2π f(x) e ikx dx 2π 0
More informationIterative Methods for Ill-Posed Problems
Iterative Methods for Ill-Posed Problems Based on joint work with: Serena Morigi Fiorella Sgallari Andriy Shyshkov Salt Lake City, May, 2007 Outline: Inverse and ill-posed problems Tikhonov regularization
More informationRecovery of anisotropic metrics from travel times
Purdue University The Lens Rigidity and the Boundary Rigidity Problems Let M be a bounded domain with boundary. Let g be a Riemannian metric on M. Define the scattering relation σ and the length (travel
More informationFoundations of Image Science
Foundations of Image Science Harrison H. Barrett Kyle J. Myers 2004 by John Wiley & Sons,, Hoboken, 0-471-15300-1 1 VECTORS AND OPERATORS 1 1.1 LINEAR VECTOR SPACES 2 1.1.1 Vector addition and scalar multiplication
More informationRich Tomography. Bill Lionheart, School of Mathematics, University of Manchester and DTU Compute. July 2014
Rich Tomography Bill Lionheart, School of Mathematics, University of Manchester and DTU Compute July 2014 What do we mean by Rich Tomography? Conventional tomography reconstructs one scalar image from
More informationSparse tomography. Samuli Siltanen. Department of Mathematics and Statistics University of Helsinki, Finland
Sparse tomography Samuli Siltanen Department of Mathematics and Statistics University of Helsinki, Finland Minisymposium: Fourier analytic methods in tomographic image reconstruction Applied Inverse Problems
More informationQuadrature Formula for Computed Tomography
Quadrature Formula for Computed Tomography orislav ojanov, Guergana Petrova August 13, 009 Abstract We give a bivariate analog of the Micchelli-Rivlin quadrature for computing the integral of a function
More informationErkut Erdem. Hacettepe University February 24 th, Linear Diffusion 1. 2 Appendix - The Calculus of Variations 5.
LINEAR DIFFUSION Erkut Erdem Hacettepe University February 24 th, 2012 CONTENTS 1 Linear Diffusion 1 2 Appendix - The Calculus of Variations 5 References 6 1 LINEAR DIFFUSION The linear diffusion (heat)
More informationON A CLASS OF GENERALIZED RADON TRANSFORMS AND ITS APPLICATION IN IMAGING SCIENCE
ON A CLASS OF GENERALIZED RADON TRANSFORMS AND ITS APPLICATION IN IMAGING SCIENCE T.T. TRUONG 1 AND M.K. NGUYEN 2 1 University of Cergy-Pontoise, LPTM CNRS UMR 889, F-9532, France e-mail: truong@u-cergy.fr
More informationHigh-Resolution. Transmission. Electron Microscopy
Part 4 High-Resolution Transmission Electron Microscopy 186 Significance high-resolution transmission electron microscopy (HRTEM): resolve object details smaller than 1nm (10 9 m) image the interior of
More informationTowards Improved Sensitivity in Feature Extraction from Signals: one and two dimensional
Towards Improved Sensitivity in Feature Extraction from Signals: one and two dimensional Rosemary Renaut, Hongbin Guo, Jodi Mead and Wolfgang Stefan Supported by Arizona Alzheimer s Research Center and
More informationBIHARMONIC WAVE MAPS INTO SPHERES
BIHARMONIC WAVE MAPS INTO SPHERES SEBASTIAN HERR, TOBIAS LAMM, AND ROLAND SCHNAUBELT Abstract. A global weak solution of the biharmonic wave map equation in the energy space for spherical targets is constructed.
More informationThe Radon Transform and the Mathematics of Medical Imaging
Colby College Digital Commons @ Colby Honors Theses Student Research 2012 The Radon Transform and the Mathematics of Medical Imaging Jen Beatty Colby College Follow this and additional works at: http://digitalcommons.colby.edu/honorstheses
More informationInversions of ray transforms on simple surfaces
Inversions of ray transforms on simple surfaces François Monard Department of Mathematics, University of Washington. June 09, 2015 Institut Henri Poincaré - Program on Inverse problems 1 / 42 Outline 1
More informationNonlinear Flows for Displacement Correction and Applications in Tomography
Nonlinear Flows for Displacement Correction and Applications in Tomography Guozhi Dong 1 and Otmar Scherzer 1,2 1 Computational Science Center, University of Vienna, Oskar-Morgenstern-Platz 1, 1090 Wien,
More informationIntroduction of Partial Differential Equations and Boundary Value Problems
Introduction of Partial Differential Equations and Boundary Value Problems 2009 Outline Definition Classification Where PDEs come from? Well-posed problem, solutions Initial Conditions and Boundary Conditions
More informationLinear Inverse Problems
Linear Inverse Problems Ajinkya Kadu Utrecht University, The Netherlands February 26, 2018 Outline Introduction Least-squares Reconstruction Methods Examples Summary Introduction 2 What are inverse problems?
More informationNon-abelian Radon transform and its applications
Non-abelian Radon transform and its applications Roman Novikov CNRS, Centre de Mathématiques Appliquées, Ecole Polytechnique March 23, 2017 1/25 Consider the equation θ x ψ +A(x,θ)ψ = 0, x R d, θ S d 1,
More informationNIH Public Access Author Manuscript Inverse Probl. Author manuscript; available in PMC 2010 April 20.
NIH Public Access Author Manuscript Published in final edited form as: Inverse Probl. 2010 January 1; 26(3): 350131 3501329. doi:10.1088/0266-5611/26/3/035013. High Order Total Variation Minimization for
More informationFunctionalanalytic tools and nonlinear equations
Functionalanalytic tools and nonlinear equations Johann Baumeister Goethe University, Frankfurt, Germany Rio de Janeiro / October 2017 Outline Fréchet differentiability of the (PtS) mapping. Nonlinear
More informationLevenberg-Marquardt method in Banach spaces with general convex regularization terms
Levenberg-Marquardt method in Banach spaces with general convex regularization terms Qinian Jin Hongqi Yang Abstract We propose a Levenberg-Marquardt method with general uniformly convex regularization
More informationDecay profiles of a linear artificial viscosity system
Decay profiles of a linear artificial viscosity system Gene Wayne, Ryan Goh and Roland Welter Boston University rwelter@bu.edu July 2, 2018 This research was supported by funding from the NSF. Roland Welter
More informationInverse Transport Problems (Generalized) Stability Estimates. Guillaume Bal
Inverse Transport Problems (Generalized) Stability Estimates Guillaume Bal Department of Applied Physics & Applied Mathematics Columbia University Joint with Alexandre Jollivet Data Acquisition in CT-scan
More informationInterpolation procedure in filtered backprojection algorithm for the limited-angle tomography
Interpolation procedure in filtered backprojection algorithm for the limited-angle tomography Aleksander Denisiuk 1 Abstract A new data completion procedure for the limited-angle tomography is proposed.
More informationHigh-order total variation minimization for interior tomography
Home Search Collections Journals About Contact us My IOPscience High-order total variation minimization for interior tomography This article has been downloaded from IOPscience. Please scroll down to see
More informationBayesian X-ray Computed Tomography using a Three-level Hierarchical Prior Model
L. Wang, A. Mohammad-Djafari, N. Gac, MaxEnt 16, Ghent, Belgium. 1/26 Bayesian X-ray Computed Tomography using a Three-level Hierarchical Prior Model Li Wang, Ali Mohammad-Djafari, Nicolas Gac Laboratoire
More informationRecovery of Sparse Signals from Noisy Measurements Using an l p -Regularized Least-Squares Algorithm
Recovery of Sparse Signals from Noisy Measurements Using an l p -Regularized Least-Squares Algorithm J. K. Pant, W.-S. Lu, and A. Antoniou University of Victoria August 25, 2011 Compressive Sensing 1 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 informationReconstructing conductivities with boundary corrected D-bar method
Reconstructing conductivities with boundary corrected D-bar method Janne Tamminen June 24, 2011 Short introduction to EIT The Boundary correction procedure The D-bar method Simulation of measurement data,
More informationAn Introduction to X-ray tomography and Radon Transforms
Proceedings of Symposia in Applied Mathematics An Introduction to X-ray tomography and Radon Transforms Eric Todd Quinto Abstract. This article provides an introduction to the mathematics behind X-ray
More informationA Double Regularization Approach for Inverse Problems with Noisy Data and Inexact Operator
A Double Regularization Approach for Inverse Problems with Noisy Data and Inexact Operator Ismael Rodrigo Bleyer Prof. Dr. Ronny Ramlau Johannes Kepler Universität - Linz Florianópolis - September, 2011.
More informationIntroduction to the Mathematics of Medical Imaging
Introduction to the Mathematics of Medical Imaging Second Edition Charles L. Epstein University of Pennsylvania Philadelphia, Pennsylvania EiaJTL Society for Industrial and Applied Mathematics Philadelphia
More informationp 1 ( Y p dp) 1/p ( X p dp) 1 1 p
Doob s inequality Let X(t) be a right continuous submartingale with respect to F(t), t 1 P(sup s t X(s) λ) 1 λ {sup s t X(s) λ} X + (t)dp 2 For 1 < p
More informationImage Reconstruction from Projection
Image Reconstruction from Projection Reconstruct an image from a series of projections X-ray computed tomography (CT) Computed tomography is a medical imaging method employing tomography where digital
More informationRegion-of-interest reconstructions from truncated 3D x-ray projections
Region-of-interest reconstructions from truncated 3D x-ray projections Robert Azencott 1, Bernhard G. Bodmann 1, Demetrio Labate 1, Anando Sen 2, Daniel Vera 1 June 13, 2014 Abstract This paper introduces
More informationA Double Regularization Approach for Inverse Problems with Noisy Data and Inexact Operator
A Double Regularization Approach for Inverse Problems with Noisy Data and Inexact Operator Ismael Rodrigo Bleyer Prof. Dr. Ronny Ramlau Johannes Kepler Universität - Linz Cambridge - July 28, 211. Doctoral
More informationImproved Radon Based Imaging using the Shearlet Transform
Improved Radon Based Imaging using the Shearlet Transform Glenn R. Easley a, Flavia Colonna b, Demetrio Labate c a System Planning Corporation, Arlington, Virginia b George Mason University, Fairfax, Virginia
More informationQualitative Properties of Numerical Approximations of the Heat Equation
Qualitative Properties of Numerical Approximations of the Heat Equation Liviu Ignat Universidad Autónoma de Madrid, Spain Santiago de Compostela, 21 July 2005 The Heat Equation { ut u = 0 x R, t > 0, u(0,
More informationSemi-discrete equations in Banach spaces: The approximate inverse approach
Semi-discrete equations in Banach spaces: The approximate inverse approach Andreas Rieder Thomas Schuster Saarbrücken Frank Schöpfer Oldenburg FAKULTÄT FÜR MATHEMATIK INSTITUT FÜR ANGEWANDTE UND NUMERISCHE
More informationStarting from Heat Equation
Department of Applied Mathematics National Chiao Tung University Hsin-Chu 30010, TAIWAN 20th August 2009 Analytical Theory of Heat The differential equations of the propagation of heat express the most
More informationLinear Diffusion and Image Processing. Outline
Outline Linear Diffusion and Image Processing Fourier Transform Convolution Image Restoration: Linear Filtering Diffusion Processes for Noise Filtering linear scale space theory Gauss-Laplace pyramid for
More informationInverse problems Total Variation Regularization Mark van Kraaij Casa seminar 23 May 2007 Technische Universiteit Eindh ove n University of Technology
Inverse problems Total Variation Regularization Mark van Kraaij Casa seminar 23 May 27 Introduction Fredholm first kind integral equation of convolution type in one space dimension: g(x) = 1 k(x x )f(x
More informationThe Levenberg-Marquardt Iteration for Numerical Inversion of the Power Density Operator
The Levenberg-Marquardt Iteration for Numerical Inversion of the Power Density Operator G. Bal (gb2030@columbia.edu) 1 W. Naetar (wolf.naetar@univie.ac.at) 2 O. Scherzer (otmar.scherzer@univie.ac.at) 2,3
More informationError Estimates and Convergence Rates for Filtered Back Projection Reconstructions
Universität Hamburg Dissertation Error Estimates and Convergence ates for Filtered Back Projection econstructions Dissertation zur Erlangung des Doktorgrades an der Fakultät für Mathematik, Informatik
More informationNONLINEAR DIFFUSION PDES
NONLINEAR DIFFUSION PDES Erkut Erdem Hacettepe University March 5 th, 0 CONTENTS Perona-Malik Type Nonlinear Diffusion Edge Enhancing Diffusion 5 References 7 PERONA-MALIK TYPE NONLINEAR DIFFUSION The
More informationworked out from first principles by parameterizing the path, etc. If however C is a A path C is a simple closed path if and only if the starting point
III.c Green s Theorem As mentioned repeatedly, if F is not a gradient field then F dr must be worked out from first principles by parameterizing the path, etc. If however is a simple closed path in the
More informationI. Introduction. A New Method for Inverting Integrals Attenuated Radon Transform (SPECT) D to N map for Moving Boundary Value Problems
I. Introduction A New Method for Inverting Integrals Attenuated Radon Transform (SPECT) D to N map for Moving Boundary Value Problems F(k) = T 0 e k2 t+ikl(t) f (t)dt, k C. Integrable Nonlinear PDEs in
More informationA COUNTEREXAMPLE TO AN ENDPOINT BILINEAR STRICHARTZ INEQUALITY TERENCE TAO. t L x (R R2 ) f L 2 x (R2 )
Electronic Journal of Differential Equations, Vol. 2006(2006), No. 5, pp. 6. ISSN: 072-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) A COUNTEREXAMPLE
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 informationTwo-parameter regularization method for determining the heat source
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 13, Number 8 (017), pp. 3937-3950 Research India Publications http://www.ripublication.com Two-parameter regularization method for
More information1 Assignment 1: Nonlinear dynamics (due September
Assignment : Nonlinear dynamics (due September 4, 28). Consider the ordinary differential equation du/dt = cos(u). Sketch the equilibria and indicate by arrows the increase or decrease of the solutions.
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 informationBessel Functions Michael Taylor. Lecture Notes for Math 524
Bessel Functions Michael Taylor Lecture Notes for Math 54 Contents 1. Introduction. Conversion to first order systems 3. The Bessel functions J ν 4. The Bessel functions Y ν 5. Relations between J ν and
More informationAn inverse elliptic problem of medical optics with experimental data
An inverse elliptic problem of medical optics with experimental data Jianzhong Su, Michael V. Klibanov, Yueming Liu, Zhijin Lin Natee Pantong and Hanli Liu Abstract A numerical method for an inverse problem
More informationTomography and Reconstruction
Tomography and Reconstruction Lecture Overview Applications Background/history of tomography Radon Transform Fourier Slice Theorem Filtered Back Projection Algebraic techniques Measurement of Projection
More informationENGI 4430 PDEs - d Alembert Solutions Page 11.01
ENGI 4430 PDEs - d Alembert Solutions Page 11.01 11. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives
More informationInverse Problems: From Regularization to Bayesian Inference
1 / 42 Inverse Problems: From Regularization to Bayesian Inference An Overview on Prior Modeling and Bayesian Computation Application to Computed Tomography Ali Mohammad-Djafari Groupe Problèmes Inverses
More informationStatistical Inverse Problems and Instrumental Variables
Statistical Inverse Problems and Instrumental Variables Thorsten Hohage Institut für Numerische und Angewandte Mathematik University of Göttingen Workshop on Inverse and Partial Information Problems: Methodology
More informationM.Sc. in Meteorology. Numerical Weather Prediction
M.Sc. in Meteorology UCD Numerical Weather Prediction Prof Peter Lynch Meteorology & Climate Centre School of Mathematical Sciences University College Dublin Second Semester, 2005 2006. In this section
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 informationarxiv: v1 [math.na] 30 Jan 2018
Modern Regularization Methods for Inverse Problems Martin Benning and Martin Burger December 18, 2017 arxiv:1801.09922v1 [math.na] 30 Jan 2018 Abstract Regularization methods are a key tool in the solution
More informationAn Inversion Method for the Exponential Radon Transform Basedonthe Harmonic Analysis of the Euclidean Motion Group
An Inversion Method for the Exponential Radon Transform Basedonthe Harmonic Analysis of the Euclidean Motion Group Can Evren Yarman and Birsen Yazıcı Department of Electrical, Computer and System Engineering,
More informationLocal Tomography in 3-D SPECT
Pubblicazioni del Centro De Giorgi Proceedings of the Interdisciplinary Workshop on Mathematical Methods in Biomedical Imaging and Intensity-Modulated Radiation Therapy (IMRT), Pisa, Italy, October 2007.
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes Fachrichtung 6.1 Mathematik Preprint Nr. 362 An iterative method for EIT involving only solutions of Poisson equations. I: Mesh-free forward solver Thorsten Hohage and Sergej
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 informationDiscretization-invariant Bayesian inversion and Besov space priors. Samuli Siltanen RICAM Tampere University of Technology
Discretization-invariant Bayesian inversion and Besov space priors Samuli Siltanen RICAM 28.10.2008 Tampere University of Technology http://math.tkk.fi/inverse-coe/ This is a joint work with Matti Lassas
More informationResearch Article Exact Interior Reconstruction with Cone-Beam CT
Biomedical Imaging Volume 2007, Article ID 10693, 5 pages doi:10.1155/2007/10693 Research Article Exact Interior Reconstruction with Cone-Beam CT Yangbo Ye, 1 Hengyong Yu, 2 and Ge Wang 2 1 Department
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 informationStatistical regularization theory for Inverse Problems with Poisson data
Statistical regularization theory for Inverse Problems with Poisson data Frank Werner 1,2, joint with Thorsten Hohage 1 Statistical Inverse Problems in Biophysics Group Max Planck Institute for Biophysical
More informationP(E t, Ω)dt, (2) 4t has an advantage with respect. to the compactly supported mollifiers, i.e., the function W (t)f satisfies a semigroup law:
Introduction Functions of bounded variation, usually denoted by BV, have had and have an important role in several problems of calculus of variations. The main features that make BV functions suitable
More informationSeismic imaging and optimal transport
Seismic imaging and optimal transport Bjorn Engquist In collaboration with Brittany Froese, Sergey Fomel and Yunan Yang Brenier60, Calculus of Variations and Optimal Transportation, Paris, January 10-13,
More information1. Partial differential equations. Chapter 12: Partial Differential Equations. Examples. 2. The one-dimensional wave equation
1. Partial differential equations Definitions Examples A partial differential equation PDE is an equation giving a relation between a function of two or more variables u and its partial derivatives. The
More informationPreconditioning. Noisy, Ill-Conditioned Linear Systems
Preconditioning Noisy, Ill-Conditioned Linear Systems James G. Nagy Emory University Atlanta, GA Outline 1. The Basic Problem 2. Regularization / Iterative Methods 3. Preconditioning 4. Example: Image
More informationDiffusion for a Markov, Divergence-form Generator
Diffusion for a Markov, Divergence-form Generator Clark Department of Mathematics Michigan State University Arizona School of Analysis and Mathematical Physics March 16, 2012 Abstract We consider the long-time
More informationThe spherical mean value operator with centers on a sphere
The spherical mean value operator with centers on a sphere David Finch Department of Mathematics Oregon State University Corvallis, OR 9733-465 Email: finch@math.oregonstate.edu Rakesh Department of Mathematics
More informationSome Topics in Stochastic Partial Differential Equations
Some Topics in Stochastic Partial Differential Equations November 26, 2015 L Héritage de Kiyosi Itô en perspective Franco-Japonaise, Ambassade de France au Japon Plan of talk 1 Itô s SPDE 2 TDGL equation
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More information