code: def cut_edges(image, numPix):
      
      nx, ny = image.shape
      if nx < numPix or ny < numPix:
          print('WARNING: image can not be resized, in routine cut_edges.')
          return image
      if nx % 2 == 0 or ny % 2 == 0 or numPix % 2 == 0:
          #pass
          print("WARNING: image or cutout side are even number. This routine only works for odd numbers %s %s %s"
                           % (nx, ny, numPix))
      cx = int((nx-1)/2)
      cy = int((ny-1)/2)
      d = int((numPix-1)/2)
      if nx % 2 == 0:
          cx += 1
      if ny % 2 == 0:
          cy += 1
      resized = image[cx-d:cx+d+1, cy-d:cy+d+1]
      return copy.deepcopy(resized)
docs: cuts out the edges of a 2d image and returns re-sized image to numPix
      center is well defined for odd pixel sizes.
      :param image: 2d numpy array
      :param numPix: square size of cut out image
      :return: cutout image with size numPix