The region parameter to grdimage can set the plot region:
def plot_map(grid, area):
"""Plot the grid map out with a mercator projection."""
fig = pygmt.Figure()
fig.grdimage(grid=grid, projection="M20c", frame="a", region=area, cmap="geo")
fig.show()