Skip to main content

heidi_seed_problem

qsm_core::pipeline::inversion

Function heidi_seed_problem 

Source 
pub fn heidi_seed_problem(config: &InversionConfig) -> Option<String>
Expand description

Why HEIDI’s cone fill cannot follow config’s inversion, if it cannot (None when HEIDI is off or the seed is allowed). The hard errors of run_heidi_cone_fill, for hosts to check at plan time:

  • the method has no local field to fill the cone from (the total-field and phase methods);
  • TKD’s or TSVD’s threshold is above HEIDI’s cone threshold: the coefficients they truncated would be kept as data (upstream: “DipoleFilter parameter must not be larger that Cone threshold”, computesusceptibility_heidi.m L922-923).

The softer advice on seeds (LSQR is the reference; regularised seeds are discouraged) is in the heidi.* rules (super::rules).