Radon (slant stack or tau-p) transform

[Zack]

Hello all...

Just thought I would see if anybody has any code that performs a
radon (also called a slant stack or tau-p) transform on a matrix. Or
some suggestions on how to code it up better.

I'm applying the transform to seismic data.

Like the name suggest, I want to sum all values along a linear line.
I want to do this for a range of slopes (equal to the number of
columns) and for each intercept time (number of rows). So a straight
line in the original matrix (time-distance domain) will become a
point in the transformed matrix (tau-slowness domain).

I have not yet found an effective way to account for the times when
the chosen linear line to stack values along does not fall exactly on
a data point. Do I interpolate between the values above and below?
Or do I just use the nearest value?

Thanks in advance...

[Zack]

*bump*

Just thought a little bump may help...thanks!!!!

[Steve Eddins]

See the algorithm description for the Image Processing Toolbox radon
function:

http://www.mathworks.com/access/helpdesk/help/toolbox/images/radon.html

Steve Eddins
image processing blog - http://blogs.mathworks.com/steve

[Lara Bertello]

Zack <z...@rock.com> wrote in message <ef35...@webx.raydaftYaTP>...

> *bump*
>
> Just thought a little bump may help...thanks!!!!

Hello Zack,

Did you find a simple code for the slant stack trasform? I am working on seismic data too. I would like to know all the steps for anaylise the field data and to obtain the dispersion curve.
Do you know something about that?

Thank you.
Lara