可以用xarray库代替netCDF4吗

是的，您可以使用xarray库代替netCDF4库来读取和处理nc文件。xarray库是基于pandas的Python库，用于处理标记数据集的N维数组。

下面是使用xarray库绘制气温等值线图的示例代码：

import xarray as xr
import numpy as np
import matplotlib.pyplot as plt

# 读取nc文件
data = xr.open_dataset('stipple.nc')
temp = data['temp'].sel(time='1963-01')

# 获取经度和纬度
lon = data['lon']
lat = data['lat']

# 绘制等值线图
plt.contour(lon, lat, temp, levels=10, colors='black')

# 标记正负距平格点
pos = np.where(temp > 0)
plt.scatter(lon[pos], lat[pos], marker='+', color='red')
neg = np.where(temp < 0)
plt.scatter(lon[neg], lat[neg], marker='-', color='blue')

plt.show()

这个代码示例使用了xarray库的open_dataset函数打开nc文件，然后使用了sel函数选择1963年1月份的数据。经度和纬度是通过data变量的索引获取的。

在绘制等值线图时，我们仍然使用了matplotlib中的contour函数，它可以根据温度值绘制等值线。我们选择了10个等级，从而将等值线分成了10个区域。标记正负距平格点时，我们使用了matplotlib中的scatter函数，分别绘制了红色的"+"和蓝色的"-"。