Update hmrR_MotionCorrectPCArecurse.m (595784b0) · Commits · github_fork / Homer3

FuncRegistry/UserFunctions/hmrR_MotionCorrectPCArecurse.m

+11 −5

Original line number	Diff line number	Diff line
		% SYNTAX:
		% [data_dN, tInc, svs, nSV, tInc0] = hmrR_MotionCorrectPCArecurse(data_d, probe, mlActMan, tIncMan, tMotion, tMask, STDEVthresh, AMPthresh, nSV, maxIter, turnon)
		% [data_dN, tInc, svs, nSV, tInc0] = hmrR_MotionCorrectPCArecurse(data_d, probe, mlActMan, mlActAuto, tIncMan, tMotion, tMask, STDEVthresh, AMPthresh, nSV, maxIter, turnon)
		%
		%
		% UI NAME:
		@@ -19,6 +19,8 @@
		% SD.MeasListAct.
		% mlActMan: Cell array of vectors, one for each time base in data, specifying
		% active/inactive channels with 1 meaning active, 0 meaning inactive
		% mlActAuto: Cell array of vectors, one for each time base in data, specifying
		% active/inactive channels with 1 meaning active, 0 meaning inactive
		% tIncMan: Data that has been manually excluded. 0-excluded. 1-included.
		% Vector same length as d.
		% tMotion: Check for signal change indicative of a motion artifact over
		@@ -52,7 +54,7 @@
		% artifacts
		%
		% USAGE OPTIONS:
		% Motion_Correct_PCA_Recurse: [dod, tInc, svs, nSV, tInc0] = hmrR_MotionCorrectPCArecurse(dod, probe, mlActMan, tIncMan, tMotion, tMask, STDEVthresh, AMPthresh, nSV, maxIter, turnon)
		% Motion_Correct_PCA_Recurse: [dod, tInc, svs, nSV, tInc0] = hmrR_MotionCorrectPCArecurse(dod, probe, mlActMan, mlActAuto, tIncMan, tMotion, tMask, STDEVthresh, AMPthresh, nSV, maxIter, turnon)
		%
		% PARAMETERS:
		% tMotion: 0.5
		@@ -63,7 +65,7 @@
		% maxIter: 5
		% turnon: 1
		%
		function [data_dN, tInc, svs, nSV, tInc0] = hmrR_MotionCorrectPCArecurse(data_d, probe, mlActMan, tIncMan, tMotion, tMask, STDEVthresh, AMPthresh, nSV, maxIter, turnon)
		function [data_dN, tInc, svs, nSV, tInc0] = hmrR_MotionCorrectPCArecurse(data_d, probe, mlActMan, mlActAuto, tIncMan, tMotion, tMask, STDEVthresh, AMPthresh, nSV, maxIter, turnon)

		nBlks = length(data_d);
		for iBlk=1:nBlks
		@@ -80,6 +82,9 @@ tInc0 = tIncMan;
		if isempty(mlActMan)
		mlActMan = cell(nBlks,1);
		end
		if isempty(mlActAuto)
		mlActAuto = cell(nBlks,1);
		end
		if isempty(tIncMan)
		tIncMan = cell(nBlks,1);
		end
		@@ -94,7 +99,8 @@ tIncMerged = cell(1,1);
		for iBlk=1:nBlks
		data_dN(iBlk) = DataClass(data_d(iBlk));

		tInc(iBlk) = hmrR_MotionArtifact(data_d(iBlk), probe, mlActMan(iBlk), tIncMan(iBlk), tMotion, tMask, STDEVthresh, AMPthresh);
		tInc(iBlk) = hmrR_MotionArtifact(data_d(iBlk), probe, mlActMan(iBlk), mlActAuto(iBlk), tIncMan(iBlk), tMotion, tMask, STDEVthresh, AMPthresh);

		tInc0{iBlk} = tInc{iBlk};

		ii=0;
		@@ -102,7 +108,7 @@ for iBlk=1:nBlks
		ii=ii+1;
		tIncMerged{1} = min([tInc{iBlk}, tIncMan{iBlk}], [], 2);
		[data_dN(iBlk), svs_ii, nSV(iBlk)] = hmrR_MotionCorrectPCA(data_d(iBlk), mlActMan(iBlk), tIncMerged, nSV(iBlk));
		tInc(iBlk) = hmrR_MotionArtifact(data_dN(iBlk), probe, mlActMan(iBlk), tIncMan(iBlk), tMotion, tMask, STDEVthresh, AMPthresh);
		tInc(iBlk) = hmrR_MotionArtifact(data_dN(iBlk), probe, mlActMan(iBlk), mlActAuto(iBlk), tIncMan(iBlk), tMotion, tMask, STDEVthresh, AMPthresh);
		data_d(iBlk).Copy(data_dN(iBlk));
		svs{iBlk}(:,ii) = svs_ii{1};
		end

Admin message