def convert_bbox_to_radec_ranges(wcs, blc, trc):
    """
    Given BLC, TRC in coordinates of numpy array return RA, DEC ranges.
    :param wcs:
        Instance of ``astropy.wcs.wcs.WCS``.
    :return:
        RA_min, RA_max, DEC_min, DEC_max
    """
    blc_deg = wcs.all_pix2world(blc[0], blc[1], 1)
    trc_deg = wcs.all_pix2world(trc[0], trc[1], 1)

    ra_max = blc_deg[0]
    ra_min = trc_deg[0]
    dec_min = blc_deg[1]
    dec_max = trc_deg[1]

    return ra_min, ra_max, dec_min, dec_max