用三次 GPI 曝光提取对比度光谱

从三个文件和原始观测夜开始

观测发生在二〇一三年十二月十日，项目 GS-ENG-GPI-COM。固定作者仓库提交为 72a7517b1a34bdc29c33211f2ef591ddff64dcf9，提交日期二〇一六年八月九日。三个文件编号 S0092、S0099、S0111，后接 _spdc_distorcorr.fits。这些文件在二〇一六年约简，观测、后来处理和论文发表属于不同日期。

每个原文件长 14,673,600 字节。下载器匿名请求原始头及十组 SCI、DQ 平面的精确字节范围，逐段检查 206 状态、Content-Range、长度和 SHA256。完整文件的 Git blob 标识也保存了，但部分下载没有被说成完整文件 SHA256。

保留的原始零基通道编号是 0、4、8、12、16、20、24、28、32、36。三份输入都是实际观测背景；之后的 PSF 注入才是单独标记的合成标定实验。

对应阅读：Bonnefoy et al. · Jason Wang and contributors · GPI Data Pipeline 1.4.0 data formats

先读单位和质量位，再看亮点

SCI 有三十七个原生波长平面，每面为二百八十一乘二百八十一像素。BUNIT 写的是每次 coadd 的 ADU，ITIME 为每次 88.74219 秒，COADDS 为一。起止窗口含有开销，直接取中点不保证等于电荷加权曝光中点；投影采用 CD 矩阵，头中 AVPARANG 用来检查视场转角跨度。

波长起点是 1.11407300105 微米，原生间距为 0.00654594640474 微米，FITS 参考像素从一开始。选定通道延伸到 1.34972707162 微米。原生采样比论文约三十五至三十九的分辨本领更密，十个保留通道不能当作十个独立光谱分辨单元。

DQ 扩展保存质量位。默认屏蔽坏像素和两种饱和位；读坡离群值剔除位、已修补或插值像素，只要没有同时标成基本坏像素，就暂保留。验证脚本还会尝试屏蔽所有非零质量位，若有效像素不足就报告测量失效。

文件没有给出校准方差的 ERR 扩展，经过约简的 ADU 不能直接当独立泊松电子。教程从处理后的比较区估条件空间光谱协方差，并保留模型欠拟合。

对应阅读：GPI Data Pipeline 1.4.0 data formats · Bonnefoy et al.

检查方向，再重测卫星斑

小天空偏移采用二维 CD 矩阵，列向量长度对应每原生像素十四点一四毫角秒。旧 CDELT 与它不符；PAR_ANG 注释还明说应使用 AVPARANG。FITS 的 CRPIX 从一计数，源读取程序的 PSFCENT 和 SATS 则按原生零基位置解释。

S0099、S0111 的 SATS 坐标字符串与 PSFCENT 完全重复。这些卡片原样保留。脚本对实际斑点像素重新拟合局部高斯和倾斜背景；自由迭代重心容易被邻近散斑吸走，因此搜索限制在头位置周围一个像素，边界解拒绝，并用独立经验 PSF 注入检查定位偏差。

独立定位检验把经验斑形加到真实光晕的其他已知位置。一百四十四次注入中，一百一十四次局部拟合通过，三十次失败；通过者的百分之九十五位置偏差分位约为零点八二五像素。失败和模板差异限制了真实卫星斑定位的精度解释。

每个通道以重测卫星斑的平均位置配准，再经原始 CD 投到北上东左网格。稀疏约简中，目标环带最强点的位置角约二百一十二点六二度，历史表一为二百一十一点六加减一点三度。这是像素尺度的独立方向对照，后来约简数据的精密测位还需自己的板尺与真北标定。

对应阅读：GPI Data Pipeline 1.4.0 data formats · Jason Wang and contributors · Bonnefoy et al.

通量标定沿用同一个有限模板

