hr_fft.pro

来自「IDL语言编写的用于天文自适应光学仿真的软件CAOS V6.0的第一部分。」· PRO 代码 · 共 31 行

; $Id: hr_fft.pro,v 1.2 2002/03/14 11:49:12 riccardi Exp $

;+
; HR_FFT
;
; high resolution FFT
;
; Result = HR_FFT(Fun, Direction, Res)
;
; INPUT
; Fun:	matrice o vettore della funzione da trasformare
; Direction: +/-1 come per FFT
; Res:	intero > 1
;
; OUTPUT
; vettore o matrice complessa con dimensioni Res volte quelle di Fun.
; contiene la FFT di una estensione di Fun:
; Fun matrice NxM, allora Result matrice Res*NxRes*M che e` la FFT
; di Fun' con Fun'(0:N-1,0:M-1) = Fun e per gli altri elementi Fun' = 0.
;+ 
function hr_fft, array, direction, resolution
 
    nx = n_elements(array(*,0))
    ny = n_elements(array(0,*))
 
    hr_array = complexarr(nx*resolution,ny*resolution)
    hr_array(0:nx-1,0:ny-1) = array
    return, fft(hr_array, direction)
end