comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » Radon forward projection problem
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
Re: Radon forward projection problem [message #59981 is a reply to message #59931] Thu, 24 April 2008 04:14 Go to previous messageGo to previous message
Wox is currently offline  Wox
Messages: 184
Registered: August 2006
Senior Member
On Wed, 23 Apr 2008 12:53:05 -0400, mmiller3@iupui.edu (Michael A.
Miller) wrote:

>>>> >> "Wox" == Wox <nomail@hotmail.com> writes:
>
>> The algorithm I was talking about (OSEM, although I think
>> it's really called MLEM, I'm not sure) goes like this
>> (BP=backprojection, FP=forward projection):
>
> Be careful that you are sure of what you are implementing. OSEM
> and MLEM differ - especially in computational speed. The
> citations that I've got for the algorithms are:
>
> "Accelerated Image Reconstruction Using Ordered Subsets of
> Projection Data," IEEE Trans Med Img, 13, 601-609, 1994.
>
> "Maximum likelihood reconstruction for emission tomography,"
> IEEE Trans Med Img, MI-2, 113-122, 1982
>
>
> Mike

I have been reading those articles, but I couldn't understand the
difference. I would really appreciate a professional opinion on this.
As far as I understand, for each iteration in MLEM:
v'=v * BP(s0/FP(v))/BP(s1)

s0: experimental sinogram
s1: sinogram with all 1's
v: tomogram of previous iteration (this is a uniform image with 1's
for the first iteration)
v': the tomogram calculated for this iteration
BP: backprojection
FP: forward projection

Now what is OSEM doing?

Btw, is it correct that SIRT is doing this:
v'=v - b*BP(s0-FP(v))
where b a relaxation factor.
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: Re: Some AVI bother...
Next Topic: idlanroigroup masking weirdness

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Fri Oct 10 07:12:07 PDT 2025

Total time taken to generate the page: 0.96293 seconds