flagmask

UNTITLED update flagging mask based on some interactive stuff
   Detailed explanation goes here

0001 function [ mask ] = flagmask( mask, LL, RR, QQ, UU )
0002 %UNTITLED update flagging mask based on some interactive stuff
0003 %   Detailed explanation goes here
0004 
0005 %apply initial mask
0006 mask=logical(mask);
0007 
0008 display('flagmask:: Initial flag channels'); find(mask) 
0009 
0010 LL(:,mask)=NaN;
0011 RR(:,mask)=NaN;
0012 QQ(:,mask)=NaN;
0013 UU(:,mask)=NaN;
0014 
0015 mask = rmsflag(LL, mask);
0016 mask = rmsflag(RR, mask);
0017 mask = rmsflag(QQ, mask);
0018 mask = rmsflag(UU, mask);
0019 
0020 end
0021 
0022 function mask = rmsflag(data, mask)
0023 happy='n';
0024 while(happy=='n')
0025 figure
0026 imagesc(data);
0027 figure
0028 plot(std(data));
0029 cut = input('flagmask::rmsflag:: rms level for cut? ');
0030 mask = logical(mask + std(data)>cut);
0031 close all;
0032 data(:,mask)=NaN;
0033 figure
0034 imagesc(data);
0035 figure
0036 plot(std(data));
0037 happy = input('flagmask::rmsflag:: Happy? (y/n) ', 's');
0038 end
0039 close all
0040 display('flagmask::rmsflag:: Flagging channels'); find(mask)
0041 end