Get S_decomp.
(
exp, all_coils, trans, coil_scale, cal, ignore_ref, grad_picks, mag_picks, mag_scale
)
| 1369 | |
| 1370 | |
| 1371 | def _get_s_decomp( |
| 1372 | exp, all_coils, trans, coil_scale, cal, ignore_ref, grad_picks, mag_picks, mag_scale |
| 1373 | ): |
| 1374 | """Get S_decomp.""" |
| 1375 | S_decomp = _trans_sss_basis(exp, all_coils, trans, coil_scale) |
| 1376 | if cal is not None: |
| 1377 | # Compute point-like mags to incorporate gradiometer imbalance |
| 1378 | grad_cals = _sss_basis_point(exp, trans, cal, ignore_ref, mag_scale) |
| 1379 | # Add point like magnetometer data to bases. |
| 1380 | if len(grad_picks) > 0: |
| 1381 | S_decomp[grad_picks, :] += grad_cals |
| 1382 | # Scale magnetometers by calibration coefficient |
| 1383 | if len(mag_picks) > 0: |
| 1384 | S_decomp[mag_picks, :] /= cal["mag_cals"] |
| 1385 | # We need to be careful about KIT gradiometers |
| 1386 | return S_decomp |
| 1387 | |
| 1388 | |
| 1389 | @verbose |
no test coverage detected