STK绘图

tle数据来源于[celestrak](https://celestrak.org)网站 # Iridium next - [tle数据链接](https://celestrak.org/NORAD/elements/gp.php?GROUP=iridium-NEXT&FORMAT=2le)    # Qianfan - [tle数据链接](https://celestrak.org/NORAD/elements/gp.php?GROUP=qianfan&FORMAT=2le)   # Oneweb - [tle数据链接](https://celestrak.org/NORAD/elements/gp.php?GROUP=oneweb&FORMAT=2le)   # Starlink-shell-53-part - [tle数据链接](https://celestrak.org/NORAD/elements/gp.php?GROUP=starlink&FORMAT=2le) 网站上Starlink的tle文件包含了所有轨道的卫星，使用下面的脚本筛选出轨道倾角为53度的卫星，有4901颗，如果全部进行模拟，不便于观察，因此只模拟了部分卫星 ```python import requests # 1. 从 CelesTrak 获取所有 Starlink TLE url = "https://celestrak.org/NORAD/elements/supplemental/sup-gp.php?FILE=starlink&FORMAT=tle" response = requests.get(url) lines = response.text.strip().split('\n') target_tle = [] # 2. 遍历 TLE（每 3 行一组：名称、Line 1、Line 2） for i in range(0, len(lines), 3): name = lines[i].strip() l1 = lines[i + 1] l2 = lines[i + 2] # 提取倾角：Line 2 的 [8:16] 位置 try: inclination = float(l2[8:16].strip()) # 筛选 53° 附近的卫星（考虑到摄动，范围设为 52.5-54.0） if 52.5 <= inclination <= 54.0: target_tle.extend([name, l1, l2]) except ValueError: continue # 3. 保存结果 with open("starlink_53_shell.tle", "w") as f: f.write('\n'.join(target_tle)) print(f"筛选完成，共提取 {len(target_tle) // 3} 颗卫星。") ```