Graduate Performing MLEM with Response Matrix for 4pi Imager and 1x3 Vector Input

  • Thread starter Thread starter a.glick
  • Start date Start date
  • Tags Tags
    Image
Click For Summary
The discussion focuses on using the MLEM algorithm with a specific response matrix for a 4pi imager, which involves a complex dimensionality of observation and source bins. The user seeks clarification on how to perform MLEM iterations with a 1x3 vector input while needing a vectorized image in a 1x(72x36) format. There is a request for more information to facilitate better responses, and a suggestion to provide a link explaining MLEM for clarity. Additionally, there is a query about the concept of 4pi kinematic imaging. Understanding these elements is crucial for effectively implementing MLEM in this context.
a.glick
Messages
4
Reaction score
0
TL;DR
I am studying 4pi kinematic imaging with a neutron scatter camera and want to perform MLEM given det1, det2, and theta values for my cones.
I have a response matrix for my 4pi imager to perform MLEM that has dimensions:

(16 x 16 x 18) x (72 x 36) which is Observation Bins x Source Bins.

I'm told that I can use MLEM to create an image by feeding in a 1x3 vector of

det1, det2, theta

into the MLEM algorithm. Where det1, det2 are the 2 detectors a particle interacts with in the kinematic imaging and theta is the opening angle of a cone of probability for the particle's incident direction. However every paper I've read says I need a vectorized image (or simple backprojection) that has a 1 x (72 x 36) dimensional form. How can I do the MLEM iteration in my case?
 
Physics news on Phys.org
You'll have a much better chance of getting answers to your questions if you provide more information.
  • What is MLEM? We shouldn't need to do a search to find what this algorithm is. A link to a description of this algorithm would be helpful.
  • What is 4pi kinematic imaging?
 
  • Like
Likes WWGD
I am studying the mathematical formalism behind non-commutative geometry approach to quantum gravity. I was reading about Hopf algebras and their Drinfeld twist with a specific example of the Moyal-Weyl twist defined as F=exp(-iλ/2θ^(μν)∂_μ⊗∂_ν) where λ is a constant parametar and θ antisymmetric constant tensor. {∂_μ} is the basis of the tangent vector space over the underlying spacetime Now, from my understanding the enveloping algebra which appears in the definition of the Hopf algebra...