pygmt.clib.Session.extract_region

Session.extract_region()[source]

Extract the WESN bounding box of the currently active figure.

Retrieves the information from the PostScript file, so it works for country codes as well.

Returns: * wesn (1d array) – A 1D numpy array with the west, east, south, and north dimensions of the current figure.

Examples

>>> import pygmt
>>> fig = pygmt.Figure()
>>> fig.coast(
...     region=[0, 10, -20, -10],
...     projection="M6i",
...     frame=True,
...     land="black",
... )
>>> with Session() as lib:
...     wesn = lib.extract_region()
...
>>> print(", ".join(["{:.2f}".format(x) for x in wesn]))
0.00, 10.00, -20.00, -10.00

Using ISO country codes for the regions (for example 'US.HI' for Hawaii):

>>> fig = pygmt.Figure()
>>> fig.coast(
...     region="US.HI", projection="M6i", frame=True, land="black"
... )
>>> with Session() as lib:
...     wesn = lib.extract_region()
...
>>> print(", ".join(["{:.2f}".format(x) for x in wesn]))
-164.71, -154.81, 18.91, 23.58

The country codes can have an extra argument that rounds the region a multiple of the argument (for example, 'US.HI+r5' will round the region to multiples of 5):

>>> fig = pygmt.Figure()
>>> fig.coast(
...     region="US.HI+r5", projection="M6i", frame=True, land="black"
... )
>>> with Session() as lib:
...     wesn = lib.extract_region()
...
>>> print(", ".join(["{:.2f}".format(x) for x in wesn]))
-165.00, -150.00, 15.00, 25.00