1

绘制随机取向的三轴椭球与极坐标平面的交点的最佳方法是什么？这架飞机在给定纬度的一系列经度上运行。

下面的代码应绘制球体数组与纬度为 plane_lat 的平面的交点。（球体的中心坐标：sphere_x,sphere_y,sphere_z；距原点的距离：sphere_dist；球半径：sphere_rad。）

for i in range(len(hole_rad)):
    deltaz = (sphere_dist[i]*np.cos(sphere_lat[i]*degtorad))*np.tan(plane_lat*degtorad)-sphere_dist[i]*np.sin(sphere_lat[i]*degtorad)
    if np.abs(deltaz)<sphere_radius[i]:
        rprime = sphere_rad[i]*np.sin(np.arccos(abs(deltaz)/(sphere_rad[i])));
        x = rprime * np.sin(newtheta)+sphere_x[i]*H0
        y = rprime * np.cos(newtheta)+sphere_y[i]*H0
        z = np.zeros(np.shape(newtheta))
        cr,clat,clon=ACD.cartesian_to_spherical(x,y,z)
        circles=ax.plot(np.rad2deg(clon),cr,c='blue',linewidth=0.1)

这是我使用 python/matplotlib 的迂回尝试。必须有更好的方法来实现这一点。关于如何为椭圆体（最好在 python 中）执行此操作的任何想法？

4

0 回答 0