first CT scan
ELECTRIC & MUSICAL INDUSTRIES
imaging
diagnostics
data-driven imagingautomatic analysis and rec.societal implicationsdata-driven imagingautomatic analysis and rec.societal implicationsdata-driven imagingautomatic analysis and rec.societal implications1920
1920
1958
1920
1958
2012
2012
2023
1920
1958
radiopedia.org
data-driven imagingautomatic analysis and rec.societal implicationsdata-driven imaging
automatic analysis and rec.
societal implications
\(x:\) phenomenon of interest
\(y = A(x) + v:\) measurements
How do we estimate \(x\) from \(y\)?
\(\mathbb{R}^n\)
\(A : \mathbb{R}^n \to \mathbb{R}^m\)
\(x\)
\(y\)
\(\mathbb{R}^m\)
inversion
Computed Tomography
inversion
\(A : \mathbb{R}^n \to \mathbb{R}^m\)
\(\mathbb{R}^n\)
\(x\)
\(y\)
\(\mathbb{R}^m\)
inversion
\(A : \mathbb{R}^n \to \mathbb{R}^m\)
\(\mathbb{R}^n\)
\(x\)
\(y\)
\(\mathbb{R}^m\)
Computed Tomography
inversion
Imaging
inversion
[Xiang et al, 2024]\(A : \mathbb{R}^n \to \mathbb{R}^m\)
\(\mathbb{R}^n\)
\(x\)
\(y\)
\(\mathbb{R}^m\)
inversion
\(A : \mathbb{R}^n \to \mathbb{R}^m\)
\(\mathbb{R}^n\)
\(x\)
\(y\)
\(\mathbb{R}^m\)
Computed Tomography
inversion
Imaging
Astronomy
inversion
[Event Horizon Telescope Collaboration, Astrophys. J. Lett. 875 (2019).]A problem is a well-posed problem if
1. (Existence) There exist at least one solution
2. (Uniqueness) There is at most one solution
3. (Stability) The solution depends continuously on the data
Inverse problem
\[\text{find }x \text{ so that }y = A x + v\]
Can we invert the operator \(A\)?
\[\text{find }x \text{ so that }y = A x + v\]
(unstable)
Error \(\|x - \hat x\|_2 \leq \|A^{-1}\|_2\|v\|_2 \propto \frac{1}{\lambda_n} \|v\|_2\)
\[\hat x = A^{-1}y = x + A^{-1}v\]
\(1~ \bullet~\) If \(n=m\) and \(A:\) non-singular:
\(=\)
\(+\)
\(2~ \bullet~\) If \(n<m\), then \(A\) cannot span \(\mathbb R^m\).
Thus, \(y \notin \text{span}(A)\), and thus \(\nexists ~ x : y = Ax\).
(non-existence)
\(=\)
\(+\)
\[\hat x = A^{-1}y = x + A^{-1}v\]
\(1~ \bullet~\) If \(n=m\) and \(A:\) non-singular:
\[\text{find }x \text{ so that }y = A x + v\]
(non-uniqueness)
\(3~ \bullet~\) If \(n>m\), and \(A:\) full-rank.
\(=\)
\(+\)
\(2~ \bullet~\) If \(n<m\), then \(A\) cannot span \(\mathbb R^m\).
Thus, \(y \notin \text{span}(A)\), and thus \(\nexists ~ x : y = Ax\).
\[\hat x = A^{-1}y = x + A^{-1}v\]
\(1~ \bullet~\) If \(n=m\) and \(A:\) non-singular:
\[\text{find }x \text{ so that }y = A x + v\]
(non-uniqueness)
\(3~ \bullet~\) If \(n>m\), and \(A:\) full-rank.
\(=\)
\(+\)
\[\min_x ~R(x) ~~ \text{s.t.}~~ y = Ax\]
Strategy: define a new criteria, \(R(x): \mathbb R^n \to \mathbb R\)
\[\text{find }x \text{ so that }y = A x + z\]
degree to which a material is magnetized when placed in a magnetic field
[from talk by S. Bollmann]
[Li et al, 2012]
[Lai, Aggarwal, van Zijl, Li, Sulam, Learned Proximal Networks for Quantitative Susceptibility Mapping, MICCAI 2020]
data-driven imaging
automatic analysis and rec.
societal implications
Multiaccuracy:
Example: Predicting abnormal findings in chest X-raysMultiaccuracy:
Example: Predicting abnormal findings in chest X-raysProblem: This is not always possible...
Problem: This is not always possible...
Problem: This is not always possible...
not computable
computable gaurantee
not computable
computable gaurantee
correcting multi-group accuracy
data-driven imaging
automatic analysis and rec.
societal implications
Formal frameworks for interpretability for decision making
Understanding social implications of algorithms in the wild
Efficient and robust generative AI
Biomarker discovery through AI
many open questions...