def get_r_s_nz(self):
        """ Return (reflectance, shading), with just the nonmasked entries """
        s_nz = self.input.image_gray_nz / self.intensities[self.labels_nz]
        r_nz = self.input.image_rgb_nz / np.clip(s_nz, 1e-4, 1e5)[:, np.newaxis]
        assert s_nz.ndim == 1 and r_nz.ndim == 2 and r_nz.shape[1] == 3
        return r_nz, s_nz
dataset/ETHPy150Open seanbell/intrinsic/bell2014/decomposition.py/IntrinsicDecomposition.get_r_s_nz/original