← All work

Research

CTsynth

What can spectral CT measurements tell us—and what remains unidentifiable?

  • Inverse problems
  • CT physics
  • Numerical methods
On this page
  1. The problem
  2. The mathematical idea
  3. Implementation and validation
  4. Identifiability and limitations
  5. Next questions

The problem

A forward model predicts measurements from a material distribution. The inverse problem asks which material distributions are consistent with the measurements. These questions are related, but a differentiable simulator does not automatically make the inverse problem well posed.

The mathematical idea

A schematic spectral transmission model integrates energy-dependent attenuation along a ray:

Ir=∫S(E)exp⁡[−∫rμ(x,E) dx]dE.I_r=\int S(E)\exp\left[-\int_r\mu(x,E)\,dx\right]dE.

Here S(E)S(E) represents the incident spectrum and detector weighting, and μ(x,E)\mu(x,E) is the attenuation field. This expression illustrates the modeling question; it is not a specification of every implementation detail.

Implementation and validation

CTsynth explores differentiable spectral CT physics, forward projection, numerical adjoints, and material decomposition. An important numerical check compares the forward operator and its adjoint through the inner-product identity ⟨Ax,y⟩=⟨x,A∗y⟩\langle Ax,y\rangle=\langle x,A^*y\rangle. Gradient checks and decomposition benchmarks address different parts of the problem.

Identifiability and limitations

A small measurement residual need not imply correct material recovery. Spectral information, measurement geometry, noise, and the chosen material basis all affect what can be recovered. Benchmark values and experiment figures are not included here because validated artifacts have not yet been supplied.

Next questions

Which ambiguities arise from the physics, and which come from discretization or optimization? Which additional measurements would distinguish otherwise similar material configurations?

← Back to WorkDiscuss this project ↗