*已沦为AI工具人*

好的，老金，判定实验全部备齐，就差你下载 24 个小文件。

已完成

1. 找到理想试验场：西沙 path/row 123/046、123/047——两星都同时有 L1GT（云覆盖致 GCP 匹配失败的日期）和 L1TP（匹配成功日期）。同景网格、地物静态，对比 MTL 的 SCENE_CENTER 直接测出“官方前向模型 vs GCP 平差”的系统差。
2. 排掉一个方法学坑：STAC 目录元数据的中心坐标被产品填充像元污染（东向 ±1.4 km、30 m 量化噪声），测不出来（+120±296 m 不显著）——必须用真 MTL。
3. 产出：
- docs/l1gt_mtl_download_list.txt — 24 条下载直链（2 path/row × L8/L9 × L1GT/L1TP 各 3，全是 2024 年档）
- scripts/check_l1gt_vs_l1tp.py — 一键分析+自动判定（已编译自检通过）

需要你做的

1. 浏览器登录 ers.cr.usgs.gov（EarthExplorer 账号）
2. 打开 docs/l1gt_mtl_download_list.txt，逐条访问 URL（登录后直链生效），每个都是几 KB 的 _MTL.txt
3. 全部放进 D:\bak_miandian\landsat\official_L1TP\
4. 告诉我，我跑：uv run --no-project scripts/check_l1gt_vs_l1tp.py

判定标准（出结果的含义）

结果: L1GT−L1TP 东向 ≈ +200 m
结论: 官方前向模型同样带这个偏差——我链与官方前向完全同级，残余纯属“官方做了逐景
GCP 平差而我们没做”
行动: 加常值 roll 修正（−64″/−68″）对齐官方口径，OLI 残余 →30~60 m
────────────────────────────────────────
结果: L1GT−L1TP ≈ 0~20 m
结论: 官方前向干净，我链确有 −65″ 差
行动: 审计 parse_landsat_ancillary.py 四元数提取口径 vs ADD §4.1.4（逐项对齐
FINE 消息、时标、缩放）

README 已记 v1.0.260917-9。