zgyin的个人博客分享 http://blog.sciencenet.cn/u/zgyin

博文

卷积核(多波段协作)

已有 3873 次阅读 2014-4-9 02:01 |个人分类:天文知识|系统分类:科研笔记

最近看到一篇文章

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


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



https://wap.sciencenet.cn/blog-783377-783215.html

上一篇:用epstopdf命令将eps转换为pdf格式
下一篇:idl 将两个直方图叠加到一块
收藏 IP: 159.226.171.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2023-2-5 08:29

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部