Python-сообщество

from mpl_toolkits.mplot3d import axes3d
import matplotlib.pyplot as plt
import numpy as np
from matplotlib import cm
import scipy.signal
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
xydata, zdata = np.loadtxt('zdata1.txt', unpack = True)
Z1, Z2 = np.meshgrid(zdata, zdata)
Z_1 = Z1 * Z2
xydata, zdata = np.loadtxt('zdata2.txt', unpack = True)
Z1, Z2 = np.meshgrid(zdata, zdata)
Z_2 = Z1 * Z2
xydata, zdata = np.loadtxt('zdata3.txt', unpack = True)
Z1, Z2 = np.meshgrid(zdata, zdata)
Z_3 = Z1 * Z2
xydata, zdata = np.loadtxt('zdata4.txt', unpack = True)
Z1, Z2 = np.meshgrid(zdata, zdata)
Z_4 = Z1 * Z2
xydata, zdata = np.loadtxt('zdata5.txt', unpack = True)
Z1, Z2 = np.meshgrid(zdata, zdata)
Z_5 = Z1 * Z2
Z = Z_1 + Z_2 + Z_3 + Z_4 + Z_5
X, Y = np.meshgrid(xydata, xydata)
t = 1 - np.abs(np.linspace(-1, 1, 21))
kernel = t.reshape(21, 1) * t.reshape(1, 21)
kernel /= kernel.sum()
R = scipy.signal.convolve2d(Z[:,:], kernel, mode='same')
R = R / np.max(zdata)
ax.plot_surface(
    X, Y, R, rstride=6, cstride=6,
    cmap=cm.jet,
    linewidth=0, antialiased=False, shade=False)
plt.show()