code: def plot_img_and_hist(image, axes, bins=256):
      """Plot an image along with its histogram and cumulative histogram."""
      ax_img, ax_hist = axes
      ax_cdf = ax_hist.twinx()
  
      # Display image
      ax_img.imshow(image, cmap=plt.cm.gray)
      ax_img.set_axis_off()
  
      # Display histogram
      ax_hist.hist(image.ravel(), bins=bins)
      ax_hist.ticklabel_format(
          axis="y", style="scientific", scilimits=(0, 0)
      )
      ax_hist.set_xlabel("Pixel intensity")
  
      xmin, xmax = dtype_range[image.dtype.type]
      ax_hist.set_xlim(xmin, xmax)
  
      # Display cumulative distribution
      img_cdf, bins = exposure.cumulative_distribution(image, bins)
      ax_cdf.plot(bins, img_cdf, "r")
  
      return ax_img, ax_hist, ax_cdf
  
file_path: scikit-image/doc/examples/color_exposure/plot_local_equalize.py
imports: import numpy
  import matplotlib
  import matplotlib.pyplot
  from skimage import data
  from skimage.util.dtype import dtype_range
  from skimage.util import img_as_ubyte
  from skimage import exposure
  from skimage.morphology import disk
  from skimage.morphology import ball
  from skimage.filters import rank
label_desc: Use image node_id 1 for calling a global function with example usage: plot_img_and_hist(image, axes, bins) and returns: ax_img,  ax_hist,  ax_cdf