我试图从netcdf数据中绘制一个4D表面图,它有4个维度:5个变量(DU001、DU002...005)的time、lat、long和lev(sample data)。我必须绘制第一个变量DU001与纬度、长度和水平(72个水平)的关系图,这样x轴是纬度,y轴是长度,z轴是水平,DU001将用颜色表示。到目前为止,我已经尝试了下面的代码,但我只得到一个表面在我的阴谋
。
我认为这只是一个层面。如何纠正?
import xarray as xr
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
import numpy as np
from matplotlib import cm
path= 'D:\\DATA\\2015\\test_data'# Open the NetCDF file
data = xr.open_dataset('D:\\DATA\\2015\\test_data\\MERRA2_400.inst3_3d_aer_Nv.20150515.SUB.nc')
# Select the DU01 variable and the lat, long, and lev dimensions
lat = data['lat']
lon = data['lon']
lev = data['lev']
DMR = data['DU001']
# Reshape the data
du001_2d = DMR[:, :, :].squeeze()
dmr_values = du001_2d.values.squeeze()
# Create meshgrid for coordinates
lon_2d, lat_2d = np.meshgrid(lon, lat)
# Create the 3D plot
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
# Plot the surface for each level
for i, level in enumerate(lev):
ax.plot_surface(lon_2d, lat_2d, dmr_values[i], cmap='viridis')
# Set labels and title
ax.set_xlabel('Longitude')
ax.set_ylabel('Latitude')
ax.set_zlabel('Level')
ax.set_title('3D Plot of Dust Mixing Ratio')
# Set the z-limits based on the valid range of the lev array
ax.set_zlim(lev[0], lev[71]) # Assuming lev is a 1D array
# Display the plot
plt.show()
我不知道我哪里做错了。我对Python还是个新手。任何帮助都将不胜感激
我尝试绘制的示例是
1条答案
按热度按时间guykilcj1#
你把所有的
dmr_values
一个标在另一个上,即你看到的只是dmr_values[-1]
。层的最大高度约为6e-8
,这就是为什么当使用0 - 72的z轴范围时,它只显示为一个平面。如果你想绘制72个彩色图层,你需要提供72个平面并自己给它们上色:
**更新以下注解:**如果要进行插值,可以为3个可见表面使用3个填充等高线图: