最近看到一篇文章

“Common-Resolution Convolution Kernels for Space- and Ground-Based Telescopes”

讲的是不同的望远镜有不同的点扩散函数（PSF），需要通过卷积核将不同望远镜观测的图像，转换相同点扩散函数的图像。这样可以实现多波段的协作，用来计算光谱能量的分布（SED）。

他们的这项工作还是非常漂亮的，给出IDL程序和一些卷积核。

使用程序convolve_image.pro，需要把

convol_fft.pro

放到目录/usr/local/itt/idl71/lib/，因为convolve_image.pro会调用convol_fft.pro

下面是convolve_image.pro的使用方法

首先我在http://hedam.lam.fr/VNGS/data.php 下载了M83的herschel光谱和光度影像接收机SPIRE的500微米波段的数据

将SPIRE 500 微米波段的数据转换到MIPS 160微米波段

我们需要在http://www.astro.princeton.edu/~ganiano/Kernels/Ker_2012_April_Spec/Kernels_fits_Files/Hi_Resolution/ 下载卷积核：Kernel_HiRes_SPIRE_500_to_MIPS_160.fits.gz

在命令行中依次执行下面的命令

.compile convolve_image

image=readfits('/root/kernel/M83-2C78-PLWmap_MS-v51-Brigade.fits',header,EXTEN_NO=1)

kernel_image=readfits('/root/kernel/Kernel_HiRes_SPIRE_500_to_MIPS_160.fits',kernel_header)

do_we_write=1

do_the_convolution,image,header,kernel_image,kernel_header,result_image,result_header,result_kernel_image,result_kernel_header,do_we_write

WRITEFITS,"/root/kernel/data/M83-SPIRE_500_to_MIPS_160.FITS", result_image, result_header

WRITEFITS,"/root/kernel/data/kernelSPIRE_500_to_MIPS_160.FITS", result_kernel_image, result_kernel_header

左边的图是原始图像，右边的图是转换后的图像。