title: Histogram equalization with SimpleITK
body: <p>On this <a href="https://stackoverflow.com/a/59950816/68571">SO answer</a> suggest me this code:</p>
  
  <pre><code>import SimpleITK as sitk
  import numpy as np
  
  # Create a noise Gaussian blob test image
  img = sitk.GaussianSource(sitk.sitkFloat32, size=[240,240,48], mean=[120,120,24])
  img = img + sitk.AdditiveGaussianNoise(img,10)
  
  # Create a ramp image of the same size
  h = np.arange(0.0, 255,1.0666666666, dtype='f4')
  h2 = np.reshape(np.repeat(h, 240*48), (48,240,240))
  himg = sitk.GetImageFromArray(h2)
  print(himg.GetSize())
  
  # Match the histogram of the Gaussian image with the ramp
  result=sitk.HistogramMatching(img, himg)
  
  # Display the 3d image
  import itkwidgets
  itkwidgets.view(result)
  </code></pre>
  
  <p>Why do I need two images to do Histogram equalization?</p>
  
  <p>Because I want to do Histogram Equalization, and this is Histogram Matching. In this <a href="https://automaticaddison.com/difference-between-histogram-equalization-and-histogram-matching/" rel="nofollow noreferrer">article</a> explain the different.</p>