Plotly中如何用go.Scattergeo地球底图绘制天体轨道

你的实现存在两个层面的硬性不兼容，直接导致类型报错和渲染失效：

• 类型嵌套错误：go.Figure是画布容器对象，不是可渲染的trace类型，你把Earth = go.Figure(go.Scattergeo())生成的Figure实例直接传入另一个Figure的data列表，会直接触发类型校验错误。
• 坐标系体系冲突：Scattergeo属于二维地理投影图层，采用经纬度（WGS84椭球）坐标系；Scatter3d属于三维笛卡尔坐标系，二者没有共享的投影变换逻辑，即使解决类型问题，轨道点也会完全错位。

方案1：在3D笛卡尔场景复刻正射纹理地球（推荐） #

你的轨道计算输出本身就是地心直角坐标，在3D场景中构建带真实地表纹理的球面作为底图，和原有3D轨道坐标系完全匹配，观感和Scattergeo的正射地球效果无差异，同时能完整保留轨道的空间三维关系，不会出现投影失真。

具体实现代码：

import plotly.graph_objects as go
import numpy as np
import pandas as pd

# 统一单位：以地球平均半径6371km为长度基准，和你轨道计算的单位保持一致
R_EARTH = 6371

# 生成3D球面网格
theta = np.linspace(0, 2 * np.pi, 100)
phi = np.linspace(0, np.pi, 100)
theta, phi = np.meshgrid(theta, phi)
x_sphere = R_EARTH * np.cos(theta) * np.sin(phi)
y_sphere = R_EARTH * np.sin(theta) * np.sin(phi)
z_sphere = R_EARTH * np.cos(phi)

# 构建带真实地表纹理的地球底图（匹配Scattergeo正射模式的视觉效果）
earth_mesh = go.Surface(
    x=x_sphere, y=y_sphere, z=z_sphere,
    showscale=False,
    colorscale='Earth',
    lighting=dict(ambient=0.8, roughness=0.9),
    opacity=1
)

# 加载你的轨道计算结果
Traj = dict()
Traj['Xc'] = DATA[0]
Traj['Yc'] = DATA[1]
Traj['Zc'] = DATA[2]
Traj['Time'] = DATA[3]
Traj = pd.DataFrame(Traj)

# 构建3D轨道trace
orbit_trace = go.Scatter3d(
    x=Traj["Xc"], y=Traj["Yc"], z=Traj["Zc"],
    marker=dict(size=3),
    line=dict(color='pink', width=1.5)
)

# 组合画布，调整为正射投影观感
fig = go.Figure(data=[earth_mesh, orbit_trace])
fig.update_layout(
    height=600,
    margin={"r": 0, "t": 0, "l": 0, "b": 0},
    scene=dict(
        # 隐藏3D坐标轴，和Scattergeo的简洁显示风格一致
        xaxis=dict(showbackground=False, showticklabels=False, visible=False),
        yaxis=dict(showbackground=False, showticklabels=False, visible=False),
        zaxis=dict(showbackground=False, showticklabels=False, visible=False),
        # 调整相机位置，模拟正射投影的平视观感
        camera=dict(eye=dict(x=1.5, y=1.5, z=1.5))
    )
)
fig.show()

方案2：基于原生Scattergeo图层叠加投影轨道 #

如果必须使用Scattergeo内置的正射地球效果，需要把三维直角坐标转换为对应经纬度，将Scatter3d替换为同体系的Scattergeo trace，才能实现同图层叠加。

该方案仅适用于近地轨道演示，因为Scattergeo没有高度维度，无法还原轨道的3D空间高差

具体实现代码：

# 地心直角坐标转经纬度转换函数
def cart2geo(x, y, z):
    r = np.sqrt(x**2 + y**2 + z**2)
    lat = np.arcsin(z / r) * 180 / np.pi
    lon = np.arctan2(y, x) * 180 / np.pi
    return lat, lon

# 批量转换轨道坐标
lat_arr, lon_arr = cart2geo(Traj["Xc"].values, Traj["Yc"].values, Traj["Zc"].values)

# 构建地理图层的轨道trace
orbit_geo_trace = go.Scattergeo(
    lon=lon_arr,
    lat=lat_arr,
    mode='lines',
    line=dict(color='pink', width=1.5)
)

# 构建纯底图的Scattergeo trace（注意是trace，不是Figure对象）
earth_geo_trace = go.Scattergeo(
    lon=[],
    lat=[]
)

# 同层组合两个trace，无Figure嵌套
fig = go.Figure(data=[earth_geo_trace, orbit_geo_trace])
fig.update_geos(
    projection_type="orthographic",
    showland=True,
    showocean=True,
    oceancolor="rgb(32, 102, 148)"
)
fig.update_layout(height=600, margin={"r":0,"t":0,"l":0,"b":0})
fig.show()

你原代码的核心错误就是把go.Scattergeo()这个trace包装成了go.Figure实例，又把Figure实例塞回另一个Figure的data列表。无论什么场景下，传入Figure data数组的元素都必须是go.Bar/go.Scatter/go.Surface这类trace对象，不能是Figure容器。

内容的提问来源于stack exchange，提问作者Alexsander Costa