Choosing dimension for filter

20 visualizaciones (últimos 30 días)
Aboltabol
Aboltabol el 28 de Mzo. de 2023
Comentada: Star Strider el 28 de Mzo. de 2023
I have voltage recordings from x sensors for y time-frames, each spanning z seconds. This is stored in a x*y*z matrix. Now, I designed an iirnotch to remove the line noise (60 Hz). When I apply the filter to the matrix using filtfilt, am I correct in assuming that the filter will run along the x dimension, and is this choice of dimension appropriate?

Respuestas (1)

Star Strider
Star Strider el 28 de Mzo. de 2023
I believe so.
From the filtfilt documentation for x:
The function operates along the first array dimension of x unless x is a row vector. If x is a row vector, then the function operates along the second dimension.
So I believe it will operate along each column (down the rows) of the matrix and across all dimensions. (I have only filtered 2D matrices with it, not 3D.)
.
  2 comentarios
Aboltabol
Aboltabol el 28 de Mzo. de 2023
Editada: Aboltabol el 28 de Mzo. de 2023
My matrix has sensors as rows, events as columns, and time as pages. Now, ought not the filter run per sensor per event along the time dimension i.e. the third dimension (page)?
Star Strider
Star Strider el 28 de Mzo. de 2023
Yes, it should, however that information did not originally appear.
Use the permute function to put them in the order you want them, for example —
A = randn(12, 5, 30); % Sensor x Event x Time
Ap = permute(A, [3 1 2]);
PermutedSize = size(Ap) % Time x Sensor x Event
PermutedSize = 1×3
30 12 5
That appears to me to be appropriate for the filtering you want to do. If it is not, changing it is straightforward.
.

Iniciar sesión para comentar.

Categorías

Más información sobre EEG/MEG/ECoG en Help Center y File Exchange.

Etiquetas

Productos


Versión

R2022b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by