CATTOPY-绘图2D Geodata作为2D经度和纬度的函数

Question

我想从以下NetCDF文件（ R90C_20P_MSG_201403250100.NC ）绘制2D数组 rfnet rfnet 作为TEH 2D数组的函数>纬度。

下载NetCDF文件： 或这里

我尝试了以下操作：

import xarray as xr
data = xr.open_dataset('R90C_20P_MSG_201403250100.nc', engine='netcdf4')
print(data)

哪个给出：

<xarray.Dataset>
Dimensions:    (image_y: 565, image_x: 1215)
Dimensions without coordinates: image_y, image_x
Data variables: (12/15)
    LONGITUDE  (image_y, image_x) float32 ...
    LATITUDE   (image_y, image_x) float32 ...
    TAUCO      (image_y, image_x) float32 ...
    RFLW       (image_y, image_x) float32 ...
    RFSW       (image_y, image_x) float32 ...
    TAUCI      (image_y, image_x) float32 ...
    ...         ...
    OLRCICO    (image_y, image_x) float32 ...
    RSRCICO    (image_y, image_x) float32 ...
    RFnet      (image_y, image_x) float32 ...
    ATD        (image_y, image_x) float32 ...
    FUEL       (image_y, image_x) float32 ...
    PC         (image_y, image_x) float32 ...

我可以提取Eg：我可以提取Eg：

LONGITUDE=data['LONGITUDE']
LATITUDE=data['LATITUDE']
RFnet=data['RFnet']

并使用一个简单的绘图访问相应的数据

LONGITUDE.data
LATITUDE.data
RFnet.data

是：

import matplotlib.pyplot as plt

fig = plt.figure(figsize=(8, 4))

ax = fig.add_subplot(111)
img=plt.imshow(RFnet.data, cmap='seismic', vmin=-10, vmax=10)
ax.set_aspect('equal')

cbar = fig.colorbar(img, location='bottom', shrink=0.6, pad=0.15, ax=ax)
cbar.set_label('RFnet')

plt.show()

​我可以将 rfnet 绘制为经度和 latitude ？

纬度和经度的函数：

Answer 1

fig = plt.figure(figsize=(12, 12))
pc = ccrs.PlateCarree()
ax = fig.add_subplot(projection=pc, extent=(-55.530094, 55.530094, 33.582264, 66.823235))

ax.set_global()
ax.coastlines(resolution='50m', # alternative: '110m', '50m', '10m' 
              color='grey')
ax.gridlines(draw_labels=True) 

x = LONGITUDE.data
y = LATITUDE.data
c = RFnet.data
norm=mcolors.Normalize(vmin=-10, vmax=10)
cmap = "seismic"

ax.set_extent(extents=(-20, 20, 32, 52)) # x0, x1, y0, y1 in degrees

im = ax.scatter(x, y, c=c, cmap=cmap, norm=norm, s=2)

cbar = fig.colorbar(im, location='bottom', shrink=1, pad=0.05, ax=ax)
cbar.set_label('RFnet')

plt.show()