国家太空安全是国家安全在空间领域的表现。随着太空技术在政治、经济、军事、文化等各个领域的应用不断增加，太空已经成为国家赖以生存与发展的命脉之一，凝聚着巨大的国家利益，太空安全的重要性日益凸显^[1]。而在信息化时代，太空安全与信息安全紧密地结合在一起。

2020年9月4日,美国白宫发布了首份针对太空网络空间安全的指令——《航天政策第5号令》，其为美国首个关于卫星和相关系统网络安全的综合性政策，标志着美国对太空网络安全的重视程度达到新的高度。在此背景下，美国自2020年起，连续两年举办太空信息安全大赛“黑掉卫星（Hack-A-Sat）”，在《Hack-A-Sat太空信息安全挑战赛深度解析》一书中有详细介绍，本文介绍了Hack-A-Sat黑掉卫星挑战赛的定位卫星Jackson这道赛题的解题过程。

题目介绍

Let's start with an easy one, I tell you where I'm looking at a satellite, you tell me where to look for it later.

主办方告诉参赛者在哪里看到了一颗卫星，需要参赛者告诉主办方在哪里还可以看到这个卫星。给出的资料有：

（1）压缩包stations.zip，其中文件就是一个stations.txt文件，是TLE文件，关于TLE文件的格式说明在前文已有介绍，为了便于读者阅读，本节会再次给出简要介绍。

（2）给出了一个链接地址，使用netcat连接到题目给的链接后，会给出进一步提示，如图7-1所示（其中的坐标是随机的，时间也是随机的）。

图7-1  jackson题目的提示信息

连接后，会告诉参赛者当前看到这个卫星的时刻、卫星的地心惯性坐标系（ECI）坐标，接着会依次给出3个新的时刻，要求参赛者给出在哪里还可以看到这个卫星，输入具体的经纬度坐标。3次都输入正确后，会给出flag值。

编译及测试

这道挑战题的代码位于jackson目录下，查看challenge、solver目录下的Dockerfile，发现其中用到的是python:3.7-slim，为了加快题目的编译进度，在jackson目录下新建一个文件sources.list，内容如下：

deb https://mirrors.aliyun.com/debian/ bullseye main non-free contrib

deb-src https://mirrors.aliyun.com/debian/ bullseye main non-free contrib

deb https://mirrors.aliyun.com/debian-security/ bullseye-security main

deb-src https://mirrors.aliyun.com/debian-security/ bullseye-security main

deb https://mirrors.aliyun.com/debian/ bullseye-updates main non-free contrib

deb-src https://mirrors.aliyun.com/debian/ bullseye-updates main non-free contrib

deb https://mirrors.aliyun.com/debian/ bullseye-backports main non-free contrib

deb-src https://mirrors.aliyun.com/debian/ bullseye-backports main non-free contrib

将sources.list复制到jackson、challenge、solver目录下，修改challenge、solver目录下的Dockerfile，在所有的FROM python:3.7-slim下方添加：

ADD sources.list /etc/apt/sources.list

打开终端，进入jackson所在目录，执行命令：

sudo make build

此时如果使用make test命令进行测试，会提示错误，如图7-2所示。

图7-2  执行make test命令时的错误信息

查询run.log，得到如下错误信息：

Traceback (most recent call last):

File "challenge.py", line 22, in <module>

ts = load.timescale()

File "/opt/venv/lib/python3.7/site-packages/skyfield/iokit.py", line 314, in timescale

data = self('deltat.data')

File "/opt/venv/lib/python3.7/site-packages/skyfield/iokit.py", line 203, in __call__

download(url, path, self.verbose)

File "/opt/venv/lib/python3.7/site-packages/skyfield/iokit.py", line 528, in download

raise e2

OSError: cannot get ftp://cddis.nasa.gov/products/iers/deltat.databecause <urlopen error ftp error: TimeoutError(110, 'Connection timed out')>

Try opening the same URL in your browser to learn more about the problem.

If you want to fall back on the timescale files that Skyfield ships with,

try `.timescale(builtin=True)` instead.

错误原因是需要到NASA的CDDIS的ftp下载其中的文件，但是无法打开该链接。CDDIS（Crustal Dynamics Data Information System，地壳动力学数据信息系统）最初是为NASA的地壳动力学项目（Crustal Dynamics Project，CDP）提供中央数据库而开发的，建立于1982年，是一个专用数据库，用于归档和分发与空间大地测量相关的数据集。

CDDIS主要归档和分发如下数据：

• 全球导航卫星系统GNSS的广播星历和精密星历：包括美国的GPS、俄罗斯的GLONASS、中国的北斗等。
• 激光测距：包括人造卫星激光测距和月球激光测距。
• 甚长基线干涉测量（Very Long Baseline Interferometry，VLBI）。
• 星基多普勒轨道确定和无线电定位组合系统（Doppler Orbitography and Radio-positioning Integrated by Satellite，DORIS）。

