zap(苏黎世大气净化)是一种高精度的天空减法工具。
zap的Python项目详细描述
zap(苏黎世大气净化)
厌倦了天空减法残差?干掉他们
zap是一个高精度的天空减法工具,可以作为完整的天空。 减法解决方案,或作为先前天空减法缪斯的增强 数据。该方法利用主成分分析(pca)分离残差天空减法特征和 从观察到的数据多维数据集中删除它们。未来的发展将包括 用于各种仪器的改装。
zap的最后一个稳定版本可以简单地通过pip安装:
pip install zap
或进入用户路径:
pip install --user zap
引文
描述原始方法的论文可以在这里找到: http://adsabs.harvard.edu/abs/2016MNRAS.458.3210S
请引用zap作为:
\bibitem[Soto et al.(2016)]{2016MNRAS.458.3210S} Soto, K.~T., Lilly, S.~J., Bacon, R., Richard, J., \& Conseil, S.\ 2016, \mnras, 458, 3210
2.1(2019-07-03)
- Zap现在需要Python 3.5或更高版本,Astropy 2.0或更高版本
- 与NFM-AO模式兼容。
- 忽略“中值”和“拟合”方法中的陷波滤波器区域 连续过滤(cftype)
- 将默认的连续过滤方法(cftype)更改为“中值”,然后 删除“权重”方法。以前的默认方法('weight')有一些 谱边处理和陷波滤波器边问题 对于ao模式,也引入了红光波长的起伏。
- 为nevals添加命令行参数
- 需要SCIPY 0.18.1或更高版本,因为 以前的版本
- 添加一个函数来屏蔽立方体边(zap.mask_nan_edges):spaxel 立方体边缘的光谱中可能有很多NaN,所以ZAP 不减去这些留下高残差的天空。
2.0(2017-09-08)
与WFM-AO模式兼容
对PCA使用Scikit-learn的实现,而不是自定义的实现。 这解决了zap引入的空间变化问题。它也是 比以前的实现快得多。但缺点是 无法再保存svd文件。
默认情况下,只使用一个sky segment,这意味着立方体是 在波长轴上不再分裂。最初Zap使用了11个片段, 我们的目标是要有一组连贯的天空发射线 每段特征值的数目。它还允许并行化 计算。但是,这些节段也负责连续体 振荡,并选择每段的特征值数目 非常困难和敏感只有一个环节 天空减法更好,因为 整个波长范围内的天空线。
因此,只使用一个片段可以大大降低 新PCA的发射线。如果 需要:
from zap.zap import SKYSEG SKYSEG[:] = [0, 5400, ..., 10000]
具有多项式拟合(cftype='fit')的新连续体筛选器类型必须是 但要小心使用,因为红色部分的配合很容易失控 在光谱中。
将median和 weight连续过滤器。以前使用的值,100和50,也是 很小,这解释了 光谱。
新参数(ncpu)用于设置使用的CPU数量
删除将SVD结果保存到FITS文件中的可能性。这是因为 PCA实现的变化仍然有可能通过SVD 由SVDoutput到process计算,传递内存中对象 直接保存到磁盘
速度的提高,主要归功于新的PCA实现
新参数保存解释的方差曲线。
1.0(2016-04-02)
首次公开发行。这是2016MNRAS.458.3210S中描述的版本 纸张。