def METHOD_NAME():
    region = CircleAnnulusPixelRegion(PixCoord(0, 0), 10, 20)
    data = np.ones((101, 101)).astype(int)
    cutout = region.to_mask().multiply(data, fill_value=np.nan)
    xypos = ((20, 20), (5, 5), (5, 35), (35, 5), (35, 35))
    for x, y in xypos:
        assert np.isnan(cutout[y, x])