comp.lang.idl-pvwave archive: archive » Radon forward projection problem

Home » Public Forums » archive » Radon forward projection problem

Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend

Radon forward projection problem [message #59931]

Tue, 22 April 2008 08:07

Wox
Messages: 184
Registered: August 2006

Senior Member

Hi All,

I'm trying to implement an alternative for the FBP(filtered
backprojection) method for reconstructing objects measured in
tomography experiments. It's supposed to give less noisy tomograms.

Anyway, it's called OSEM and it's some iterative procedure using
forward and backward projection until the real sinogram and the
calculated sinogram are close.

For the projection, I use IDL's radon function. But I noticed
something strange with the forward projection. Try the code below. It
calculates the sinogram of a tomogram which is an image with all
pixels equal to 1. If you look at the result, something strange is
going on in the corners of the sinogram image. Does anyone know what
causes it and whether it is an intrinsic radon transform problem?

I would like to get rid of it, because this "estsinogram" is
calculated in each iteration of the OSEM (only in the first iteration
on an image with 1's) and used to normalize the measured sinogram
before adapting the tomogram. The resulting tomogram has some
artifacts because of it.

Thanks in advance,

Wout

pro test
; Detector
N=80
projcen=(N-1)/2.

; Angles
anglestart=0.
anglerange=180.
NAngles=anglerange/2.

angleinc=anglerange/(NAngles-1)
angles=anglestart+angleinc*findgen(NAngles)
angles*=!pi/180

; Reconstructing an object with 1's
tomogram=replicate(1.,N,N)
estsinogram=radon(tomogram,theta=angles,xmin=-projcen,RMIN=- projcen,drho=1,NRHO=N,/LINEAR)

loadct,0
window
tvscl,not bytscl(rebin(estsinogram,NAngles*3,N*3,/sample))
end

Report message to a moderator

[Message index]

		Radon forward projection problem By: Wox on Tue, 22 April 2008 08:07
		Re: Radon forward projection problem By: Wox on Thu, 24 April 2008 04:14
		Re: Radon forward projection problem By: mmiller3 on Thu, 24 April 2008 09:22
		Re: Radon forward projection problem By: mmiller3 on Wed, 23 April 2008 09:53
		Re: Radon forward projection problem By: Vontobel Peter on Wed, 23 April 2008 03:42
		Re: Radon forward projection problem By: Wox on Wed, 23 April 2008 03:02
		Re: Radon forward projection problem By: Vontobel Peter on Wed, 23 April 2008 00:42
		Re: Radon forward projection problem By: Wox on Sun, 27 April 2008 04:16
		Re: Radon forward projection problem By: mmiller3 on Mon, 28 April 2008 06:33

Previous Topic:	Re: Some AVI bother...
Next Topic:	idlanroigroup masking weirdness

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Sat Nov 29 04:27:04 PST 2025

Total time taken to generate the page: 2.24376 seconds