带有matplotlib的netCDF数据的怪异图

Question

我一直在寻找可悲的长时间，因此，我将不胜感激。

我正在尝试在南极海中绘制一些海冰干舷数据（netCDF，网格总干舷），但是应该很好地绘制在南极周围的数据位于图像的底部。 NetCDF和matplotlib对我来说还很新，因此错误可能出在例如处理尺寸或投影时。

from scipy.io.netcdf import netcdf_file as Dataset
import numpy as np
import matplotlib.pyplot as plt

FB = Dataset('./datasets/fb-0217-0320.nc', 'r')
f = FB.variables['f'][:,:]
lat = FB.variables['lat'][:,0]
lon = FB.variables['lon'][0,:]
masked_fb = np.ma.masked_where(np.isnan(f), f)
mtx_lon, mtx_lat = np.meshgrid(lon, lat)
m = Basemap(projection='spstere',boundinglat=-50, lon_0=180., resolution='l')
m.bluemarble()

plt.figure()
m.pcolormesh(mtx_lon, mtx_lat, masked_fb, latlon=True)
plt.show()

ncdump给出：

dimensions:
x = 79 ;
y = 83 ;
variables:
float lat(y, x) ;
    lat:standard_name = "latitude" ;
    lat:long_name = "latitude coordinate" ;
    lat:units = "degrees_north" ;
float lon(y, x) ;
    lon:standard_name = "longitude" ;
    lon:long_name = "longitude coordinate" ;
    lon:units = "degrees_east" ;
float f(y, x) ;
    f:long_name = "total_freeboard" ;
    f:units = "mm" ;
    f:coordinates = "lat lon" ;

我注意到的一件奇怪的事是最小纬度为-5156.6201，但我不知道如何计算其中的多少...

编辑：如尼尔建议的那样，将代码格式化为适合通用的方式。

Answer 1

好的，我从matplotlib获得帮助，并认为如果其他人有时也遇到类似的问题，我应该在这里分享。 问题出在网格上。 由于netCDF文件中的纬度和经度已经是2D格式，因此不需要网格。 对我有用的解决方案是：

from scipy.io.netcdf import netcdf_file as Dataset
import numpy as np
import matplotlib.pyplot as plt

FB = Dataset('./datasets/fb-0217-0320.nc', 'r')
f = FB.variables['f'][:,:]
lat = FB.variables['lat'][:,:]
lon = FB.variables['lon'][:,:]
masked_fb = np.ma.masked_where(np.isnan(f), f)
m = Basemap(projection='spstere',boundinglat=-50, lon_0=180., resolution='l')
m.bluemarble()

plt.figure()
m.pcolormesh(lon, lat, masked_fb, latlon=True)
plt.show()

Answer 2

首先，通常的做法是将netcdf模块读为

from scipy.io.netcdf import netcdf_file as Dataset

然后，您可以读取文件并以以下方式访问变量

FB = Dataset('./datasets/fb-0217-0320.nc', 'r')
f = FB.variables['f'][:,:]
lat = FB.variables['lat'][:,:]
lon = FB.variables['lon'][:,:]

您确定lat[:,0]和lon[0,:]在网格坐标中正确读取吗？ ncdump表示它们是2D变量，我怀疑问题是根据lat[:,0]和lon[0,:]创建了一个网meshgrid 。