code: def check_grid_aligned(vol, img, offset):
    
    shape = Vec(*img.shape)[:3]
    offset = Vec(*offset)[:3]
    bounds = Bbox( offset, shape + offset)
    alignment_check = bounds.expand_to_chunk_size(vol.underlying, vol.voxel_offset)
    alignment_check = Bbox.clamp(alignment_check, vol.bounds)
    is_aligned = np.all(alignment_check.minpt == bounds.minpt) and np.all(alignment_check.maxpt == bounds.maxpt)
    return (is_aligned, bounds, alignment_check)
docs: Returns (is_aligned, img bounds Bbox, nearest bbox inflated to grid aligned)