37,719
社区成员
发帖
与我相关
我的任务
分享
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
from mpl_toolkits import mplot3d
import scipy.interpolate as ip
import matplotlib.mlab as ml
from matplotlib.colors import LightSource
from matplotlib import cm
path = "E:/Model_building/A/"
s1 = pd.read_excel(path+"sheet1.xlsx")
x = np.array(s1.iloc[:,2])
y = np.array(s1.iloc[:,3])
z = np.array(s1.iloc[:,4])
xx = np.linspace(x.min(), x.max())
yy = np.linspace(y.min(), y.max())
X, Y = np.meshgrid(xx, yy)
Z = ip.griddata((x, y), z, (X, Y),method='cubic')
contr = plt.contour(X, Y, Z, 30) # 等高线图
plt.figure()
ax = plt.axes(projection='3d') # 3D地形图
ax.invert_xaxis()
ax.plot_surface(X,Y, Z,rstride=1,cstride=1,cmap="viridis" )
plt.show()