- T1w and T2*w thresholds are chosen with different RDs
- Added script for refining T2*w thresholds with fminsearch()
- Modularized pSegment_us and validate()