INVLIB::InvTrapMusic Class Reference

TRAP-MUSIC (Truncated RAP-MUSIC) source localization. More...

#include <inv_trap_music.h>

Public Member Functions
	InvTrapMusic (int iMaxSources=5, double dThreshold=0.85)
	Construct TRAP-MUSIC scanner.
QList< TrapMusicDipole >	compute (const Eigen::MatrixXd &matLeadField, const Eigen::MatrixXd &matData, const Eigen::MatrixXd &matSourcePos, int iNOrient=3) const
	Compute TRAP-MUSIC source localization.

Static Public Member Functions
static Eigen::VectorXd	scanCorrelations (const Eigen::MatrixXd &matLeadField, const Eigen::MatrixXd &matSignalSubspace, int iNOrient)
	Compute the MUSIC-type subspace correlation for all source locations.

Detailed Description

TRAP-MUSIC (Truncated RAP-MUSIC) source localization.

Extends RAP-MUSIC by adding a truncation step to the signal subspace at each iteration, improving robustness to correlated sources. Uses SVD of the measurement data to estimate the signal subspace, then iteratively scans the lead field for the best-matching dipole, projects it out, and truncates the subspace dimension.

Reference: Makela et al., NeuroImage 197, 616-626, 2019.

Definition at line 77 of file inv_trap_music.h.

Constructor & Destructor Documentation

InvTrapMusic()

InvTrapMusic::InvTrapMusic ( int iMaxSources = 5, double dThreshold = 0.85 )

explicit

Construct TRAP-MUSIC scanner.

Parameters

[in]	iMaxSources	Maximum number of sources to find (default 5).
[in]	dThreshold	Correlation threshold to stop scanning (default 0.85).

Definition at line 56 of file inv_trap_music.cpp.

Member Function Documentation

compute()

QList< TrapMusicDipole > InvTrapMusic::compute	(	const Eigen::MatrixXd &	matLeadField,
		const Eigen::MatrixXd &	matData,
		const Eigen::MatrixXd &	matSourcePos,
		int	iNOrient = 3 ) const

Compute TRAP-MUSIC source localization.

Parameters

[in]	matLeadField	Lead field matrix (n_channels × n_sources*n_orient).
[in]	matData	Measurement data (n_channels × n_times).
[in]	matSourcePos	Source positions (n_sources × 3).
[in]	iNOrient	Number of orientations per source (1=fixed, 3=free; default 3).

Returns

List of found dipoles, ordered by descending correlation.

Definition at line 64 of file inv_trap_music.cpp.

scanCorrelations()

VectorXd InvTrapMusic::scanCorrelations ( const Eigen::MatrixXd & matLeadField, const Eigen::MatrixXd & matSignalSubspace, int iNOrient )

static

Compute the MUSIC-type subspace correlation for all source locations.

Parameters

[in]	matLeadField	Lead field matrix (n_channels × n_sources*n_orient).
[in]	matSignalSubspace	Signal subspace (n_channels × n_signal_dims).
[in]	iNOrient	Number of orientations per source (1 or 3).

Returns

Correlation vector (n_sources).

Definition at line 158 of file inv_trap_music.cpp.

---

The documentation for this class was generated from the following files:

• src/libraries/inv/rap_music/inv_trap_music.h
• src/libraries/inv/rap_music/inv_trap_music.cpp