docs: def lon_lat_bins(bb, coord_bin_width):
      """
      Define bin edges for disaggregation histograms.
  
      Given bins data as provided by :func:`collect_bin_data`, this function
      finds edges of histograms, taking into account maximum and minimum values
      of magnitude, distance and coordinates as well as requested sizes/numbers
      of bins.
      """
      west, south, east, north = bb
      west = numpy.floor(west / coord_bin_width) * coord_bin_width
      east = numpy.ceil(east / coord_bin_width) * coord_bin_width
      lon_extent = get_longitudinal_extent(west, east)
      lon_bins, _, _ = npoints_between(
          west, 0, 0, east, 0, 0,
          numpy.round(lon_extent / coord_bin_width + 1))
      lat_bins = coord_bin_width * numpy.arange(
          int(numpy.floor(south / coord_bin_width)),
          int(numpy.ceil(north / coord_bin_width) + 1))
      return lon_bins, lat_bins