comp.lang.idl-pvwave archive: archive » Hankel (Fourier-Bessel) Transform

Home » Public Forums » archive » Hankel (Fourier-Bessel) Transform

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

Re: Hankel (Fourier-Bessel) Transform [message #25914 is a reply to message #25846]

Tue, 24 July 2001 12:15

thompson
Messages: 584
Registered: August 1991

Senior Member

Georg.Pabst@nrc.ca (Georg Pabst) writes:

> Hi,

> I'm looking for the Hankel (Fourier-Bessel) Transform, i.e.,
> int(0,infinity) f(t)*BesselJ(t*r)t dt being implemented in IDL.

> There is a paper "Siegman A. 1980. Quasi fast Hankel transform. Opt.
> Lett. 1, 13-15" and one can also find the code in Fortran or C...

> Thanks,
> Georg

Here's an old program that I think might be what you need.

Bill Thompson

FUNCTION HANKEL,F
;
; This function returns the Hankel transform of the argument.
;
S = SIZE(F)
IF S(0) NE 1 THEN BEGIN
PRINT,'*** Variable must be a one-dimensional array, name= F, routine HANKEL.'
RETURN,F
ENDIF
;
X = INDGEN(F)
K = ( 2. * !PI / FLOAT(N_ELEMENTS(X)) ) * X
SC = 0.*X + 1.
IF N_ELEMENTS(SC) GT 3 THEN BEGIN
SC(0) = 3.D0 / 8.D0
SC(1) = 7.D0 / 6.D0
SC(2) = 23.D0 / 24.D0
ENDIF
;
H = BES0( K # X ) # ( K * F * SC )
;
RETURN,H
END

Report message to a moderator

[Message index]

		Hankel (Fourier-Bessel) Transform By: Georg.Pabst on Thu, 19 July 2001 10:53
		Re: Hankel (Fourier-Bessel) Transform By: Georg.Pabst on Wed, 25 July 2001 10:23
		Re: Hankel (Fourier-Bessel) Transform By: Craig Markwardt on Tue, 24 July 2001 12:46
		Re: Hankel (Fourier-Bessel) Transform By: thompson on Tue, 24 July 2001 12:15
		Re: Hankel (Fourier-Bessel) Transform By: Mats Löfdahl on Wed, 26 March 2014 09:25
		Re: Hankel (Fourier-Bessel) Transform By: Mats Löfdahl on Thu, 27 March 2014 08:12
		Re: Hankel (Fourier-Bessel) Transform By: Helder Marchetto on Thu, 27 March 2014 08:41
		Re: Hankel (Fourier-Bessel) Transform By: Mats Löfdahl on Thu, 27 March 2014 10:24

Previous Topic:	CGMAP_GRID: Lattitude lines
Next Topic:	IDL Indexing 2D->3D

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Wed Oct 08 16:15:00 PDT 2025

Total time taken to generate the page: 0.00185 seconds