随着工作学习的深入，和原始天文照片的接触越来越多，发现对天文专业来说数字图像处理还是一门很有用的基础课程，连CMB功率谱之类看上去高深莫测的技术，其实都写在信号处理专业的本科教材里，把这段流程完全交给现成软件或者编程人员是没办法真正理解观测数据的。开始补课了，先从找星星开始～
比如我们有上面这张照片, 图中的亮点就是实际拍摄的星空，但左中和右下两个最亮的白点都是打在CCD上的宇宙线，它们能量很集中，没有扩展的形状。那要怎么让程序找到这些天体的位置，并识别出正确的星点呢？ 继续阅读
标签存档： code
用Matlab做星点识别
比如我们有上面这张照片, 图中的亮点就是实际拍摄的星空，但左中和右下两个最亮的白点都是打在CCD上的宇宙线，它们能量很集中，没有扩展的形状。那要怎么让程序找到这些天体的位置，并识别出正确的星点呢？ 继续阅读
天文Latex模版
Latex可以精准的控制输出效果，早已成为科技界印刷排版的通用语言，但它并不像Word那样直观，初学者要熟悉不少的命令才能获得理想的效果。而我们要熟悉的软件和语言已经太多了，没有多少精力再来研究格式排版，但是如果不知道如何使用式样模版，总是自己调整字体、分栏、页边距……，或者不了解BibTex，总要花半天时间逐个调整参考文献格式，这样无疑失去了用Latex的意义，因此有人宁愿先用Word编辑输入，再通过Scientific Workplace转成Latex文件。其实Latex还是非常友好的，虽然参数选项众多，但都有详尽的帮助说明；虽然采用命令行编译，但有完整的日志记录；就是安装配置稍微麻烦一些，人家也做好了CTex安装包……还有什么可抱怨呢？如果这些理由仍不能说服你放弃doc，请看ArXiv FAQ： Why Tex ？
当前的主流天文学术期刊有：美国天文学会主办的AJ(The Astronomical Journal)、ApJ(The Astrophysical Journal)、欧南台主导的A&A (Astronomy & Astrophysics)，英国皇家天文学会主办的MNRAS (Monthly Notices of the Royal Astronomical Society)，以及美国物理学会出版的PRD(Physics Review D)，PRL(Physics Review Letters)等；至于受众较广的Nature、Science和Scientific American 之类则有专门的美编统一版面，不需作者费心。因此Nature仍采用Word编排（由于兼容性问题不接受Office2007格式），其它则都有公布官方的Latex模版文件。AJ、APJ及同系列的 APJL、ApJS都采用AASTex模版；PRD、PRL则采用美国物理学会统一的RevTex模版；至于Science、MNRAS和AA则也都有自己的格式。除了官方指定的模版，还有相当多用户自行制作的式样，比如针对ApJ的emulateapj，CTAN模版库里还有很多，有需要就自己去看吧。
由度规算曲率标量
上周被某人说成是“不务正业”，这次就讲讲本行吧(忽悠新手，外行可以无视，同行尽管鄙视)。
当代宇宙学的理论基石是在Robertson-Walker度规下由Einstein场方程推导出的Fridemman方程，现在的绝大部分工作都是以此为起点，而且也得到了可靠的观测支持。从给定度规得出场方程的具体形式是广义相对论的基础内容，但是由于方法繁复，教科书中都不会给出具体的计算过程，而结果又不是一望便知的，学到这里谁都少不了课下的一番推导验证。我当年偷懒跳过，现在却发现自己怎么都算不对了……
怅惘之际在论文库里发现2000年《上海天文台年刊》第21期中有一篇《利用Mathematica软件表示真空Einstein场方程》的文章，但又不想为此学门新语言，便用Matlab仿写了一个，代码如下：
