我正在使用python读取netcdf文件,需要使用python中的matplotlib库绘制图形。netcdf文件包含3个变量:u
v
和w
组件。我必须在垂直比例上绘制这三个组件。由于这些数据将用于天气预报雷达。我需要在拐角处绘制带有平滑曲线的图形。现在情节看起来并且代码是...
from netCDF4 import Dataset
from pylab import *
import numpy
from scipy import interpolate
from scipy.interpolate import spline
#passing the filename
root_grp=Dataset('C:\\Python27\\MyPrograms\\nnt206rwpuvw.nc')
#getting values of u component
temp1 = root_grp.variables['ucomponent']
data1 = temp1[:]
u=data1[0]
#getting values of v component
temp2 = root_grp.variables['wcomponent']
data2 = temp2[:]
v=data2[0]
#getting values of w component
temp3 = root_grp.variables['wcomponent']
data3 = temp3[:]
w=data3[0]
#creating a new array
array=0.15*numpy.arange(1,55).reshape(1,54)
#aliasing
y=array[0]
#sub-plots
f, (ax1, ax2, ax3) = plt.subplots(1,3, sharey=True)
ax1.plot(u,y,'r')
ax2.plot(v,y,'g')
ax3.plot(w,y,'b')
#texts
ax1.set_title('u component',fontsize=16,color='r')
ax2.set_title('v component',fontsize=16,color='g')
ax3.set_title('w component',fontsize=16,color='b')
show()
一种简单的解决方案是用样条线进行插值,如nice Scipy tutorial中所述。
这里是一个例子:
<< img src =“ https://image.soinside.com/eyJ1cmwiOiAiaHR0cHM6Ly9pLnN0YWNrLmltZ3VyLmNvbS9LS0hDMi5wbmcifQ==” alt =“在此处输入图像描述”>