Rotated rectangle using azimuth

The pygmt.Figure.plot method can plot rotated rectangles based on a given azimuth (in degrees east of north) as well as length and width. We can define the required parameters in a numpy array or use an appropriately formatted input file. Such representations are often used to e.g. display results of shear-wave splitting analysis.

Out:

<IPython.core.display.Image object>

import pygmt
import numpy as np

fig = pygmt.Figure()

# generate a basemap around Big Island (Hawai'i) showing coastlines, land, and water
fig.coast(
    region = [-156.5, -154.5, 18.5, 20.5],
    projection = "M6c",
    land = "grey",
    water = "lightblue",
    shorelines = True,
    resolution = "f",
    frame=["x1", "y1"]
)

# store parameters for rotated rectangle in a numpy
# array (lon, lat, azimuth in degrees east of north, lenght, width)
data = np.array([[-155.533, 19.757, 45, 60, 5]])

# pass the data to the plotting function in addition to the corresponding
# style shortcut for rotated rectangles ("J") as well as set color and pen for
# the rectangle
fig.plot(data = data,
         style = "J",
         color = "red3",
         pen = "black")

fig.show()