Broken axis example: uneven subplot size
matplotlib, python, subplot
Solution
UPDATE 2018
There is a github project https://github.com/bendichter/brokenaxes which probably is more convenient to use.
My own solution looks like (using gridspec, assuming that the units of the two y-axis should be equal):
"""
Broken axis example, where the y-axis will have a portion cut out.
"""
import matplotlib.pylab as plt
import matplotlib.gridspec as gridspec
import numpy as np
pts = np.array([ 0.015, 0.166, 0.133, 0.159, 0.041, 0.024, 0.195,
0.039, 0.161, 0.018, 0.143, 0.056, 0.125, 0.096, 0.094, 0.051,
0.043, 0.021, 0.138, 0.075, 0.109, 0.195, 0.05 , 0.074, 0.079,
0.155, 0.02 , 0.01 , 0.061, 0.008])
pts[[3,14]] += .8
ylim = [0.82, 1.0]
ylim2 = [0.0, 0.32]
ylimratio = (ylim[1]-ylim[0])/(ylim2[1]-ylim2[0]+ylim[1]-ylim[0])
ylim2ratio = (ylim2[1]-ylim2[0])/(ylim2[1]-ylim2[0]+ylim[1]-ylim[0])
gs = gridspec.GridSpec(2, 1, height_ratios=[ylimratio, ylim2ratio])
fig = plt.figure()
ax = fig.add_subplot(gs[0])
ax2 = fig.add_subplot(gs[1])
ax.plot(pts)
ax2.plot(pts)
ax.set_ylim(ylim)
ax2.set_ylim(ylim2)
plt.subplots_adjust(hspace=0.03)
ax.spines['bottom'].set_visible(False)
ax2.spines['top'].set_visible(False)
ax.xaxis.tick_top()
ax.tick_params(labeltop='off')
ax2.xaxis.tick_bottom()
ax2.set_xlabel('xlabel')
ax2.set_ylabel('ylabel')
ax2.yaxis.set_label_coords(0.05, 0.5, transform=fig.transFigure)
kwargs = dict(color='k', clip_on=False)
xlim = ax.get_xlim()
dx = .02*(xlim[1]-xlim[0])
dy = .01*(ylim[1]-ylim[0])/ylimratio
ax.plot((xlim[0]-dx,xlim[0]+dx), (ylim[0]-dy,ylim[0]+dy), **kwargs)
ax.plot((xlim[1]-dx,xlim[1]+dx), (ylim[0]-dy,ylim[0]+dy), **kwargs)
dy = .01*(ylim2[1]-ylim2[0])/ylim2ratio
ax2.plot((xlim[0]-dx,xlim[0]+dx), (ylim2[1]-dy,ylim2[1]+dy), **kwargs)
ax2.plot((xlim[1]-dx,xlim[1]+dx), (ylim2[1]-dy,ylim2[1]+dy), **kwargs)
ax.set_xlim(xlim)
ax2.set_xlim(xlim)
plt.savefig('broken_axis-mod.png')
Problem
I haven't found a solution to adjust the height of the bottom and top plot of the broken axis example of matplotlib. BTW: The space between the two plots can be adjusted by: ``` plt.subplots_adjust(hspace=0.03) ``` UPDATE: I've almost figured it out using gridspec: ``` """ Broken axis example, where the y-axis will have a portion cut out. """ import matplotlib.pylab as plt # NEW: import matplotlib.gridspec as gridspec import numpy as np pts = np.array([ 0.015, 0.166, 0.133, 0.159, 0.041, 0.024, 0.195, 0.039, 0.161, 0.018, 0.143, 0.056, 0.125, 0.096, 0.094, 0.051, 0.043, 0.021, 0.138, 0.075, 0.109, 0.195, 0.05 , 0.074, 0.079, 0.155, 0.02 , 0.01 , 0.061, 0.008]) pts[[3,14]] += .8 # NEW: gs = gridspec.GridSpec(2, 1, height_ratios=[1, 3]) ax = plt.subplot(gs[0]) ax2 = plt.subplot(gs[1]) plt.subplots_adjust(hspace=0.03) ax.plot(pts) ax2.plot(pts) ax.set_ylim(.78,1.) ax2.set_ylim(.0,.22) ax.spines['bottom'].set_visible(False) ax2.spines['top'].set_visible(False) ax.xaxis.tick_top() ax.tick_params(labeltop='off') ax2.xaxis.tick_bottom() d = .03 kwargs = dict(transform=ax.transAxes, color='k', clip_on=False) ax.plot((0-d,0+d), (0-d,0+d), **kwargs) # top-left diagonal ax.plot((1-d,1+d), (0-d,0+d), **kwargs) # top-right diagonal kwargs.update(transform=ax2.transAxes) # switch to the bottom axes ax2.plot((0-d,0+d),(1-d,1+d), **kwargs) # bottom-left diagonal ax2.plot((1-d,1+d),(1-d,1+d), **kwargs) # bottom-right diagonal plt.show() ``` The remaining issues are: The two parallel line segments (marker for the broken axis) are not parallel anymore because of the different length of the y-axis. A further issue is how to position conveniently the ylabel. Finally I don't see the option `sharex=True` in gridspec. Does this matter? UPDATE: added ylim and ylim2 parameters do determine the height ratios in such a way that the data units are equal: ``` ylim = [0.8, 1.0] ylim2 = [0.0, 0.3] ylimratio = (ylim[1]-ylim[0])/(ylim2[1]-ylim2[0]+ylim[1]-ylim[0]) ylim2ratio = (ylim2[1]-ylim2[0])/(ylim2[1]-ylim2[0]+ylim[1]-ylim[0]) gs = gridspec.GridSpec(2, 1, height_ratios=[ylimratio, ylim2ratio]) ax = plt.subplot(gs[0]) ax2 = plt.subplot(gs[1]) ax.set_ylim(ylim) ax2.set_ylim(ylim2) ```