我将一些观察研究视为“每日”。其中有一个计划
import numpy as np
import pylab as pl
import math
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
x = np.array([-93.690, -97.839, 105.193, -18.550, 166.095, -67.075]) # longitude
y = np.array([14.591, 2.011, -7.380, 68.075, -10.757, -24.041]) # latitude
degtorad = math.pi/180.
alpha = np.arange(-180,180.,1.)
alpha *= degtorad
x1 = np.sin(192.25*degtorad - alpha)
x2 = np.cos(192.25*degtorad - alpha)*np.sin(27.4*degtorad)
yy = np.arctan2(x1, x2)
longitude = 303*degtorad - yy
x3 = np.cos(27.4*degtorad) * np.cos(192.25*degtorad - alpha)
latitude = np.arcsin(x3)
xrad = x * degtorad
yrad = y * degtorad
fig1 = pl.figure(1)
ax1 = fig1.add_subplot(111, projection="mollweide")
ax1.scatter(xrad,yrad)
ax1.plot([-math.pi, math.pi], [0,0],'r-')
ax1.plot([0,0],[-math.pi, math.pi], 'r-')
pl.title("before 27 Feb 12:00-after 28 Feb 24:00")
pl.grid(True)
pl.show()
enter code here
plt.figure(figsize=(6, 3))
ax = plt.axes(projection=ccrs.Mollweide())
ax.coastlines(resolution='110m')
ax.gridlines()
plt.show()
将图形分开进行。我想知道注释表示点的值，并得出两个图形的等效图形结论。enter image description here enter image description here
