Объект 'QuadContourSet' не имеет атрибута 'set_title'

Question

Я использую контурную карту с участками. Когда я хочу использовать ax3.set_title ('Title'), он возвращает следующую ошибку: объект 'QuadContourSet' не имеет атрибута 'set_title'. Я пытался разными способами, также установить метки заголовка и осей для каждого отдельного графика отдельно, но всегда выдает одну и ту же ошибку ... Как я мог установить метки заголовка и осей? Был бы признателен за любую помощь:)

df = pd.read_csv('energy_2D_map1[enter image description here][1].dat', header = None, sep = "\s+")
x = np.arange(0, 2.6, 0.1)
y = np.arange(0, 4.4, 0.1)
X, Y = np.meshgrid(y, x)
z = df[2]
z_array=z.to_numpy()
Z = np.reshape(z_array, (26, 44))

#Plotting 0leg

fig, (ax0, ax1, ax2, ax3) = plt.subplots(nrows=1, ncols=4,figsize=(9.6,4.4), sharey=True)
subplots_adjust(left=0.16, bottom=0.08, right=1, top=0.97, wspace=0.05,
                hspace=1)
ax0 = ax0.contourf(Y, X, Z, 20, cmap='RdGy',vmin=0,vmax=0.017);



#1leg
df = pd.read_csv('energy_2D_map2.dat', header = None, sep = "\s+")
x = np.arange(0, 2.6, 0.1)
y = np.arange(0, 4.4, 0.1)
X, Y = np.meshgrid(y, x)
z = df[2]
z_array=z.to_numpy()
Z = np.reshape(z_array, (26, 44))

#Plotting 1leg

ax1 = ax1.contourf(Y, X, Z, 20, cmap='RdGy',vmin=0,vmax=0.017);

#2leg

df = pd.read_csv('energy_2D_map3.dat', header = None, sep = "\s+")
x = np.arange(0, 2.6, 0.1)
y = np.arange(0, 4.4, 0.1)
X, Y = np.meshgrid(y, x)
z = df[2]
z_array=z.to_numpy()
Z = np.reshape(z_array, (26, 44))

#Plotting 2leg

ax2 = ax2.contourf(Y, X, Z, 20, cmap='RdGy',vmin=0,vmax=0.017);

#flat

df = pd.read_csv('energy_2D_map4.dat', header = None, sep = "\s+")

x = np.arange(0, 2.6, 0.1)
y = np.arange(0, 4.4, 0.1)
X, Y = np.meshgrid(y, x)

z = df[2]
z_array=z.to_numpy()
Z = np.reshape(z_array, (26, 44))

#Plotting flat

ax3 = ax3.contourf(Y, X, Z, 20, cmap='RdGy',vmin=0,vmax=0.017);
ax3.set_title('Title')

#Cbar settings
cbar_ax = fig.add_axes([1.02, 0.08, 0.05, 0.89])
cb = plt.colorbar(ax0,cax=cbar_ax)
tick_locator = ticker.MaxNLocator(nbins=10)
cb.locator = tick_locator
plt.savefig('graphene_pes.pdf', transparent=True, dpi=100,bbox_inches='tight')
plt.show()

Это то, что я имею в качестве вывода