https://surfer.nmr.mgh.harvard.edu/fswiki/AutomatedDefacingTools

[osotolongo@detritus facehbi]$ /nas/usr/local/bin/mri_deface_linux bids/sub-0001/anat/sub-0001_T1w.nii.gz /nas/usr/local/share/mri_deface/talairach_mixed_with_skull.gca /nas/usr/local/share/mri_deface/face.gca deface/sub-0001_T1w_defaced.nii.gz
logging results to sub-0001_T1w_defaced.nii.log
reading '/nas/usr/local/share/mri_deface/talairach_mixed_with_skull.gca'...
reading 'bids/sub-0001/anat/sub-0001_T1w.nii.gz'...
changing type of input volume to 8 bits/voxel...
MRIchangeType: Building histogram
bounding unknown intensity as < 10.3 or > 1461.6
total sample mean = 111.8 (0 zeros)
spacing=8, using 2055 sample points, tol=1.00e-03...
resetting wm mean[0]: 143 --> 144
resetting gm mean[0]: 89 --> 89
input volume #1 is the most T1-like
using real data threshold=7.0
skull bounding box = (48, 17, 11) --> (211, 255, 216)
using (102, 96, 114) as brain centroid...
mean wm in atlas = 144, using box (82,67,89) --> (122, 125,139) to find MRI wm
before smoothing, mri peak at 133
robust fit to distribution - 134 +- 8.9
after smoothing, mri peak at 134, scaling input intensities by 1.075
scaling channel 0 by 1.07463
initial log_p = -4.7
 
after initial translation: (-5.7, 19.3, -5.7): log p = -4.3
 
scale 1.000: max=-3.9, old_max =-4.3 (thresh=-4.2)
 
scale 1.000: max=-3.8, old_max =-3.9 (thresh=-3.9)
 
scale 1.000: max=-3.8, old_max =-3.8 (thresh=-3.8)
reducing scale to 0.2500
....
....
....
anonymizing volume...
using wm = 107.0, gm = 93.4
wm covar:
 502.39935;
gm covar:
 659.34979;
.........................
1058443 face voxels erased, 250588 ambiguous voxels retained
resampling to original coordinate system...
writing anonymized volume to deface/sub-0001_T1w_defaced.nii.gz...
deidentification took 2 minutes and 54 seconds.