原来发过一个Python下的换桌面脚本,由于相关网站改版,脚本失效.。我也相应更新一下。移除了图片质量一般的EPOD。充分测试了windows(XP) 和 linux(Ubuntu gnome)双平台. 同时将原有的拉伸全屏模式,更改为保持长宽比的居中模式(其他模式的使用参数也已写在注释中)。 对于APOD,可以手动指定更新任意一天的图片(只支持JPG), 而NGPOD的网站地址与时间无关,只能获取最新图片。下载前会检查文件名,如果图片已经存在,则自动跳过。默认下载大图,然后根据设定分辨率生成桌面图片。各项参数要在源码中自行调整。仍需要PIL模块支持,这个模块的详细介绍可以看这里。新脚本点此下载。
继续阅读
Sherpa表面亮度拟合
与致力成为X射线数据分析通用工具的HEASOFT不同,Ciao (Chandra Interactive Analysis of Observations) 是专为Chandra卫星数据开发的。它也包含了一个数据分析工具Sherpa。基于现有的脚本语言,从而可以方便地移植扩充。开始有Python 和S-Lang 两个版本,不过现在已经停止了对S-Lang的支持。下面我们介绍如何用它来拟合星系团的面亮度分布。
为了减少计算量,我们先从去除点源的图像中裁出星系团所在区域.使用 dmstat处理前面生成的图像文件.可以得到图像中的极大极小值坐标,如果目标形状不规则,还可以设置 centroid =yes 计算图像重心.
继续阅读
XSPEC光谱分析
Chandra图像定标
这是2008年Paolo Tozzi 来北京做的X射线数据处理讲座的笔记,一直拖到现在才开始整理。这里是第一部分,主要是些准备工作。高能数据处理要用HEASOFT,Chandra数据还要用到Ciao和CALDB,按照说明文件一步步安装并不困难,只是2G多文件下载起来要费些时间。附加脚本也经常用到,解压到安装路径下就可以了。安装好之后,先用 heainit 和ciao初始化运行环境,若没有错误提示就可以运行下面的命令了。
如果手头上没有数据。可以先到Chandra数据中心找个自己喜欢的目标,数据通常会在保密期过后自动公开,既避免了对同一目标的反复观测,也能最大限度地发挥已有数据的价值。 我这里选择了红移0.3的 MACSJ0257.6-2209 (即Abell 402, Obs ID 3267 )。输入到左上角的Target Name中,点Search,会得Chandra已有的观测数据列表。我们选个曝光时间短点的,数据文件小些,加到获取列表中(Add to Retrieval List)。点击retrieval之后,服务器会自动将文件打包放到匿名ftp的临时目录下。也有一个基于python的自动下载脚本。 官方详尽的下载介绍由此跳转。
继续阅读