HMI Showcase: Cutout#

This example demonstrates how to plot a cutout region of a Map with connector lines that indicate the region of interest in the full-disk image.

Since this example deals with the creation of a specific style of image, there are multiple lines that deal directly with matplotlib axes.

import matplotlib.colors
import matplotlib.pyplot as plt
from matplotlib.patches import ConnectionPatch

import astropy.units as u
from astropy.coordinates import SkyCoord

from import HMI_LOS_IMAGE

First, we use the sample HMI LOS image and focus the cutout over an active region near the solar center.

magnetogram =
left_corner = SkyCoord(Tx=-142*u.arcsec, Ty=50*u.arcsec, frame=magnetogram.coordinate_frame)
right_corner = SkyCoord(Tx=158*u.arcsec, Ty=350*u.arcsec, frame=magnetogram.coordinate_frame)

We clean up the magnetogram by masking off all data that is beyond the solar limb.

hpc_coords =
mask =
magnetogram_big =, magnetogram.meta, mask=mask)

We create the figure in two stages. The first stage is plotting the full-disk magnetogram.

fig = plt.figure(figsize=(7.2, 4.8))

We create a nice normalization range for the image.

norm = matplotlib.colors.SymLogNorm(50, vmin=-7.5e2, vmax=7.5e2)

Plot the full-disk magnetogram.

ax1 = fig.add_subplot(121, projection=magnetogram_big)
magnetogram_big.plot(axes=ax1, cmap='RdBu_r', norm=norm, annotate=False,)
magnetogram_big.draw_grid(axes=ax1, color='black', alpha=0.25, lw=0.5)
<CoordinatesMap with 2 world coordinates:

  index aliases    type   unit    wrap   format_unit visible
  ----- ------- --------- ---- --------- ----------- -------
      0     lon longitude  deg 180.0 deg         deg     yes
      1     lat  latitude  deg      None         deg     yes


These lines deal with hiding the axis, its ticks and labels.

We draw the rectangle around the region we plan to showcase in the cutout image.

magnetogram_big.draw_quadrangle(left_corner, top_right=right_corner, edgecolor='black', lw=1)
<astropy.visualization.wcsaxes.patches.Quadrangle object at 0x7f57c0944700>

The second stage is plotting the zoomed-in magnetogram.

magnetogram_small = magnetogram.submap(left_corner, top_right=right_corner)
ax2 = fig.add_subplot(122, projection=magnetogram_small)
im = magnetogram_small.plot(axes=ax2, norm=norm, cmap='RdBu_r', annotate=False,)

Unlike the full-disk image, here we just clean up the axis labels and ticks.

Now for the finishing touches, we add two lines that will connect the two images as well as a colorbar.

xpix, ypix = magnetogram_big.wcs.world_to_pixel(right_corner)
con1 = ConnectionPatch(
    (0, 1), (xpix, ypix), 'axes fraction', 'data', axesA=ax2, axesB=ax1,
    arrowstyle='-', color='black', lw=1
xpix, ypix = magnetogram_big.wcs.world_to_pixel(
    SkyCoord(right_corner.Tx, left_corner.Ty, frame=magnetogram_big.coordinate_frame))
con2 = ConnectionPatch(
    (0, 0), (xpix, ypix), 'axes fraction', 'data', axesA=ax2, axesB=ax1,
    arrowstyle='-', color='black', lw=1

pos = ax2.get_position().get_points()
cax = fig.add_axes([
    pos[0, 0], pos[1, 1]+0.01, pos[1, 0]-pos[0, 0], 0.025
cbar = fig.colorbar(im, cax=cax, orientation='horizontal')

For the colorbar we want it to have three fixed ticks.

cbar.locator = matplotlib.ticker.FixedLocator([-1e2, 0, 1e2])
cbar.set_label("LOS Magnetic Field [gauss]", labelpad=-40, rotation=0)
hmi cutout

Total running time of the script: (0 minutes 2.659 seconds)

Gallery generated by Sphinx-Gallery