def __contains__(self, g): geometry = ops.transform(self.__geo_transform__.rev, g) img_bounds = box(0, 0, *self.shape[2:0:-1]) return img_bounds.contains(geometry)