def mk_crop_np(product, x_start, y_start, size=256):
    """ crop a squared region size^2
        provide upper-left corner with (x_start, y_start)
    """
    return product[y_start:y_start+size, x_start:x_start+size]