FUNCTION	PFFT,h,m

;Return FFT, padded with zeros to magnification factor m
;Useful for image interpolation, precise cross-correlation, etc.
;Create AHB 03/12

  if not KEYWORD_SET(m) then m=1
  sz = SIZE(h)
  fh = COMPLEXARR(sz[1]*m,sz[2]*m)
  fh(sz[1]*(m-1)/2 : sz[1]*(m+1)/2 - 1, $
     sz[2]*(m-1)/2 : sz[2]*(m+1)/2 - 1) = SHIFT(FFT(h),sz[1]/2,sz[2]/2)

RETURN,SHIFT(fh,sz[1]*m/2,sz[2]*m/2)
end