def rotate_voxel(xmin, ymin, xmax, ymax):
    """
    Given center position, rotate to the first quadrant
    
    Parameters
    ----------
        
    xmin: float
        low point X position, mm
    ymin: float
        low point Y position, mm
    xmax: float
        high point X position, mm
    ymax: float
        high point Y position, mm
        
    returns: floats
        properly rotated voxel in the first quadrant
    """
    
    xc = 0.5 * (xmin + xmax)
    yc = 0.5 * (ymin + ymax)
    
    if xc >= 0.0 and yc >= 0.0: # no rotation
        return (xmin, ymin, xmax, ymax)
    
    if xc < 0.0 and yc >= 0.0: # CW 90 rotation
        return (ymin, -xmax, ymax, -xmin)
        
    if xc < 0.0 and yc < 0.0: # CW 180 rotation
        return (-xmax, -ymax, -xmin, -ymin)
        
    # xc > 0.0 && yc < 0.0: # CW 270 rotation
    return (-ymax, xmin, -ymin, xmax)