def compute_image_median(image, quadrant=None): xmin, xmax, ymin, ymax = util.slice_indices_for_quadrant(quadrant) return np.median(image[ymin:ymax, xmin:xmax])
def compute_mad_std(image, quadrant=None): xmin, xmax, ymin, ymax = util.slice_indices_for_quadrant(quadrant) return mad_std(image[ymin:ymax, xmin:xmax])
def compute_std(image, quadrant=None): # non-robust standard deviation xmin, xmax, ymin, ymax = util.slice_indices_for_quadrant(quadrant) return mad_std(image[ymin:ymax, xmin:xmax])
def n_non_finite(image, quadrant=None): xmin, xmax, ymin, ymax = util.slice_indices_for_quadrant(quadrant) return np.sum(np.logical_not(np.isfinite(image[ymin:ymax, xmin:xmax])))