Drawing a rectangle on a map

This example will demonstrate how to draw a rectangle on a map using draw_quadrangle().

import matplotlib.pyplot as plt

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

import sunpy.data.sample
import sunpy.map

Let’s start with a sample AIA image.

aia_map = sunpy.map.Map(sunpy.data.sample.AIA_171_IMAGE)

Here are four different ways to draw a rectangle. The first three ways directly calls the SkyCoord class. The fourth way converts pixel coordinates to the equivalent SkyCoord objects using the pixel_to_world().

fig = plt.figure(figsize=(5, 5))
fig.add_subplot(111, projection=aia_map)
aia_map.plot(clip_interval=(1, 99.99)*u.percent)

Out:

<matplotlib.image.AxesImage object at 0x7f40baf403d0>

Specify two opposite corners of the rectangle as a single, two-element SkyCoord object.

coords = SkyCoord(
    Tx=(100, 500) * u.arcsec,
    Ty=(200, 500) * u.arcsec,
    frame=aia_map.coordinate_frame,
)
aia_map.draw_quadrangle(
    coords,
    edgecolor="blue",
    linestyle="-",
    linewidth=2,
    label='2-element SkyCoord'
)

Out:

<astropy.visualization.wcsaxes.patches.Quadrangle object at 0x7f40bd09d450>

Specify two opposite corners of the rectangle as separate SkyCoord objects.

bottom_left = SkyCoord(-500 * u.arcsec, 200 * u.arcsec, frame=aia_map.coordinate_frame)
top_right = SkyCoord(-100 * u.arcsec, 500 * u.arcsec, frame=aia_map.coordinate_frame)
aia_map.draw_quadrangle(
    bottom_left,
    top_right=top_right,
    edgecolor="green",
    linestyle="--",
    linewidth=2,
    label='two SkyCoords'
)

Out:

<astropy.visualization.wcsaxes.patches.Quadrangle object at 0x7f40ba950450>

Specify one corner of the rectangle and the rectangle’s width and height.

bottom_left = SkyCoord(-500 * u.arcsec, -500 * u.arcsec, frame=aia_map.coordinate_frame)
width = 400 * u.arcsec
height = 300 * u.arcsec
aia_map.draw_quadrangle(
    bottom_left,
    width=width,
    height=height,
    edgecolor="yellow",
    linestyle="-.",
    linewidth=2,
    label='width/height'
)

Out:

<astropy.visualization.wcsaxes.patches.Quadrangle object at 0x7f40bd07ae10>

Draw a desired rectangle in pixel coordinates by first converting to SkyCoord objects.

bottom_left = aia_map.pixel_to_world(600 * u.pixel, 350 * u.pixel)
top_right = aia_map.pixel_to_world(800 * u.pixel, 450 * u.pixel)
aia_map.draw_quadrangle(
    bottom_left,
    top_right=top_right,
    edgecolor="red",
    linestyle=":",
    linewidth=2,
    label='pixel_to_world()'
)
plt.legend()

plt.show()
AIA $171 \; \mathrm{\mathring{A}}$ 2011-06-07 06:33:02

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

Gallery generated by Sphinx-Gallery