code: def __ms_npenalty_fcn(self, axis, mask, orig_shape):
          
          maskz = zoom_to_shape(mask, orig_shape)
  
          maskz_new = np.zeros(orig_shape, dtype=np.int16)
          maskz_new[maskz == 0] = self._msgc_npenalty_table[0, axis]
          maskz_new[maskz == 1] = self._msgc_npenalty_table[1, axis]
          # import sed3
          # ed = sed3.sed3(maskz_new)
          # import ipdb; ipdb.set_trace() #  noqa BREAKPOINT
  
          return maskz_new
docs: :param axis: direction of edge
          :param mask: 3d ndarray with ones where is fine resolution
  
          Neighboorhood penalty between small pixels should be smaller then in
          bigger tiles. This is the way how to set it.