Hi!
Perhaps I should in more detail say, what I did:
b.file+orig .... bucket output from 3dDeconvolve
err.file+orig .... output od 3dDeconcolve using the "-error" option
3dTstat -MAD -prefix temp err.file+orig
3dcalc -datum short -a temp+orig -expr "astep(a,2.5)" -prefix temp2
3dcalc -datum short -a temp2+orig -expr "rect(a)" -prefix errmask
3drefit -fim errmask+orig
3drename b.file t.file
3dcalc -a errmask+orig -b t.file+orig -fscale -expr "a*b" -prefix b.file
rm -f temp+orig.* temp2+orig.* errmask+orig* t.file+orig.*
Please, let me know, you can see any problems in this procedure?
Thanks,
Lukas