beaff02d91a80c2cf12e626d47b3f1678188472d,mne/source_space.py,,_get_volume_label_mask,#Any#Any#Any#,1852

Before Change

    vox_bool = mgz_data == vol_id

    // Get the 3 dimensional indices in voxel space
    vox_ijk = np.array(np.where(vox_bool)).T

    // Transform to MRI coordinates (where our surfaces live)
    _, vox_mri_t, _, _, _ = _read_mri_info(mri)
    rr_voi = apply_trans(vox_mri_t, vox_ijk)  // mri voxels -> MRI surface RAS
    // Filter out points too far from volume region voxels
    dists = _compute_nearest(rr_voi, rr, return_dists=True)[1]
    // Maximum distance from center of mass of a voxel to any of its corners
    maxdist = linalg.norm(vox_mri_t["trans"][:3, :3].sum(0) / 2.)
    return dists <= maxdist

After Change

    rr_vox = apply_trans(mri_vox_t, rr)
    good = (rr_vox >= -.5).all(-1)
    idx = np.empty(rr.shape[::-1], np.int64)
    for ii in range(3):
        good &= rr_vox[:, ii] < mgz_data.shape[ii] - 0.5
        idx[ii] = np.clip(np.round(rr_vox[:, ii]).astype(np.int64),
                          0, mgz_data.shape[ii] - 1)
    good &= mgz_data[tuple(idx)] == vol_id
    return good

Instances