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 #60013 is a reply to message #59931] Wed, 23 April 2008 03:02 Go to previous messageGo to previous message
Wox is currently offline  Wox
Messages: 184
Registered: August 2006
Senior Member
On Wed, 23 Apr 2008 00:42:07 -0700 (PDT), VP <peter.vontobel@psi.ch>
wrote:

> Hi
>
> please compare your sinogram with the following:
>
> estsinogram=radon(tomogram,rho=rho,theta=theta,ntheta=nangle s)
>
> compare the rho and theta values and read the IDL radon help pages.
>
> HTH
>
> Peter

So basically what you're saying is: don't undersample.

However, now there are more zero's in the sinogram. Let me explain
what I want to do. 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):

========Pseudo-code========
sino1=sinogram with all 1's
tomo1=BP(sino1)

tomo=tomogram with all 1's
for i=0,niter-1 do begin
estsino=FP(tomo)
tomo = tomo * BP(sino_measured/estsino)/tomo1
endfor
========================

You see that "estsino" and "tomo1" can't have zeroed pixels. Off
course, I tried replacing the zeroed pixels by 1, max(rest),
min(rest), etc... But this gives some artifacts in the resulting
tomogram. Any ideas?
[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 10:49:37 PDT 2025

Total time taken to generate the page: 0.72110 seconds