但是，从2020年10月31日起，因美国政府安全要求不再允许CDDIS通过传统的未加密匿名ftp提供数据，所有数据仍然可用，但是必须通过HTTPS或ftp-ssl进行访问，所以上述代码会报错。本节为了简化，直接修改challenge、solver两个python文件中的所有：

load.timescale()

将其改为

load.timescale(builtin=True)

这里使用的是Python的Skyfield库，不带参数时，将从上述ftp地址下载国际地球自转服务（International Earth Rotation Service，IERS），这个服务的内容很多，其中一项是世界时，参考run.log中的错误提示，这里将其参数改为内置的，就表示不再从NASA的ftp上读取数据。

再次使用make test命令进行测试，会顺利通过，输出信息如图7-3所示。

图7-3  jackson挑战题测试输出

相关背景知识

1．卫星星历TLE文件介绍

TLE是两行轨道根数（TLE），覆盖了气象卫星、海洋卫星、地球资源卫星、教育卫星等应用卫星。以北斗的某颗卫星TLE数据为例，如下：

BEIDOU 2A

1 30323U 07003A  07067.68277059  .00069181  13771-5  44016-2 0  587

2 30323 025.0330 358.9828 7594216 197.8808 102.7839 01.92847527  650

第一行主要元素解析如下：

（1）30323U：30323是北美防空司令部给出的卫星编号，U代表不保密，我们看到的都是U，否则我们就不会看到这组TLE了。

（2）07003A：国际编号，07表示2007年，003表示这一年的第3次发射，A表示这次发射编号为A的物体，其他还有B、C、D等。国际编号就是2007-003A。

（3）07067.68277059：表示这组轨道数据的时间点，07表示2007年，067表示第67天，即3月8日。

（4）68277059：表示这一天里的时刻，大约是16时22分左右。

（5）58：表示关于这个空间物体的第58组TLE。

（6）7：最后一位是校验位。

第二行主要元素解析如下：

（1）30323：北美防空司令部给出的卫星编号。

（2）025.0330：轨道倾角。

（3）358.9828：升交点赤经。

（4）7594216：轨道偏心率。

（5）197.8808：近地点幅角。

（6）102.7839：平近点角，表示在给出这组TLE时，卫星在轨道的什么位置。

（7）01.92847527：每天环绕地球的圈数。其倒数就是周期。可以看出，该北斗卫星目前的周期大约是12h。

（8）65：发射以来飞行的圈数。

（9）0：校验位。

2．地心惯性坐标系ECI介绍

地心惯性坐标系（Earth Center Inertial Coordinates，ECI），原点是地球质心，z轴是地球平均自转极点，x轴是春分点（每年春分点均会发生变动，参考J2000.0），y轴由右手系决定。

题目解析

这道题目的解法还是比较直观的，使用Python提供的Skyfield、NumPy库，可以分为两步：

（1）已知在某个时刻目标卫星的ECI坐标，依据此信息，从给出的station文件中找到目标卫星对应的TLE。

（2）已知目标卫星的TLE，那么就可以计算任意时刻的ECI坐标。

关键代码如下：

from pwn import *

import numpy as np

from skyfield.api import load

import astropy.units

# 加载TLE文件，读出所有的卫星信息，保存在satellites 中

satellites = load.tle_file('./stations.txt')

......

# 下面代码中的t就是题目中给出的观察到目标卫星的那个时刻；eci_coords就是题目中给出的t时刻目标
# 卫星的ECI坐标；通过遍历给出的TLE，取出在t时刻与给定坐标最接近的卫星

match = satellites[np.argmin([np.linalg.norm(s.at(t).position.km-eci_coords) for s in satellites])]

......

# 通过Skyfield可以获取在题目给出的新时刻new_t时的卫星ECI坐标

x,y,z = match.at(new_t).position.km

其中match存储的就是目标卫星的TLE，其计算过程如下：

（1）遍历station中的所有卫星。

（2）对其中的每颗卫星计算其在时刻t的坐标，坐标系是ECI，单位是km。

（3）将上一步得到的坐标与目标卫星在t时刻的坐标相减，然后调用np.linalg.norm函数计算结果的范数，默认就是x、y、z轴坐标差值的平方和再开根号。假设当前从station中取出的卫星坐标是（x,y,x），目标卫星坐标是（x₀,y₀,z₀），那么实际计算的就是如下：

（4）将station中所有卫星进行上述运算，取出范数最小的卫星，这个卫星就是目标卫星。

（5）知道了目标卫星的TLE，就可以通过Skyfield可以获取在题目给出的新时刻new_t时目标卫星的ECI坐标，将该坐标输入终端即可。