逐曝光、逐通道取实际卫星斑小图，扣局部背景，对齐后平均形状。模板与卫星斑亮度使用同一个半径三像素的圆孔径和局部背景约定归一，归一量是积分流量，不拿 SATF 头中的峰值顶替。

论文第二节采用 J 波段卫星斑比恒星暗九点三六星等，误差零点一五星等。脚本在选定约简通道间采用这个常数公共比值，没有测得精密的色标定；来源与误差未明确的离轴透射表也未应用。零点一五星等约对应百分之十四的整体尺度误差，相关地作用于各波段，与经验噪声协方差分别保存。

得到的是采用经验模板标定后的约简通道伴星恒星对比度。要换成伴星绝对流量，还需要恒星谱形与标定；原论文采用恒星模板，教程在对比度处停止。

对应阅读：Bonnefoy et al. · Greco & Brandt

处理、注入，再读实际结果

主约简取三张配准图的平均作角差参考，残差各经自己的 CD 矩阵投影再叠加。明确采用历史角距四百三十毫角秒、位置角二百一十一点六度定义局部测量区，已经利用了伴星的先验位置知识，不是盲搜发现。

把经验点源模板逐通道注入原选定数据，整套减除和投影重跑，测全部输出。另跑一次先谱差后角差，展示通道混合。这些确定性注入标定的是冻结观测背景上的处理响应，没有重新生成光子观测。

下表列中心孔径的代表结果；误差是条件经验噪声估计。公共星等标定和各单独卫星模板敏感性分别保存。缺少增益、原始噪声信息时，不能从约简 ADU 精确重建源光子方差。

对应阅读：Bonnefoy et al. · Pueyo · Greco & Brandt

把空间拟合失败放在光谱旁边

十个通道对比度能恰好拟合十个孔径和，因此脚本还逐通道加入中心孔径、径向差和切向差。四十七个比较区估三十维协方差，再在扩大的四十九点位置网格及新增二十四个局部细化点上拟合十个对比度；测量区域和协方差保持不变。

扩大的角差剖面在网格内部得到最小二次残差约一百七十五，对应二十个残差维度；位置约为四百三十毫角秒、二百一十三点一度。谱差加角差的最小残差仍约三百八十九，局部细化也仍碰到边缘。简单空间模型没有拟好这组数据，单独卫星模板也会改变提取颜色。先保留条件诊断光谱，改善 PSF 和背景模型后再拟合大气参数，更符合这些结果。

角差谱差诊断的最后通道还出现约负零点零零零零零八三六的估计，结果原样保留。处理选择和经验模板明显改变光谱，也是误差的一部分，即使合成图像仍有亮点。

输入只有三份公开文件，不是论文完整十九次曝光约简；上游微透镜光谱串扰和曝光内运动也没有完整重建。论文提供历史对照与阅读目标，保存的失败则说明这里实际完成的小规模计算。

对应阅读：Bonnefoy et al. · Samland et al. · Greco & Brandt

在本地运行并保留出处

把下方科学材料下载到同一目录。公开包只含自有脚本、精确字节来源清单、原始头与数值结果；GPI 科学像素和派生图像没有内嵌，数据专用重分发许可仍未明确。下载器直接从原作者固定公开地址读取，用具名 User-Agent 和逐段字节校验。

在干净环境里执行命令。可选 local-images 参数会把派生像素写到本地，供自己查看；让来源清单跟着这些产物走，作图时区分注入标定源与实际观测伴星。

            python3.12 -m venv .venv
. .venv/bin/activate
python -m pip install numpy==2.3.5
python download-gpi-ranges.py --manifest data-source-manifest.json --output gpi-cache
python reproduce-gpi.py --manifest data-source-manifest.json --cache gpi-cache --output local-results
python reproduce-gpi.py --manifest data-source-manifest.json --cache gpi-cache --output local-results --algorithm asdi
python validate-gpi-registration.py --manifest data-source-manifest.json --cache gpi-cache --output registration-validation.json
          

对应阅读：Jason Wang and contributors · GPI Data Pipeline 1.4.0 data formats