Recent content by MaximeMusterFrau

The MLEM is to be seen in the picture uploaded, where aij are the elements of the matrix A and y is to be considered b in our case,.

okey, a more accurate situation: we have A = [[ 5. 2. 3. ] [ 0.4 0.1 5. ] [ 5. 0.4 0.2] [ 0.3 0.2 10. ] [ 1. 0.3 2. ] [ 0.3 0.7 0.4]] b = [[ 2. ] [ 0.1] [ 0.4] [ 0.2] [ 0.3] [ 0.7]] and we are looking for x that satisfies Ax =b . In this situation the MLEM works very good and the x is...

Yes it is..

humm.. but I am a bit confused now.. because in the abstract it is cited that the elements of the system matrix also have to be positive. In the paper Pij>=0, which is equivalent to Aij >= 0 in the system : Ax = b . In my case Xj >=0 but Aij can be positive as well as negative .. Does this...

the MLEM is supposed to work for Ax = b, where : A: the probability matrix (it is given) x: the image (which we are trying to figure out, real values) b : given data (mesaered data) In my case I have a system where A is not really a probability matrix, because it has negative values (they are...

I am dealing with a different image processing procedure. The matrix is measured and obtained from a FFT Transformation. And I am trying to find a way to transform the system (A x = b ) so that the the MLEM works. In this case it doesn t work because the second derivative of the likelihood...

This equation is a modiefied version of the one I am working on. the C matrix is the one you find here : https://arxiv.org/ftp/arxiv/papers/1504/1504.06889.pdf it is the :probability of detecting an emission from the pixel j in projection’s bin i. So my question is, what to do if some...

@ohwilleke @Stephen Tashi : Thank you very much for your answer. with bild I meant image, I am sorry. It is about an image-reconstruction algorithm. http://ieeexplore.ieee.org/document/925290/ here is a link to a paper explaining what I am trying to do. In my case however, the c matrix is...

Hello everyone, I'have implemented a Maximum-Likelihood-Expectation-Maximization Algorithm in order to reconstruct a bild. let say, we have such a system Ax=b, where A is a complex matrix, b is a complex vector. A and b are known and we will iterately try to find the best x (which should be...