我正试图在潮汐流域的流量上使用模型输出。该模型使用曲线网格。我的第一个任务是绘制最高水层速度的一个分量。我根据名称下的问题编写了一些代码：Matplotlib Streamplot for Unevenly（curvilinear）Grid。
现在据我所知，除了与早先提到的问题相比的数字之外，我没有改变任何必要的东西，但数字仍然是空的。我把代码和一些数字放在下面。
import numpy as np
from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
Lat = np.array([[ 30.40098833, 30.40103752, 30.40108727, 30.40113704],
[ 30.40140046, 30.40145021, 30.40149997, 30.40154973],
[ 30.40186559, 30.40191478, 30.40196453, 30.4020143 ],
[ 30.40239781, 30.402447, 30.40249676, 30.40254652]])
Lon = np.array([[-86.51729818, -86.51794126, -86.5185871, -86.51923603],
[-86.51725858, -86.51790149, -86.51854717, -86.51919595],
[-86.51721383, -86.51785659, -86.51850228, -86.51915089],
[-86.51716242, -86.51780518, -86.51845087, -86.51909948]])
Xvel = np.array([[ 0.0325774, -0.02811189, -0.04972513, -0.07736091],
[ 0.00592685, -0.00043959, -0.00735147, -0.05015078],
[-0.03365543, -0.03183309, -0.03701356, -0.07232581],
[-0.09578606, -0.10139448, -0.11220678, -0.13221299]])
plt.ion()
fig,(ax1) = plt.subplots(1,1)
m = Basemap(llcrnrlon=Lon.min(),llcrnrlat=Lat.min(),
urcrnrlon=Lon.max(), urcrnrlat=Lat.max(),
projection='merc',resolution='i',ax=ax1)
m.contourf(Lat,Lon,Xvel,latlon=True)
m.drawcoastlines()
m.drawrivers()
m.plot(Lat,Lon,'-k',alpha=0.3,latlon=True)
m.plot(Lat.T,Lon.T,'-k',alpha=0.3,latlon=True)
有人能告诉我是什么导致情节保持空白吗？
我还有另一个关于Basemap使用的问题：我的数据表还包含很多NaN（没有信息的网格点）。我想知道如何让Basemap知道我没有关于这些位置的任何信息，而且我不想在那里进行任何绘图。在当前的代码中，它导致“LinearRing的点”不形成闭合的线串＆＃39;错误。