def _expand_spatial_bounds_to_fit_axes(bounds, ax_width, ax_height):
b = bounds<EOL>height_meters = util.wgs84_distance(b['<STR_LIT>'], b['<STR_LIT>'], b['<STR_LIT>'], b['<STR_LIT>'])<EOL>width_meters = util.wgs84_distance(b['<STR_LIT>'], b['<STR_LIT>'], b['<STR_LIT>'], b['<STR_LIT>'])<EOL>x_per_y_meters = width_meters / height_meters<EOL>x_per_y_axes = ax_width / ax_height<EOL>if x_per_y_axes > x_per_y_meters:  <EOL><INDENT>width_meters_new = (height_meters * x_per_y_axes)<EOL>d_lon_new = ((b['<STR_LIT>'] - b['<STR_LIT>']) / width_meters) * width_meters_new<EOL>mean_lon = (b['<STR_LIT>'] + b['<STR_LIT>'])/<NUM_LIT><EOL>lon_min = mean_lon - d_lon_new / <NUM_LIT><EOL>lon_max = mean_lon + d_lon_new / <NUM_LIT><EOL>spatial_bounds = {<EOL>"<STR_LIT>": lon_min,<EOL>"<STR_LIT>": lon_max,<EOL>"<STR_LIT>": b['<STR_LIT>'],<EOL>"<STR_LIT>": b['<STR_LIT>']<EOL>}<EOL><DEDENT>else:<EOL><INDENT>height_meters_new = (width_meters / x_per_y_axes)<EOL>d_lat_new = ((b['<STR_LIT>'] - b['<STR_LIT>']) / height_meters) * height_meters_new<EOL>mean_lat = (b['<STR_LIT>'] + b['<STR_LIT>']) / <NUM_LIT><EOL>lat_min = mean_lat - d_lat_new / <NUM_LIT><EOL>lat_max = mean_lat + d_lat_new / <NUM_LIT><EOL>spatial_bounds = {<EOL>"<STR_LIT>": b['<STR_LIT>'],<EOL>"<STR_LIT>": b['<STR_LIT>'],<EOL>"<STR_LIT>": lat_min,<EOL>"<STR_LIT>": lat_max<EOL>}<EOL><DEDENT>return spatial_bounds<EOL>
Parameters
----------
bounds: dict
ax_width: float
ax_height: float

Returns
-------
spatial_bounds