Existe alguma maneira de fazer gráficos aplpy de multiplot compartilharem dinamicamente eixos para que, quando um é movido ou ampliado, ele se move e amplia os outros?

Posso conseguir o efeito usando as rotinas de imshow e subplot do matplotlib pyplot, mas usando esses limites alguns outros aspectos importantes da minha plotagem, enquanto o aplpy fornece todas as ferramentas necessárias para minhas imagens.

Tentei usar os comandos matplotlib cid e uma função para atualizar todas as imagens com base nas localizações dos cliques, mas só posso aumentar ou diminuir o zoom, não as duas, e não posso clicar e arrastar ainda.

Meu MWE do meu código de plotagem está abaixo:

from astropy.io import fits
import matplotlib.pyplot as plt
import aplpy

root = '/my/data/directory/'
data = '3d_image.fits'

hdu = fits.open(root + data)[0]
hdr = hdu.header

fits1 = fits.PrimaryHDU(data = hdu.data[4,:,:], header = hdr)
fits2 = fits.PrimaryHDU(data = hdu.data[6,:,:], header = hdr)

fig = plt.figure(figsize=(15, 15))

f1 = aplpy.FITSFigure(fits1, figure=fig, subplot=[0.1,0.1,0.8,0.35])
f1.show_colorscale(cmap = 'coolwarm', vmin = 8., vmax = 10.5)

f2 = aplpy.FITSFigure(fits2, figure=fig, subplot=[0.1,0.5,0.8,0.35])
f2.show_colorscale(cmap = 'coolwarm', vmin = 1.2, vmax = 1.6)

fig.show