如何用Python计算卫星3D简单多边形几何体的阴影投射

问题背景 #

需构建Python库计算卫星太阳能板的受阴影面积（用于发电量估算），卫星结构为立方体+侧面突出矩形面板，太阳建模为三维点光源。核心需求是分析特定太阳能板多边形上的投射阴影区域，而非仅可视化阴影效果。

现有库的局限 #

调研PyVista、PyOpenGL、trimesh、pyrender后，发现两类核心问题：

• 侧重可视化而非分析：仅能渲染阴影，无法提取阴影在目标多边形上的精确几何区域
• 3D多边形分析支持不足：缺乏直接计算阴影多边形与目标面重叠区域的原生能力

自定义实现优化方案 #

如果没有原生支持阴影投射分析的库，可以基于基础几何计算工具组合实现，优化你提出的步骤如下：

步骤1：定义卫星部件的多边形 #

用三维坐标数组定义所有卫星部件的面（立方体6个面+侧面矩形面板），每个面存储为按顺时针/逆时针顺序排列的顶点列表（保证面法线方向正确）。

步骤2：点光源建模 #

将太阳位置存储为(x,y,z)三维坐标点，作为射线发射的原点。

步骤3：目标面的阴影计算（优化版） #

针对每个待分析的太阳能板面：

1. 构建目标面平面方程：通过目标面的三个顶点计算平面法向量，推导平面方程ax + by + cz + d = 0
2. 生成阴影轮廓射线：无需向其他面的每个点发射射线，仅向遮挡面的顶点发射即可——点光源下，多边形的阴影轮廓由顶点的投射射线决定
3. 计算射线与目标平面的交点：对每条从光源到遮挡面顶点的射线，求解与目标平面的交点，得到阴影多边形的顶点集合
4. 构建闭合阴影多边形：将交点按遮挡面顶点的顺序连接，剔除在目标面范围外的交点，形成闭合多边形
5. 计算重叠区域：将目标面和阴影多边形投影到二维平面（比如转换为目标面的局部坐标系），用几何计算库的多边形交集功能计算重叠面积

推荐工具库补充 #

• numpy：用于三维向量、平面方程、射线交点的数值计算
• shapely：高效处理二维多边形的交集、面积计算（三维问题降维后使用）
• pygeos：可选，比shapely更快的几何计算库，适合批量处理

内容的提问来源于stack exchange，提问作者Sander Voss