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

Home » Public Forums » archive » RADON/RIENMANN TRANSFORMATION
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/RIENMANN TRANSFORMATION [message #39835 is a reply to message #39764] Mon, 21 June 2004 10:05 Go to previous messageGo to previous message
Timm Weitkamp is currently offline  Timm Weitkamp
Messages: 66
Registered: August 2002
Member
Max:

> -> is it correct to say :
> -> the sinogram will be 128 images of 128*360 matrix

Yes.

> - > From the discussion the Theta should be an array of 360 elements
,(one
> element per angle )

Yes. If the angles are equidistant with 1 degree spacing, then you might
set

IDL> theta = FINDGEN(360) * !DTOR

> -> For every Theta there should be an assosiated RHO , which means it
should
> also be an array of 360 elements.

No! RHO is the spatial coordinate and is therefore, in your case, a vector
of 128 elements. For example:

IDL> rho = FINDGEN(128) - (127 / 2.0)

For the backprojection, you can then call RADON with the previously
defined variables "theta" and "rho" and the sinogram data. For example:

IDL> rec = RADON(/BACKPROJECT, sinogram $
, THETA=theta, RHO=rho, NX=128, NY=128).

(where "sinogram" should be the array containing your sinogram [which
you will probably have to filter before]).

Good luck
Timm

On 19.06.04 at 16:47, 555777555 wrote:

> Thanks for everyones answers ans suggestions ,i appeciate it very much
> 1. Assuming the data is aquired in 128*128 matrix , over 360 degrees at 1
> image per angle
> (for simplicity)
> -> is it correct to say :
> -> the sinogram will be 128 images of 128*360 matrix
> - > From the discussion the Theta should be an array of 360 elements ,(one
> element per angle )
> -> For every Theta there should be an assosiated RHO , which means it should
> also be an array of 360 elements.
>
> I have tried using the above but IDL respond that the dimentions do not
> agree
> ? What is the correct way / Where am i wrong
> Max

--
Timm Weitkamp <http://people.web.psi.ch/weitkamp>
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: serial port handling
Next Topic: Re: serial port handling

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

Current Time: Thu Oct 09 23:03:27 PDT 2025

Total time taken to generate the page: 1.43951 seconds