Information Geometry with Information Geometry with Applications in Components Applications in Components Analysis and X Analysis and X - - Ray CT Imaging Ray CT Imaging Joseph A. O Joseph A. O ’ ’ Sullivan Sullivan Electronic Systems and Signals Research Laboratory Department of Electrical and Systems Engineering Washington University [email protected]http://essrl.wustl.edu/~jao Supported by: ONR, ARO, NIH ONR, ARO, NIH
64
Embed
Information Geometry with Applications in Components ...jao/Talks/InvitedTalks/Duketalk033004.pdfInformation Geometry J. A. O’Sullivan. Duke Seminar, Mar. 30, 2004 21 More Information
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Information Geometry with Information Geometry with Applications in Components Applications in Components
Analysis and XAnalysis and X--Ray CT ImagingRay CT ImagingJoseph A. OJoseph A. O’’SullivanSullivan
Electronic Systems and Signals Research LaboratoryDepartment of Electrical and Systems Engineering
•• Lower rank Markov approximation Lower rank Markov approximation X X M M YYMM in a set of cardinality in a set of cardinality KK
•• Factor analysis, contingency tables, economicsFactor analysis, contingency tables, economics•• Problem: Approximation of one matrix by another of Problem: Approximation of one matrix by another of
lower ranklower rank•• C. C. EckartEckart and G. Young, and G. Young, PsychometrikaPsychometrika, vol. 1, pp. , vol. 1, pp.
211211--218 1936.218 1936.•• SVD SVD IVDIVD
AP Φ=
)||(min ASIA ΦΦ
∑=
=K
kkk yaxyxp
1)()(),( φ
J. A. O’Sullivan. Duke Seminar, Mar. 30, 2004Information Geometry
8
Spectral Components AnalysisSpectral Components Analysis• Data spectrum for each element (pixel)
• Model: linear combination of constituent spectra,
• Problem: Estimate constituents and proportions subject to nonnegativity; positivity of S assumed
• Comments: Radiometric Calibration; Constraints Fundamental
∑=
=K
kkjkj as
1φ AS Φ=
Ijs +ℜ∈
0≥kja
J. A. O’Sullivan. Duke Seminar, Mar. 30, 2004Information Geometry
9
Idealized Problem Statement:Idealized Problem Statement:MaximumMaximum--Likelihood Likelihood Minimum IMinimum I--divergencedivergence
•• Poisson distributed data Poisson distributed data loglikelihoodloglikelihood functionfunction
•• Maximization over Maximization over ΦΦ and A equivalent to and A equivalent to minimization of Iminimization of I--divergencedivergence
AS Φ=ˆ
∑∑ ∑∑= = == ⎪⎭
⎪⎬⎫
⎪⎩
⎪⎨⎧
−⎥⎥⎦
⎤
⎢⎢⎣
⎡=Φ
I
i
J
j
K
kkjik
K
kkjikij aasASl
1 1 11ln)|( φφ
∑∑ ∑∑= =
== ⎪⎭
⎪⎬⎫
⎪⎩
⎪⎨⎧
+−⎥⎥
⎦
⎤
⎢⎢
⎣
⎡=Φ
I
i
J
j
Kk kjikijK
k kjik
ijij as
a
ssASI
1 11
1
ln)||( φφ
•• Information Value Decomposition ProblemInformation Value Decomposition Problem
J. A. O’Sullivan. Duke Seminar, Mar. 30, 2004Information Geometry
10
HyperspectralHyperspectral ImagingImagingat Washington Universityat Washington University
Donald L. SnyderDonald L. SnyderWilliam H. SmithWilliam H. SmithDaniel R. Daniel R. FuhrmannFuhrmannJoseph A. OJoseph A. O’’SullivanSullivanChrysantheChrysanthe PrezaPreza
WU Team
Digital Array Scanning Interferometer (DASI)
J. A. O’Sullivan. Duke Seminar, Mar. 30, 2004Information Geometry
11
HyperspectralHyperspectral ImagingImaging• Scene Cube Data Cube• “Drink from a fire hose”• Filter wheel, interferometer,
tunable FPAs• Modeling and processing:
- data models- optimal algorithms- efficient algorithms
J. A. O’Sullivan. Duke Seminar, Mar. 30, 2004Information Geometry
12
CT Imaging in Presence of High CT Imaging in Presence of High Density AttenuatorsDensity Attenuators
J. A. O’Sullivan. Duke Seminar, Mar. 30, 2004Information Geometry
13
Filtered Back ProjectionFiltered Back Projection
Truth FBP
FBP: inverse Radon transform
J. A. O’Sullivan. Duke Seminar, Mar. 30, 2004Information Geometry
14
Transmission TomographyTransmission Tomography• Source-detector pairs indexed by y; pixels indexed by x• Data d(y) Poisson, means g(y:µ), log likelihood function
• Mean unattenuated counts I0, mean background β• Attenuation function µ(x,E), E energies
• Maximize over µ or ci; equivalently minimize I-divergence
)(),(),(exp),():(
):():(ln)()):(|(
0 yExxyhEyIyg
ygygydgdl
E x
y
βµµ
µµµ
+⎟⎟⎠
⎞⎜⎜⎝
⎛−=
−=⋅
∑ ∑
∑
∈
∈
X
Y
∑=
=I
iii ExcEx
1)()(),( µµ
J. A. O’Sullivan. Duke Seminar, Mar. 30, 2004Information Geometry
15
Maximum Likelihood Maximum Likelihood Minimum IMinimum I--DivergenceDivergence
Difficulties: log of sum, sums in exponent
)()()(),(exp),():(
):()():(
)(ln)()):(||(
):():(ln)()):(|(
10 yExcxyhEyIyg
ygydyg
ydydgdI
ygygydgdl
E x
I
iii
y
y
βµµ
µµ
µ
µµµ
+⎟⎟⎠
⎞⎜⎜⎝
⎛−=
+−=⋅
−=⋅
∑ ∑ ∑
∑
∑
∈ =
∈
∈
X
Y
Y
J. A. O’Sullivan. Duke Seminar, Mar. 30, 2004Information Geometry
16
OutlineOutline• Applications
- Components Analysis Information Value DecompositionHyperspectral Imaging
- Transmission Tomography- Maximum Likelihood Mixture
• Define problem as minq φ(q)• Derive variational representation: φ(q) = minp J(p,q)• J is an auxiliary function p is in auxiliary set P• Result: double minimization minq minp J(p,q)• Alternating minimization algorithm
J. A. O’Sullivan. Duke Seminar, Mar. 30, 2004Information Geometry
27
Alternating Minimization AlgorithmsAlternating Minimization Algorithms::II--Divergence, Linear, Exponential FamiliesDivergence, Linear, Exponential Families• Special Case of Interest: J is I-divergence• Families of Interest:
Linear Family L(A,b) = {p: Ap = b}Exponential Family E(π,B) = {q: qi = πi exp[Σj bij νj]}
)||(minarg
)||(minarg
)1(
)
)1(
)(
),(
)1(
qpIq
qpIp
l
B(q
l
l
bAp
l
+
∈
+
∈
+
=
=
,E
L
π
Csiszár and Tusnády; Dempster, Laird, Rubin; Blahut; Richardson; Lucy; Vardi, Shepp, and Kaufman; Cover;Miller and Snyder; O'Sullivan
J. A. O’Sullivan. Duke Seminar, Mar. 30, 2004Information Geometry