如何解决访问python底图或cartopy中的底层投影函数
我想使用底图或 cartopy 来绘制卫星绕地球的轨道。这两个库都提供了绘制从太空中的观察者(例如正交或地球静止等)看到的地球的能力,并在给定纬度和经度的情况下绘制地球上的某些东西。鉴于此功能,在给定纬度和经度的情况下,在图中生成一个点的基础投影函数似乎应该是可调用的,其中可以指定不同的高度(海平面除外)。
因此,在给定经度、纬度和高度的情况下,难道不能使用底图或 Cartopy 绘制卫星轨道吗?
注意 #1:我知道可以绘制卫星的地面轨迹,但这忽略了高度,我想看看高度依赖性。
注意 #2:我已经看到您可以将 Basemap 添加到 Matplotlib 3D 轴,但这似乎不适用于“正交”投影。
注意 #3 我正在寻找使用(纬度、经度、高度)坐标制作类似于 3D image of Earth with Satellite Trajectories 的东西的能力
解决方法
您不需要地图投影即可获得所需的环绕地球的卫星图。 3D 绘图可用于创建地球,在其周围叠加地球图像,然后在地球表面上方绘制点。
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
import numpy as np
from matplotlib.cbook import get_sample_data
from matplotlib._png import read_png
# Use world image with shape (360 rows,720 columns)
pngfile = r"./eo_base_2020_clean_720x360.png"
fn = get_sample_data(pngfile,asfileobj=False)
img = read_png(fn) # get array of color
# Some needed functions / constant
r = 5
pi = np.pi
cos = np.cos
sin = np.sin
sqrt = np.sqrt
# Prep values to match the image shape (360 rows,720 columns)
phi,theta = np.mgrid[0:pi:360j,0:2*pi:720j]
# Parametric eq for a spherical globe
# phi=latitude; theta=longitude
x = r * sin(phi) * cos(theta)
y = r * sin(phi) * sin(theta)
z = r * cos(phi)
fig = plt.figure()
fig.set_size_inches(9,9)
ax = fig.add_subplot(111,projection='3d',label='axes1')
# drape the image on the globe's surface
sp = ax.plot_surface(x,y,z,\
rstride=2,cstride=2,\
facecolors=img)
# plot some points above/around the globe
phi,theta = np.mgrid[0:pi:50j,0:2*pi:50j]
r2 = 1.5*r
x = r2 * sin(phi) * cos(theta)
y = r2 * sin(phi) * sin(theta)
z = r2 * cos(phi) + 0.8* sin(sqrt(x**2 + y**2)) * cos(2*theta)
ax.scatter(x,s=1,color='red')
ax.set_axis_off()
ax.set_aspect(1)
plt.show()
地图投影是坐标变换
所有地图投影都是从地理定位到平面坐标系的一组转换函数。实际实施需要一定的规则和范围,以避免误用。我上面的代码可以被认为是一个地方的 2 个地图投影。第一个地球使用 r=5 作为地球半径,该投影用于在投影平面上绘制 bluemarble 图像的像素。在公共位置的另一个地球仪使用半径 r2 = 1.5*r。它代表另一个投影,用于绘制可以是轨道地球卫星的点。
在一个 axis
或 Cartopy
的 Basemap
上实现 2 个叠加地图投影在技术上是可行的,在绘图过程中可以使用一些特殊的坐标转换。
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。